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THE PROBLEM OF AN INFINITE PLATE UNDER AN 
INCLINED LOADING, WITH TABLES OF THE 
INTEGRALS OF Ai(+x) AND Bi(+x) 

By M. ROTHMAN (The Northern Polytechnic, Holloway, N.7) 
Received 4 September 1951; revised 8 May 1952] 


SUMMARY 


The present paper deals with the problem of an infinite plate uniformly loaded at 
right angles to the infinite edges which are freely pivoted, the other two being free. 
The solution is found to depend ultimately on tables of the functions { Ai( +2) dx 


and | Bi vr) dx, all of which have been tabulated over the range |. 10 except 


) Bil ) dx which is given for the range a 2. 


Ta of the functions f Gi ix and { Hi(—-) dz are also given to a 2, these 


being necessary for the complete solution of the problem. 


1. Introduction 

THE action of wind or rain on a large door introduces the problem of an 
inclined loading; this may be assumed constant over the surface of the 
door, taken to be a rectangular plate, which may be subject to a variety 
of boundary conditions but which is commonly freely pivoted along two 
opposite edges and free along the other two as in a sliding door. 

This paper considers the simplest case in which the load is assumed to 
be at right angles to the freely pivoted edges, which are in turn assumed 
infinitely long, and at a constant angle to the plane of the plate. This is 
a suitable approximation for a long door. 

Other possible cases are those of a finite plate subject to the same loading, 
and both cases subject to a skew loading. 


2. Infinite plate under constant inclined loading 
We assume the plate to be freely pivoted along the infinite edges and to 
be free along the other two, so that the problem becomes two-dimensional. 
Also if the loading be p per unit width and be inclined at an angle @ to the 
plane of the plate, then, taking axes as in Fig. 1, we are led to the equation 
dM Y _dy 
dx ‘ dx’ 
where X, Y are the horizontal and vertical components of force and M the 
bending moment at the point (v,Yy). 
If D is the plate modulus, this equation becomes 
d3y , dy 
da? : dx - 


{Quart. Journ. Mech. and Applied Math., Vol. VII, Pt. 1 (1954)] 
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and as X = pxcos@+S and Y = pxrsin@—}plsin@, it can be written as Thu 

dy | ely 

dé ~ ae 

where € = X, and A, B, C are constants given by the relations 
A (Dp? cos*6)-}, B = tan 6(Dp* cos?@)~", 


C = (C,—C, tan @)(Dp* cos*0)-1, C,= 8S, and C, pl sin 6, St 


BE+C, (2.1) 


S being the horizontal reaction at the hinge O. 








ipl sin@ ipl sin@ 
« V+dY , 
aif i.e. 
sd . rad K+6X 
Kp S, y —_ A +02 7” 
/ M+dM = 
depe: 
y 
Y 
Fic. 1. 3. T 
The complementary function of (2.1) is given by ™~F¢,; 
W = Ea, (3 Ate!) + BI_, (3A 4S?) I, W = dy/dé, on ta 
or W aAi(—£)-+-bBi(—€&), «, 8, a, b being constants, and the particular Th 
integrals corresponding to Bé and C by reduc 
BEIA 1 CO < 9 Avr! £9r48)(9¢-4-2)! Mille: 
BE/A anc > [| 3A)r! (3r+3)!] of Bi 
respectively, so that the complete solution of equation (2.1) is advic 
. , Thi 
y x | €40,(2 AE?) dé 4 B | Et J_,(3A*E#) dé + far as 
xr of 1 t 
| SAasLOy ¢ ry! £3r4+3/(Q2y_1 Q\I1_L FW 99 ; 
BE|A+C ¥ [(—BA yr! &749/(31 3)'+F. (2.28 pe, 
The last but one term on the right-hand side of (2.2) can also be written the in 
¢ 301 
{; p~ 7 
al Gi(t) dt, = 
7. | Ait 
0 . 
Gi(t) being the function defined by Scorer (1), which can be shown to satisfy 
the equations 
; : 
Gi(t) = 4Bi(t)+Ai(t) | Bi(u) du — Bi(t) | Ai(w) du, 
0 0 = 
t t V3 
Gi'(t) = 4 Bi'(t)+ Ai’‘(t) | Bi(w) du — Bi’(t) | Ai(w) du, — 
0 0 
, , ' The 
Gi’ (t)—t Gi(t) = —1)z. 
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Thus (2.2) can be written 


‘ti | ( Ms 
y a | Ai(t)dt —b | Bi(t) dt Gi(t) dt 4 4 +F. (2.3) 


Substituting the boundary conditions 


7/2 
a-y 
y 0, QO when2z 0 and J, 
dx? 
dy " y . 
i.e. y= 0, 73 0 when € = C, and plcosé+C, 
dr2 


we then have four equations for the four constants a, b, F, and C (which 
depends on the unknown S 


3. Tables of the integrals of Ai(+-x) and Bi(+-2) 


~ ~#HXis seen from (2.3) that the complete solution of the problem is dependent 


on tables of | Ai(a) dx and | Bi(x) dz. 

These integrals were computed at intervals of 0-1 using the method of 
reduced derivatives described in (2), and the author is grateful to Dr. J.C. P. 
Miller for lending him the manuscript of his tables of reduced derivatives 
of Bi(w), his unpublished table of same for Ai(x), and for some very helpful 
advice and criticism. 


This integral of Ai(x) was taken to |x| < 10, but is given here only as 
far as x = 7-5 in the positive sense, having reached its asymptotic value 
of } to 7 decimal places. The integral of Bi(a) was also taken to x 10 
but only to 2 2 as it tends to infinity in the positive sense and it is 


the integral of Bi(—x) which is of primary importance. 
Both integrals were also computed from their ascending power series 
10 
2 x a 
\i(+2) da ; Xu Se ee 
3°13) 2.3.4 2.3.5.6.7-- 2.3.5.6.8.9.10 





) 


ior |a 2. at intervals of 0-1. 


The results here given are to 7 decimals, though the tables have been 


computed to 8 decimals throughout and for intervals of 0-01 to x = 2 


ay 
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It might be noted that 


( Ai(t) dt = a} Ai’ (%)Gi(a) —Ai(a)Gi' (x)! f 1 


0 


{ Ai(x) Hi’ (x) — Ai’ (x) Hi(x)} —, 


( Bi(t) dt = w{Bi'(x)Gi(x) — Bi(x)Gi'(x)} 77{ Bi(a) Hi’ (a) — Bi’ (a)Hi(a)} 

0 (3.1) 
so that the tables of Gi(a), Gi'(a), or Hi(x), Hi’(a2) could have been com- 
puted as alternatives. In general this would be preferable and, in fact, 
tables of Gi(a) and Hi(z), but not their derivatives, are available in (1), as 
these are better behaved functions than the ones computed; but the 
solution of the problem under consideration would necessitate the applica- 
tion of formulae (3.1) so frequently that it was thought advisable to keep 
to the integrals of Ai(--x) and Bi(-+_-2). 

It is hoped to publish tables of integrals of Gi(x) and Hi(—2) later. 

[ also wish to express my thanks to my former colleague Mr. R. D. Lord 
for his helpful criticism and suggestions. 


,EFERENCES 
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THE EFFECT OF HEAT TRANSFER ON THE 
SEPARATION OF A COMPRESSIBLE LAMINAR 
BOUNDARY LAYER 


By C. R. ILLINGWORTH 
(Department of Mathematics, The University, Manchester) 
[Received 13 March 1952; revised 19 September 1952] 


SUMMARY 

The paper treats problems of laminar boundary layers with non-uniform density, 
due either to large temperature differences or to large Mach numbers. An approxi- 
mate method is developed for cases when the stream velocity and the wall tempera- 
ture are non-uniform. In the examples to which the method is applied, only uniform 
wall temperatures are considered. Attention is concentrated on the effect of a given 
pressure gradient in producing thinning, thickening, or separation of the boundary 
layer at a wall at various uniform temperatures. 

For zero Mach number and a small temperature difference, and various given 
main-stream velocity distributions, the calculations show fairly satisfactory agree- 
ment with known exact results. Under other conditions the effect of the pressure 
gradient is found to be enhanced or reduced according as the ratio of stream 
temperature to wall temperature is less or greater than 1. Thus, in uniformly retarded 
flow the relative fall in velocity before separation varies as the square root of 
this ratio, showing little independent variation with Mach number. This result is 
checked against previous calculations at high Mach numbers for zero heat transfer. 
The general conclusion mentioned above is explained by the fact that pressure 
gradient produces thinning, thickening, or separation principally by its action on the 
fluid very near the wall, and so the effects are enhanced if and only if this fluid is 
hotter (and so lighter) than that in the main stream. 


1. Introduction 

THE study of the laminar boundary layer of a gas in steady two-dimensional 
flow past a rigid wall is now well developed in the case when the main stream 
is uniform. The work of Crocco (1) deals with the boundary-layer flow when 
the temperature of the wall is uniform, both with and without heat transfer 
between the wall and the gas, and the work of Chapman and Rubesin (2) is 
concerned with the flow when the wall temperature is non-uniform, in which 
case there is heat transfer. 

However, when the main-stream velocity is non-uniform, the calculation 
of the boundary-layer flow is in general a formidable task, which may require 
the employment of an electronic calculating machine along the lines initiated 
by Cope and Hartree (3). This task is very much simplified under three 
conditions, viz. (i) that the Prandtl number of the gas is 1, (ii) that its 
viscosity is proportional to the absolute temperature, and (iii) that there 


[Quart. Journ. Mech. and Applied Math., Vol. VII, Pt. 1 (1954)] 
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is no heat interchange between the wall and the gas. The first and third 
of these conditions together imply that the wall temperature is uniform. 
The simplicity effected under these conditions lies in the fact that the 
compressible boundary-layer flow can be easily expressed in terms of a 
related incompressible flow. 

The object of the present paper is to examine compressible boundary 
layers with non-uniform main streams and non-uniform wall temperatures. 
The condition of zero heat transfer at the wall is of course discarded, but 
only at the cost of introducing approximations in the boundary-layer 
equations: the other two conditions concerning the Prandtl number and 
the viscosity : 


re retained. These approximations are made by means of an 
artifice used by Lighthill (4) in his discussion of heat transfer in boundary 
layers. After solutions of the approximate equations have been given in 
section 2 for cases when the main-stream velocity and the wall-tempera- 
ture distributions can be expressed as power series, some particular cases 
are discussed. These all concern uniform wall temperature for which the 
computation is relatively simple. In section 4 the boundary layer of a 
blunt-nosed symmetrical cylinder is considered, and in section 5 the 
boundary layer corresponding to a linearly decreasing main-stream velocity 
is discussed in some detail. This is a suitable test case as regards separation 
of the boundary layer. For in this case the point of separation for incom- 
pressible flow has been calculated with 1 per cent. accuracy by Howarth (5), 
and Stewartson (6) has solved the same problem in compressible flow, under 
the three simplifying conditions already referred to. It has been found 
convenient, therefore, to use this particular flow in looking for alterations 
in the position of separation when the wall temperature is changed from the 


uniform value corresponding to no heat transfer. 
2. The approximate equations 


The exact equations of steady two-dimensional boundary-layer flow when 
the Prandtl number is 1 are 


pu)/ex- O(pv)/cy 0), (1) 
o(u eu/ex+-v eu/ey) p, U, du, /dx+-o(p Ou ey)/ey, (2) 
p(u 0t/Ox+-v Ci/Cy)+u p, U, du,/dx O(p ot/ey) Cy +-p(Cu/ey)?, (3) 


where x, y are orthogonal curvilinear coordinates measured along and 
perpendicular to the wall, u, v are the corresponding velocity components, 
p, i are the density and specific enthalpy, p is the viscosity, and the suffix 1 
refers to the main stream, where conditions depend only on x. From (2) and 
(3) it follows that 


p(u ch/dx+-v oh/ey) O(u ch/oy)/ey, (4) 
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where h is the specific total energy i+4u?. In the main stream h, is a 
constant (equal to the stagnation enthalpy) by Bernoulli's equation, and 
h = constant = h, is the particular integral of (4) appropriate for the case 
of no heat transfer at the wall. 

In terms of von Mises’s coordinates 2. us. where us is the stream function 


defined by pu éys/Cy, pv = —Ors/ex, (2) and (4) become 
pu ou/ex p, Uy du,/da-+- pu O(upu Cu/ ox) / eps, (5) 
oh/ex O(upu Ch/ Ors) /Oxb. (6) 


It is now convenient to introduce non-dimensional dependent variables f, ¢ 
instead of u, h by means of the equation 


u i,(b—f), _ (7) 


where b = u?/i, (y—1)M}, M, being the Mach number of the main stream, 


and the equation h = h,(1—g). (8) 
With the help of Bernoulli’s equation in the form di,/dx+-u, du,/dx = 0, 
(5) becomes 

u O(up Of Gb) / xb — Of Cx +db/dx-- (di,/i, dx)| 2(p,/p)+b—f | 0. (9) 
Now the pressure p in the boundary layer at the station x is the same as 
the main-stream pressure p,(x), and therefore, with attention restricted to 
a perfect gas with constant specific heats for which p = (y—1)pi/y, it 
follows that 


P1/P t/ty. (10) 
From (7) and (8), i/i, = (1+4b)(1—g)—4(b-f), (11) 
and, since 7,(1+4b) = A, = constant, 
di,/i,dx —(db/dx)/(2+-b). (12) 
Equations (9)—(12) give 
u O(up Of es) es —Of/ex+(db/dx)g = 0, (13) 
and (6) and (8) give 
O(upu 6g/exbs)/Os—Og/ex = 0. (14) 


Equations (13) and (14) are our final forms for the exact boundary-layer 
equations when the Prandtl number is 1, and they are completed by the 
addition of boundary conditions. The latter are derived from the facts that 
u and h must be equal to uw, and h, at the edge of the boundary layer, and 
that at the wall (quantities relating to which will be characterized by the 
suffix 0) uw = 0 and either the enthalpy distribution i,(7) or the rate q,(x) 
of heat transfer per unit area from the wall to the gas is prescribed. Hence 
the boundary conditions are: 

f=g = 0 at x = 0, for all values of %, and at % = oo for all positive 
values of x; 
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f = b, and either g l—ij/h, or —poi/ey = ppuh, ég/eb = qy at 
0 for all positive values of x. 


j 


Following Lighthill (4), approximations are made in (13) and (14) by 
replacing the factor wu which occurs in each equation by an expression 
obtained as follows. By identifying the shearing stress at the point (x, #) 


with the stress 7)(x) at the wall at the same station, the equation 


LOu/cy = ppu 0u/oys T p(X) (15) 
is obtained. Now the remaining condition to be imposed, that viscosity is 
proportional to absolute temperature, implies that 

| [Lp/ Py Py pr) py ry 3 (16) 
Hence (15) can be integrated to give 
27,(x)y |4 ie 
uw = | 2eole h (17) 
by(%) py (x) 
This approximate expression for u is substituted in (13) and (14). It is most 
accurate near the wall (according to (15) it gives the slope of the velocity 


profile correctly at yy = 0), but it overestimates u in the outer parts of the 


: boundary layer. However, the inaccuracy of u becomes less serious the 
greater the distance from the wall, because the cofactors of uw in (13) and 
(14) tend to zero as { — co. In making this substitution it is timely to recall 
Lighthill’s quantitative statement concerning the errors which the substi- 

tution involves at low speeds. For main streams of the form u, = cz”, 
) 


Lighthill shows that the shearing stress at the wall is underestimated by 
6 per cent. when the index m = 0 and by nearly 8 per cent. for all yalues 


) of m greater than 4, and that the heat transfer to the wall, when it is at 
a uniform temperature, is overestimated by 3 per cent. when m = 0, by 
, 18} per cent. when m 1, and by 24 per cent. as m > o0. 
When the substitution is made, (13) and (14) become, with the help 
of (16) e — 
}) ust 0*f/ us? — of /et+-(db/dt)g = 0, (18) 
2 ¢ (ys? Og Cys) Cds og ot Q, (19) 
ie 
t where ¢ = | | 2u,(&)p,(€)79(€) |! d€ is the new independent variable. The 
0 
d Heaviside transforms of these equations are 
L€ = - 
d? z db 
r) ed nh Og mt, (20) 
difs* dt 
“ 
d dg 
yi |) — ng = 0, (21) 
re aI dis Pg 


since f = g 0 at t 0, where bars denote transforms and p denotes 
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the operator @/ét.+ These equations are formally the same as equations 
(66) and (21) of Lighthill’s paper, save for the extra term (db/dt)g in (20) 
arising from the Mach-number gradient in the main stream. The solution 
of (21) satisfying the conditions g = 0 at 4 = 0 and g = 1—i,(t)/h, at 
ub 0 is , 

- 7 9. 4738\1—EZ wee 29(t) 99 
g ——— -(§ pty) Kh, (4p??)| 1 ———}, (22) 

\—#s)) h, 
where Kh,,(z) is the Bessel function [J_,,(z)—J,,(z)|/sinnz. If qo. and not 7%, 
is specified, (22) will still stand as the operational form of 7 with 7,(t) to be 
determined from the relation 
do [upuh, og/eb|y-o [ (241 pr to )'hy Og/eyf]y,-0, 

which gives 


t 
(_2_\ige_ (2) _ _@*/1—isl+oyhy_ 1 aileya 
My Pit) hy \edbt) yo (4)! t h,. 
0 


G—€P 

(23) 
corresponding to Lighthill’s equation (26). The problem is clearly more 
straightforward if the wall temperature is specified, and in what follows 
only this alternative will be considered. 

At this stage it is also convenient to introduce restrictions on the forms 
of the main-stream velocity and wall-temperature distributions, namely 
that 0 o 
b(t) = > B,t™, = 14, (t)/hy = ¥ a, t™, (24) 

=0 


n=0 n= 
where m is a positive number and where the «’s and f’s are constants to 
be determined later. Then (22) becomes 


gJ — (3)! S (mn)! x, pm" 4AKh, (4 pips), 
( 3): Rogen 


and so, according to Jeffreys and Jeffreys (Methods of Mathematical Physics, 
2nd edition, p. 613), 


g = 2% 9n()t™, (25) 
n=0 
where » = 4:3/9¢ and 
mn)! : eae l Bi " ; 
J,() A sro 1U(mn-+-&, 3, n) sy | en 8(E—n)mmE-m"-4 dé, (26) 
(—#)! Let. 
: ” 
in which 
(—y)! xz y—2)! ; 
U(a,y,n) = — Fy (a y; 7) ay <h YF (l+a—y, 2—y, 7), 
( ): (aE)! 


+ Since the pressure will not be required again, there should be no confusion in this new 
use of the letter p. 
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in the usual notation for confluent hypergeometric functions. The function 


fis given by . 
un len 
f > h,(n)tm", (27) 
n 0 
where h,(0 3. h(a 0, and 
n—1 
yh; h mnh',+m > (n—r)a,B,_,g, 0 (28) 
r=0 
where dashes denote differentiation with respect to 7), as can easily be 





verified by substit m in (18). The solution of (28) for the given 
boundary conditions 
1 7 
h -, 3m > (n—r)a,B,, fat) | Fy (mn-+- 1, 3. €)9,(€) dé 
0 t 
0 
sn 5 es tere 29 
ey Fi (mn 1, 3. 7) E—Sesf fg) ae}. (29) 
7 
where . 
tbe l l , Jf mné—mn—s ¢ 
f -e-1U)(mn+4, 4, ») — e$(E—y)mng—mn—t dé, 
) ( eS 
7] 
The solution is completed by expressing x in terms of t, because it is then 
straightforward matter to find the «’s and f’s from the data, which give 
ind s functions of To do this, 7, must be expressed in terms of t 
by means of the equation 
T (pu CU) ¢ Spy Py (Of Cys) 0 
x 
5) 91 ° ; ‘ : a2,’ nn 
zZ~ 309 fly Py tz ¢ > { lim n h,, \t" ; 
n=0 ‘n--0 
Now, fo1 mal 
(mn)! (—8)! 
f | ania a Y O(n) 
1\) V7] ] 
( x). (mn 5): 
n 1 
] Qa , Q 
S tl ll / kB, pa 
r=0 
. (mn)! ( z)! 
wnere A we it ¥ ss ’ 
(mn 5) ( L) 
f » £2 ss . 
und m(n—r) 3 eff (E)g,(E) dé 
0 
since is a solution of “fe (y 1) fn mnf » 0 
AGL YT) i +) te ’ be 
| S| T[eteffn(élg(€) dé 
mn . dé 
0 
@ 
Vt / tp . ti’ € 
k,— [ Betfalear(é) ae. (30) 
n 
0 











( 


14 Cc. R. ILLINGWORTH 
It is convenient at this juncture to evaluate k,,, which will be required F 
later. Since 
7) mn F Ps 
frl€) ae iy] | e-“(u &)™n ly mn—-k dy. 
é 
 eaote (£\q'(¢ eye egy ae 
Ete f, (E)go(E) dé = ay | § fn(€) dé 
0 aor n 
mn . i “ /mn—1 
—_—_— | é-4 | eu S (—£)'u-*- du} dé 
1)!(—2)! — r 
3 g)+ : A 
0 é 
x @ 
mn = (mn— 1 ] 
——. =< > ( \ 1)’ —. [ uu-i du 
Penal or r+. 
0 
after integration by parts, and this is equal to 
1 
ws ! 
mn | & 4(] -£)mn 1 dé _ (mn): 
gy ‘ : (mn—14)! 
0 
mn— 2)! Ly! 
Hence k= bk y — ee: 3): | 3): 
n n 1\! 2\1 
(mn—-+)! (—)! 
. 31) 
(mn)! (—2)! 
where k, moh ay : 
(mn— 4)! (—4)! 
now be expressed as 
~ B. —, 


The skin friction may 


To(t) 
8 fmn ( 
of sf 


7 


tu. petal (1 LJ > 


» 


at at oe 
t. 


Us Pais Pet 
It follows that 


t 
x = | [2m (€)pr(E)ro(€)]* 4 


4 
$ 


0 
t 
34(24, pit) f e(144 


0 





Dn, 


o n—1 
1)i(y—1) % . . kL 
™ > (k,, B,.— 2 a 

n=0 r=0 


“y 


s—1){(v—-1) 
_ (1+ 1h) (2y-Di(y 


where the suffix s refers to stagnation conditions in the main stream, so that 
—1) 


2 
> B,, gmn) ay 1)/Ay—) y 
n=0 


x[E (ba Bu—"S ne teBn-n)emn] de 


(32) 


(33) 


d 
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Finally, from (23), the rate of heat transfer per unit area from the wall is 


x 
a [9s38(1)!'|-1y 9 Git i(] | 1 a B pms (By—)N/Ay—-D y 
} s é - — a” 
n 





x 
; + ~- (mn)! (—4)! ’ 
> (kt, B,— > k,,o,B, jen | S —_——- 3) mo. (54) 
a 7 ‘ 1\! U 
a4 +=(mn—-)! 
%= 0 i 
The boundary-layer flow. has now been expressed in terms of the para- 
meter ¢. The velocity is given by 
i(b f) 
1+1b 
1l+if 
the enthalpy by i i ae 
P. j l+1ib ° 
wad the value of x by (33). The y-coordinate is recovered from 
* duds 3 P 4 +o f— (4 4-36 
y u ey 1h\y 1)/2Ay—1)pi | J +: 29)9 | 7 
pu 2 Pots s p 7) (b | yi 


0 0 





Finally. the shearing stress and heat transfer at the wall are given by (32) 
and (34). 

So far the number m has not been specified. It is clear, from (33), that 
if 8. is the first non-vanishing f, then 2 o t!-?"”"" near t = 0, so that the 
series for b(x) begins with the term constant x 2”"-!”™” and the powers of 
successive terms in the series differ by m/({—4mr). Hence, if u,(a) begins 
with the term x and continues with the terms 2+”, a2+2?,,..., b(a) will 
contain the terms 22), 24+?, a224+2P__.. and so m $P/(N+1),7r = 2N/P 
with the restriction that r is an integer. It should be noticed here that the 
power series for i,(7) must be of the same kind as that for b(#), because 
of (24). Therefore the choice of NV and P automatically determines the 


form of the power series for i,(2). In particular, for the case when 


which corresponds to low-speed flow past a symmetrical blunt-nosed 
cylinder, m = 4, r= 1. Also, for u,(z) = B,—B,x, m = 3,r = 0. Both 
these examples are considered further in following sections. The case 
m 1, which is particularly simple mathematically but not very important 
in applications, is discussed in the Appendix. 

Before proceeding to the examples themselves, it is desirable to make 
some general remarks about the application of the foregoing method. 
There is no doubt that it is best suited to the construction of solutions, 
and their associated main streams and wall-temperature distributions, 
starting from chosen forms for i,(t) and b(t), i.e. from chosen values of the 
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x’s and §’s. If, on the other hand, the series 7,(x), b(a) are given, it is 
necessary to invert them in order to obtain the series i,(¢), b(t), whose 
coefficients are the «’s and §’s. As the examples will verify, this inversion 
is a drawback to the method from the point of view of the computer. 
Furthermore, in order to describe the distributions of velocity and 
temperature through the layer, it is necessary to use confluent hyper- 
geometric functions of the types ,/,(a, 4, 7), .A\(«, 3, 7), ete., and these do 
not appear to have been tabulated. If, however, only the skin friction and 
the heat transfer at the wall are required, the main computational task is 
the calculation of the k,,.’s (equation (30)). Now, k,» is easily found 
(equation (31)), and for uniform wall temperatures, «a, = 0 (7 > 0) so that 
no other k,, is needed. For this reason, only uniform wall temperatures 
are considered in the subsequent examples. Otherwise, to determine an 
adequate number of /;,..’s would involve the determination of, say, f, 
(n = 1, 2...., 5) and q, (7 ie ee 5) by integrating the differential equa- 
tions which they satisfy. Once these functions were determined, it would 
be relatively simple to calculate the /,,,’s and so extend the application 
of the method to non-uniform wall temperatures. Incidentally, this pro- 
cedure would also provide the velocity and temperature distribution 


functions. 


3. The solution for low-speed flow 

If M, ~ 0, but 1—7)/7, is not small, it is necessary to take account of 
compressibility, not because high speeds are involved but because of 
temperature variations. The boundary-layer equations (1) and (2) remain 
unaltered, but the term up, u,du,/dx and the dissipation term p(éu/éy)? 
disappear from (3) because they are O(.V?) compared with the surviving 
terms. The resulting equations are not so easy to solve as the better-known 
equations for incompressible flow (1, =~ 0, 1—1t)/i, = 0), because the 
momentum and energy equations are still linked through the density and 
the viscosity. However, at least for problems of forced convection in which 
the wall temperature is uniform, the low-speed compressible flow can be 
found, even without the restrictions that the Prandtl number is | and the 
viscosity is proportional to the temperature, for those main-stream distribu- 
tions for which the incompressible flow solutions have been obtained (7). 
The problem is the same and involves the solution of two ordinary differen- 
tial equations. The extra difficulty in the compressible case is that these 
differential equations are simultaneous, whereas they are separate in the 
incompressible case. 

In the present work it will be seen that 7 h, since u?/i O(.M}) is being 
neglected. Also, for the same reason, 7, = constant = i,in the mainstream. 
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The quantity 6 is therefore the square of the main-stream velocity divided 
by the constant 7,. 

The solution developed in section 2 applies to the case M, = 0 with one 
or two small modifications, which lead to considerable simplification as 
regards computation. The velocity is given by 


and the enthalpy by i i(1 q), 
where f, g are given by (27) and (25) respectively. The expressions (32), (33), 
and (34) for tr), x, and q, still apply, provided the factors involving 

l >) a 


n 0 


re removed in each case 


As an ex imple consider the problem in which 
U,(2) Bx, 
l—2,(x)/2, A (| Xo). 
) 


Then all the B’s vanish except £,, and m SN /(N+-1), so that 


: 38[2? w, p03 Bi(N+1)(ky Kg A)* |=, 
and rol ——,)|* = 2-431(N+1)'(k,—hy9 A)!, 
} 


( 


i‘. [ 2#3!(4)!]-"(u, p,)'i, A BY(N +.1)8(k, —hyg A)ietO-®, 
The heat transfer can also be expressed in terms of the Nusselt number 
defined by 


Nu = [ qo(€) dé /ug(ig—i,). 


With the Reynolds number defined by R U,(x)xpP,/ Me, 
Nuk 2134(4)!]-1(N +-1)-*(k, —ky_ A)!. 
In these expressions 
SN \, 2\ ,) || bem ( 1\,) 
b= | 1, 
; (ax 3} | 3) }/ \\3N+3 ; } 
((6N—2 | AE 1 ( 2\ ,) 
/ oy |e —-—}! {| —=]}!)], 
” ’ 5 5 3) |/ \\3N-+-3 3) J 


which follow from (31). 





The effect of wall temperature on skin friction and heat transfer is 
illustrated in Tables 1 and 2. The numerical values shown there are based 


on i,/¢,, = 288° K. (60° F.), and i9/c,, = (i) 273° K. (f.p.), (i) 373° K. (b.p.). 
When the wall temperature is at freezing-point A = 0-0521, and when it 
5092.25 


Cc 
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is at boiling-point A = —0-2951. When A = 0 the results apply to both & 
the incompressible case and the compressible case without heat transfer. al 
rn 01 
TAB eE | (for low-speed flow) 
al 
x \% 
Values of 7) {———,)*, where u, = Ba m 
Ms Ps Uy 
N 00904 | —0'0654| 0 3 | |) 2 | 4 large 
Aalsies sia MMi | : 
A=o0521 . : : 0° 109 | o180 | 0-312 0-468 | 0686 |: 112 | 2°165 | 1:073N* 
A=o. ‘ Oror | O175 | 0-312 | 0473 | 0°698 | 1°136 | 2-218 1-100N# re 
A = —o'2951. ‘ . 0045 | O145 | 07312 | O-502 | 0°764 | I°270 | 2°511 | -252Nt 
Exact values ford =o . | ° 0°157 | 0°332 | o512 | 0-757 | 1233 | 2-405 | 1'193Nt 
F¢ 
The second and fourth lines, referring to A = 0, are quoted from Light- 1- 
hill’s paper (4, p. 376). 
TABLE 2 (for low-speed flow) 
Values of Nu R-}, with u, = BaN 
N 00904 | 00654] o | 3} | 3 | I | 4 | large 
A=oogat . «| 0498 | 01578 | 0:663 | 0-708 | o-712 | 0-638 | 0-433 | oor F It 
1 
A=o. ‘ ‘ ‘ 0485 | 0572 0663 | 0-710 | 0-716 | 0°643 | 0°436 | Soren t 
4 0'2951. : ‘ 0°371 0538 0°663 | 0°725 | 0°738 | 0°667 | 0-455 | r-054N : 
Values for 4 = o deduced 
from exact values for | 
Prandtl number o77 0°493 | 0609 0°659 o0°671 | 0649 | 0°557 | 0° 366 0835N-4 i wh 
The values in the fourth line are deduced from the exact values for f 
Prandtl number = 0-7 given in Lighthill’s paper (4, p. 367). The deduction 
is made by using the rule that Nu R- varies as the cube root of the Prandtl 
number. The resulting figures in the fourth line should be close to the exact 
values for A = 0 and Prandtl number = 1; for instance, for N = 0 the 
exact value is 0-664, and for VN = 1 it is 0-570 (see Modern Developments in 
Fluid Dynamics, pp. 624, 632). 
The tables show that the non-dimensional coefficients for the skin friction 
; ‘ : = whe 
and the heat transfer are not affected by changes in the uniform wall ; 
In t 


temperature for flow past a flat plate (V = 0). This result is not peculiar 
to the present approximate method of solution, but is valid for the exact 
solution of the low-speed flow. In both tables the trends, for varying A, 
which appear for the non-zero values of V, are doubtless more reliable than | whic 
the absolute values. They show that when the wall is heated the skin friction If 
and heat transfer are increased for accelerated flow and decreased for re- | lengy 
tarded flow. In other words, the effect of pressure gradient, in respectively 


thinning or thickening the boundary layer in these two cases, is enhanced | 
tT N 


when the wallis heated. Now it has been realized for some time that pressure me 
j 1e 
' 
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gradient produces such thinning or thickening mainly by speeding up or 
slowing down the slowly moving fluid very near the wall (and so contracting 
or dilating the stream tubes), and this makes the results easily understand- 
able. For when the wall is heated the fluid near it is lighter and is therefore 
more readily accelerated or retarded by a pressure gradient. 
4, The boundary layer of a symmetrical blunt-nosed cylinder 

The main stream round a cylinder whose section is symmetrical with 
respect to the direction of the undisturbed stream is given by 

U(x) B, «+ B,2°+ B,2°+... 

For simplicity, we again assume uniform wall temperature, i.e. that 
1—i,(2)/h, A | Xo). 


Then, according to section 2, m 4 and r 1, so that 
t 
91 f s 
34(2%u, pt)? | €4(1+ 4B, 4+ 48, 8+ ...)* > 





(ky — hy A)B, + (kg—kegg A) Bo E48 + (kg ky9 A)Bs €4+-...]-* dé. 


It follows that B, 34 (Qty. pi.) (k, — kip A) By, 


 ipe B31, 
B,/ BF (1 +2 3s) /0 +13); 
where 7; A (kg —K a9 A)/(ky—ky) A), and 
‘ l hos , >. «nt, B,i,\? 
Q /R3 73 2) 10 | 3y = ee be a 
Y3/P4 (] r,)2| Yy— ie + (3 3") BS ( 3) | 
TIT Ett artnet 





where rs (kg— kao A)/(ky—ky9 A). 
In these expressions, 
Ly ' 2\1 4 2 ! = ' 
by = ep bu = Raf — ae 
(37 3) ( (3n a)! ( 5)! 


which are obtained from (31). 
If the non-dimensional coordinate 6 = x/a, where a is a representative 
length, is introduced, the shearing force on the cylinder is 


an ) 


ry = 2-434 (un, p,)'(ky—hy) A)? BY a0(1 +7, 62+-T, 04-+-...), (35) 


+ Note that the magnitude of the said pressure gradient is determined by the acceleration 
of the fluid in the main stream, whose density is not reduced. 
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here T, _ Ba, 1 Bi 
where 1 l+r, B, ltr, 4, 
and 
p — 12 Bsat [iri t+re(—33r,+ 15rj)] (Bsa)? 
© +r, B, (1+7,)?(1+ 7.) _ 





2 
‘. ry 


“603 296, 655.2 » (7571 35% 2 
- |i 1511 8/1 l'o( Sag 3714 va) (7c) 
— > ‘ 
(1+7,)°(1+7.) { ‘. 


Similarly, the local rate of heat transfer from the cylinder is 
Vo [223#(24)!] \(u, py)tt.(ky— hy A)? A Bi(] +Q, 6+ Q, 64+...), 


3 
Q» 2 3 » 2 2 
0) B. a* 7 “ B a~ 
where Q, u _ Aaa Ele Boas 


53 8 Pe | - (263 91, 635 2 2 
+ [oe “S0M + 37) tre — En +$ as Bi a* B,a* 





and 





(4+ 187, +4777) +7,(52— 10r, +12r?)] B2 a? Ba? 
i, 3B, 


(55 — Sr, — 437?) ] ( B? a? 
I T 1s) i, ): 


The case when the cylinder is circular is particularly interesting, because 





it has been studied for incompressible flow by various authors beginning 
with Hiemenz (see Modern Developments in Fluid Dynamics, pp. 148-50). 
In this case, a being the radius of the cylinder, 6 is the angle at the axis 
measured from the forward stagnation point, and Hiemenz’s calculations, 
based on the exact boundary-layer equations, give 
T) 1-2326(u, p,)? By aé > 
x ! 19.351 23% g21 (3.19 25™ 4 0: 55( 2") lo bon 
' B, | B, B, 

The present method in the incompressible case (B?a?/i, = 0, A = 0) gives, 
from (35), 


T, 1-1360(y, p,)* Bi ad . 


4 2) 
x 4-2-3 5 Bs a? 624 ly 19 Ps™ 4 0-5 i(=3 “) \ga4 
B, \" B, By) | 


to the same order of accuracy. Hence, although the present approximate 


method underestimates the skin friction for incompressible flow every- | 


where by about 8 per cent., it predicts very satisfactorily the shape of the 


function 7)(7), and in particular the point of separation (7, = 0). 
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With this in mind, it is tempting to examine the effect of wall temperature 
and Mach number on the position of separation. For this purpose it is 
assumed that the non-dimensional parameters of the main-stream velocity 
distribution, B,a*/B, and B,a‘/B,, remain fixed while B?a?/i, and A 
are varied. In the low-speed experimental results used by Hiemenz, 
a 4-87 em.. 

B, = 7°15! B, 0-04497, Bb, 0-0003300, 
so that B,a*/B, 0-1491(5), B,a* B, 0-02596. 
It is easy to derive, from Bernoulli’s equation, that 


bois (y—1)M? 


Bya*/a ~\ 2 , een e a/5\4 
(=) | Baa (5) 4 Bsa (5) (2 1 y—l We 
2 B, \2 B, \2 2 
La 


where VW is the Mach number of the flow at the shoulder (@ = 47), and i, is 
the enthalpy at the forward stagnation point (@ = 0). 





TABLE 3 


Value S o] 4 at the point of SE paration 


fe) ol os 
I*402 1°404 1°459 
1°429 1°43] I°Sit 
I*550 I*556 


From Table 3 it is clear that there would be no appreciable difference in 
the position of separation between the cases when the cylinder has tem- 
perature 0° C, (A 0-0521) and 100° C. (A 00-2951). 

It is interesting to examine these results in the light of the remarks at 
the end of section 3. From those one would expect that a retardea flow 


TABLE 4 


(te)4=—1— (oa =0 
4 By a (To) 4=0 
‘ 
° 38°5 
0198 38°4 
0° 390 3o°2 
0'566 37°6 
O75 26°5 
0°825 34°0 
0°860 31°5 
0°575 26°6 
0°570 I 3°D 
o'osi 55°O 
0850 100 (separation for A = —1) 
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would be not negligibly influenced by wall heating and would separate 
earlier (owing to its general reduction of skin friction), and this supposition 
will be confirmed in section 5. In the present problem the effect is only 
small because the flow over the first 90 per cent. of the unseparated boundary 
layer is accelerated. Hence, as Table 4 shows, the skin friction has time to 
reach a considerably larger value at the pressure minimum than for the 
unheated cylinder; and though thickening subsequently takes place more 
rapidly than for the unheated cylinder it is insufficient to cause much 
change in the point of separation. 


5. The retarded main stream uw, = B,—B,x 

The boundary layer for an incompressible flow with the main stream 
given by u,(x) = B,—B,2 has been studied by Howarth (5), who found 
that there would be separation at X = B,2/B, = 0-120. In order to 
obtain this result, accurate to 1 per cent., Howarth found it necessary to 
employ nine terms in the series for 7(x). 

In the compressible case, with 


u, = B,—B,x = By+f,t-*+f,t-F+..., 


1—i,/hy = Ag +A,x+Agx?+... = apa, tt+a, #4 
“ 
ax = 34(2*u, ptt)! | €4(1+48,+48, +48, é8+...)!x 





x [hig Bo+ (ky — Fx %0)B1 $4 + {(h2— heap %)Bo— hay % By fE* +...) dE. 


The first term in the expansion in series on the right-hand side is Ztt, where 


i= sige pst8)*(1+3B)*(ky Bo) 





Bi (—4)1)[2* 4 (0)04(0)f8,(0)}4 85)" 
With the substitution z = B, Ett/B, it follows that X may be expressed 
as a power series in z, starting with the term z. Therefore, let 
u, = B(1—X) = B,(1—z+a,2?+a,23+...), (36) 
where the a’s are to be determined. The expression for may be rewritten 
x= EK fat > raz") (14 ¥ 3,2") : (37) 
0 
where 
n—1 
7 :-™ ( B, )" ; ky, $,,- 2 bne 06 rPa- “(3 B, ie (38) 
——$——___— — \t 
sia +By B, E ‘ ko Bo B, E 
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Equation (37) as well as (36) gives x as a power series in z, and by comparing 


the two series it follows that 





Li an pitt 10% 4 
a2 o(C,+D, Ss 3(C,+C;, D,+ Dy) 
a, 1(0,+C, D,+C, D,+ Ds) 
as 3(C,+C3 D,+-C, D,+-C, Ds+- D,) 
Ag —$(C; | C, D, T C3 D, +C, Ds+ C; Dd, T Ds) 
where 
; } 2 ’ bf re 
C, fy; C, fo Oy7, Cs; 4y5+ l2y, y2+ 41 
7” ),,2 | 
Cy = 444+ 12y, yg + 6734+ lyf vet yi 
— ‘ Dny2 te. 
CU; fy 5+ 12y, ya l2ye73+ l2yi ys 271 Y2t+ 41 ¥2 
and 
" 1 392 15 35 5 
D; 304 D, Og 31 D; 203 + 49) 02— 16 
3 392 5925 35 94 
D, 554+ $5, 53+ 353— 457 5, + Hf 
l 3 5S2S8 BS $2 3583 $3 
D; 995 0,0 £9293 1857 53— 485, 45 320] O 32504 
Moreover. 
1. 
] Lh r ~ > > 
i.—tut =. 1, —$. BA(1—z+a,2?+...)? 
With z 0, this becomes 
1. 
1+4B, s _, 
2 1 R2 
8s SBS 
2 = 
- . _ a ree Cfo © 2 
o that 1+46 l—- z (l1—z+ a, 2*-+ 52 | ; 
="T Fe 
@ 
A 1 ¥ alt 
But 1+36 = (1+ bBo)(1-+ 2 Yn?" ). 
n 


Comparison of (40) and (41) gives 
Y1 Bo, Y2 12 Bo B5: Y3 Es Bo - 2F, B5- 
ya = E,By—(— E}+ 2K3)B5+ 3E, Bo+ Bo 
E; By —(—2E, E+ 2E,)63+ (3£,—3E2)p3—4E, pé 


where : ' 
E, = a,+4, E, = ag—a, 
7 | 192 } \ 
E, = a,—a,+ 4a, E, = d;—A4+_43 


>. (39) 
1 
(40) 
(41) 
] 


—Bor. (42) 


| 


J 


Similarly, by substituting for x in terms of z from (37) into the identity 


v { yn >* By i myn 
n=0 \B,E) ©’ 
n= 
the values of «,(B,/B, £)” (n = 1, 2, 2 are found in terms of the y’s, 


5’s, and A’s. Using this last procedure together with the relations (39) and 
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(42) we can determine a,,, y,,. 5,,, «,(B,/B, 2)" for all n when the values of 
8B, = B, and the A’s are known. 


Next, if 1) = FX-#(1+ > 7X"), 
r=] 


the substitution X = z—a,z?—a,z°—... provides a power series for 7, in 
terms of z which must be the same as the expansion 


ro = 213 -Ap, pyigt (144 S By) * Y (ey Bu— >, ney) t 
2-§3-414(0)p,(O)ky Big) - ‘(14 2 yn2") *(14 > 5,2"). 

Hence, F = 2-*34[(—4)!]-4[n,(0)p,(0) |! BE By. 
T, 
T; 
T; 


BS SS 
SY 
— 
~ 


Mb +34Fe—hG.—G 


9. 


= 
SS 
oO 


F, 6,+-63+ 1, a3+-2T, a,—(T,—T, a,)G@,—T, G,—G, 

T, = F,+ F36,+F,6,+-F, 6,+8,4+ 7, a,+ 7, (2a,—a3)+ 37, a, 
(73;—2T, a,—T, a3)@, —(T,—T, a,)G@,—T, G,— Gy, 

T;, = F,+-F,5,+F5,.+/8,4+F 5,+-5,+7, a,+27,(a,—a,a3)+ 
3T;(a,—a3)+4T, a,- 





gal 97 Me 2 Y 
T",—3T; a.—T,(2a,—a3)—T, a,|G,— 
(T;,—2T,a,—T, a,)G,—(T,—T, a,)G,—T, G,—G;, 


where 


F, BY 1: F, 37274 evi. F; 33+ m1 Y2 wee 

F, Vat Py Yat Sve evi Y2— ee YL 

Fy = Syst Py vat Pye Va evi 3 — 6-1 Yat Sei Yo Beek 
and 

Gy $A, Gy 03+ 3p, Gs 44+ $aya3+ fs 3 

Gy, 305+ $a 4+ 05+ 4.05 as+ io05 

Gs; 3g + $y 45+ $a, 04+ Baza, Ba, a§+ 3303 ag+ $408. 


The computation of results is simplest when the wall temperature is 
uniform, in which case the a’s are all zero except a) (= A). Then 


5, — hected B Ba)" _ Fahad 24, 
,; ¥ ko Bo B, K ko Bo ’ 
where 
(4n)! (—2)! (4n)!(—#)! (4n)! 
k —tme ann Od &. s a A EO sil SA : 
~ (ga—gyl(—4)! m (4n—8)!(—4)!_ (4n—})! 


This is the case for which numerical results have been obtained. 
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3efore proceeding to these results it is appropriate to mention the choice 
of the constant A. The stagnation enthalpy of the main stream is 
h, v,(1+ $b) 1,(0)(14 $Bo). 
where 7,(0) is the enthalpy of the main stream at the beginning of the 
boundary layer, and A | —i,/h,, where i, is the enthalpy at the boundary. 
. ‘ 0 1 0 . . 

Hence A l a , 

; i 

(,(0)(1 + 2Bo) 
Suppose now that the boundary wall is initially at rest at the temperature 
i,(0)/c,, of the surrounding air. If the wall is accelerated so as to attain a 
steady speed B, relative to the air at the station x = 0, there will be a 
period of unsteady flow, longer than the period of acceleration, during 
which the temperature of the wall increases from 7,(0)/c,, to 7,(0)(1+ 3B5)/c, 


The latter is the final temperature of the wall, at which heat transfer between 


,° 


the wall and the air ceases. If the wall is the surface of a body with a large 
heat capacity, the unsteady régime may last much longer than the period 
of acceleration. In fact, as an extreme example of this kind of behaviour, 
one may envisage a quasi-steady state in which the wall temperature is 
still at its initial value 7,(0)/c,, when the final speed has been reached. 
Under these extreme conditions, 


A _ 
1+ Bo 2+ Bo 


and particular values of A are: 


° I VIO 10 

Q lo o-4 , 40 
} 1 2 20 

° 6 3 21 





These are examples of a cold wall. As an example of a hot wall, one may 
0) with 7,(0)/c, 
tained at the temperature of steam. Then A 


take low-speed flow (/,(0) 288° K., and the wall main 

0-2951, as in section 3. 
The case A = 0, with no heat transfer, will be used to introduce the 
numerical results, because exact results are available in this case for pur- 
poses of comparison. The first five 7).’s have been calculated, and they are 
given in the following table 


TABLE 5 








Coefficients in the series for t,, when A = 0 
1 C I NIO 10 
; 183 | 020 1550 10°152 
7 | 10°762 74°23 1918 
7 71 445°5 5010 x 108 
108°C 5389 3*106 x 10® 
1851 II 7°539 x 104 2'021 x 108 
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The results given in the first column of Table 5 may be expressed by means 
of the series 
e 315 8 / >] 5 
To = 0°312[ 4, (0)p,(0) Box]! x 
x | 1—0-7728Y —0-0792Y?—0-0727 Y3— 0-0612Y4— 0-0565Y°—...], (43) 
where Y = 8X, whereas Howarth (5) gives 
ae (226 3/y]3 
Ty) = 0°332[4,(0)p,(0) B3/x]* 
< [1 —0-7683Y —0-0521Y?—0-0422Y3—0-280Y4—0-205Y>—...] (44) 
followed by further terms, in the same low-speed case. With either series, 
the position of separation, viz. the smallest positive root of 7)(Y) = 0, may 
be circumscribed, at the stage where n terms have been retained in the 
series, by finding the roots of tr, = 0 first under the assumption that 


1=0 (>) (45) 
and then under the assumption that 
T=—T* (¢ >), where TF = 8-T.. (46) 


With n = 3. 6 the roots are 








Assumption (45) Assumption (46) 





| 
| 
‘ . . ! 
Equation (43) | X = 0°145, 0-122 | O°105, O°109 
Equation (44) O'150, 07134 | oI Io, o-116 





respectively. Now, Howarth (5) finds X = 0-120 for the exact root of (44), 
which, it will be noticed, is about one-quarter of the way from the lower 
bound to the upper bound of the root whether m = 3 or 6. On the rather 
crude assumption that there is similar behaviour in connexion with (43), the 
present method estimates that separation will occur at X = 0-112 for low- 
speed flow with no heat transfer (incompressible case). This gives an error 
of —7-5 per cent., which may be compared with —15 per cent. for the 
Karman—Millikan method and + 30 per cent. for Pohlhausen’s method (see 
Modern Developments in Fluid Dynamics, p. 171). 

In turning to high-speed flow, the procedure just outlined has to be 
modified slightly. The merit of the figure 8 occurring in the definition 
Y = 8X, is that it leads to values for the 7¥*’s which slowly decrease in 
magnitude beyond a certain stage (r = 4 in Howarth’s series), and it thereby 
justifies assumption (46) as a means of setting a lower limit to the position 
of separation. The figure 8 is not, however, appropriate for all the cases 
considered here. Therefore the following has been adopted as a standard 
procedure: 


(1) calculate X,, the lowest positive root of 7)(X) = 0, under the assump- 
tion that T, = 0,r > 5; 
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(2) with Y = X/X,, 7* = X$7,, calculate the lowest positive root Y, of 
i+ > Tey’ = 0, where T? = TF ir > 5); 
r=1 
, Xo - : , ; , 
(3) let X,., r (1+-3Y,) be the estimated point of separation. 


This leads to the results given in Table 6. 


TABLE 6 (A 0) 




















M4(0) ° I | VIO 10 
. si - | 
Xo 122 orlI2 | 0°0675 00162 
[? 0°7543 06742 0°3071 +o°1644 
TS 0°07542 0°1350 | 0°3382 | —0*5034 
Te 0°06759 0°07546 0°1370 | —O'2130 
4 0°05558 006429 | —ollIg§ | —o'2139 
[* 0°05003 -0°05835 o'1050 =| —0°2255 
mA 0°107 0°097 00578 | 00135 
iti, ssctiaiaeaieaaal <n, CNS 
K sep oli oIol | ©0602 | O°0142 
X (Stewartson) 0°120 o’11o | 0°07 0°02 
SE p ‘ 





The last line in Table 6 is quoted from Stewartson’s paper (6). He ob- 
tained these results by correlating the compressible flow u, = B,—B,2x 
with an incompressible flow, and solving the boundary layer of the latter 
flow by Howarth’s approximate method (5). The only approximation 
involved is that of Howarth’s method, and, as Stewartson says (6, p. 98), 
the true values of X,,,, should be lower than those quoted, except at 
M,(0) = 0. The case (0) = 0 suggests that the values in the penultimate 
line of Table 6, obtained by the present approximate method, are too low, 
so perhaps the true values of X,,,, are between the values given in the last 
two lines of Table 6, except for (0) = 0. Although the proportional 
deviation of X,,,, from X,,,, (Stewartson) increases fairly consistently with 
increasing Mach number, the reliability of the values of X,,,, decreases at 
the same time. For a glance at the 7*’s in the last column of Table 6 shows 
that more terms would have to be calculated in order to achieve the required 
decrease of 7* with increase in r. 

In this case of no heat transfer, the distribution of velocity and enthalpy 
through the boundary layer can be found when the functions f,(7) are 
known. For both the linearly decreasing main stream and certain other 
main streams, the required functions are 


(4n)! ,(/4n-+1 1 
n) ae 6 nU ——-.,-, 
Jn\ (—})! ( : = 7 


3 


and values of f,, fo, fs, f,, and f, are given in Table 7. 
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Two examples of flow at low speed past a hot wall are given in Table 8. 


TABLE 8 














M,(0) 0 
0°295 1 I 0°2951 
7175 I? 0°707 0°7355 
ie O° 10014 008571 
° T¥ | —0:07636 0°07129 
ie 0°06308 0°05875 
42.4 T? 05211 004804 
071025 | Xo Yo 0°0649 0°0900 

ser 0°093 
Table 8 shows that separation is estimated to occur at X 0-067 for a 


‘ 


wall maintained at the uniform temperature of 303° C., whereas the same 
method of estimation gives X 0-111 for the wall at 15° C., the tempera- 
ture of the main stream. The table also shows that raising the temperature 
of the wall from room temperature to boiling-point (A 0:2951) would 


cause the separation distance to decrease by about 16 per cent. Table 9 


TABLE 9 








0-2 2 
A 2 A= 
I | NIO | VIO 
5°233 3 2°174 
o°197 23°40 
59 65-2 
4 210°d 128°2 
_ _— 4 eins 4 
1393 4006 x 10 1°167 X 10 
0130 0°0827 0°140 
06787 0*2094 +-0°3037 
? 0°13788 0° 3815 0°4578 
: o’o7III 0°1444 0°1778 
07061 36 o110g 0°04888 
: O’OSIII | 01549 | 06211 
} 
oO'll4 00090 obs @ fe) 
o7118 0°0729 O'1I7 
o'lol 0°0602 00602 





shows an effect of similar degree concerning the greater separation distance 
for a wall colder than the main stream. The effect of wall temperature on 
separation is much more marked in this case of a linearly decreasing main- 
stream velocity than in the previous case of flow round a cylinder, and it 
is perhaps sufficiently strong to be observed in experiments. 

As regards relatively cold boundaries, attention has been restricted to 
one or two values of A, and different Mach numbers. Table 9 refers mainly 
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to A = 0-2, chosen as representative of the flow at M,(0) = 1 when the wall 
is kept at the temperature of the incident stream. 

The results for M,(0) = 10 have been omitted because it is clear that 
already at M,(0) = v10 the coefficients 7¥ are no longer decreasing in 
magnitude as r increases, and therefore that more terms are required in the 
series for 7, in order to obtain a reliable estimate of X,,,. The last column 


0-13 


* M, (o)=vi0} 


Ringed points are 


Cewartsons 
06. 56-2 
08 05 0-7 09 11 13 
t, (0) _ \nitial stream temperature 
to Wall temperature 
Fic. 1 
illustrates the trend of the results when A is still larger. The value A = 


is appropriate for a flow at 1,(0) = v10 with the wall at the temperature 
of the incident stream. Again it is obvious from the values of 7'* that more 


terms would be required in reaching a reliable answer for X Even so, 


sep’ 
there does not seem to be any noticeable inconsistency in the various values 
of } 

Acceptance of the figure 0-117 for X,,,, in the last column leads to the 
following result. If the temperature of the wall is always kept at that of 
the incident stream, the position of separation will not change very much 
as the Mach number increases, for at M,(0) = 0, 1, v10 the corresponding 
O-111, 0-115, 0-117. 


This, if accepted, strongly confirms the theory of section 3 that the degree 


values of A are then 0,'3, 3, and X 


sep 


of influence of a given pressure gradient on a boundary layer depends con- 


siderably upon the wall temperature, because this fixes the density of 


the fluid very near the wall, and through it the effectiveness of pressure 
gradient in slowing down or speeding up the fluid, and hence the thickening 
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or thinning of the boundary layer. It is now seen that this conclusion is 
almost independent of the Mach number. 

Still greater confirmation is obtained by plotting all the values of X,,,, in 
Tables 6, 8, and 9 (excluding cases with M,(0) 10) against the ratio 


2,(0) - ] 


' lo - (1- A)(] T $B) 


of the initial stream temperature to the wall temperature. 





It is seen in Fig. 1 that the relationship may be represented by a curve 
nearly independent of Mach number (the broken line), especially since the 
point X,.., = 0-0729 is ess trustworthy than neighbouring points (see 
Table 9). The figure also shows Stewartson’s points. By means of these, 
the relationship between X,,,, and 7,(0)/ig which has just been obtained 
can be corrected for the tendency of the present method to underestimate 
X..p- Taking into account that Stewartson’s lower values are probably 
slight overestimates we can suppose that the full line (which has been 
drawn roughly parallel to the broken line) gives a reasonably good descrip- 
tion of the effect of wall temperature on separation in uniformly retarded 


flow. This full line has *ae equation 


xX 0-120[7,(0)/t9 |}. 


sep 


or, in other words, 


relative fall in velocity before separation (per cent.) 





2( mercer stream temperature) 4 

“\ wall temperature ) , 
One may reasonably expect this rule to be equally good for Prandtl numbers 
other than 1.7 

At the same time, Table 9 shows that this result has not been verified for 
large Mach number. To test it at high speeds would require the introduction 
of further terms in the series. 

The general conclusion that the separation distance increases as the 

temperature of the wall is lowered has previously been obtained by 
Dr. N. Hoskin at Manchester University in unpublished work. Of the 
present quantitative information, perhaps the result to be stressed, as being 
the most trustworthy, is the result in low-speed flow that the separation 

+t Brown and Donoughe (7) find, when calculating flows with large temperature differ- 
ence, small Mach number, Prandtl number 0-7, and u, 0 2~¥, that the separation values 
of N are 0-0600, 0-0904, 0-1178, 0-1351 for 7,/i) = 4, 1, 2, 4 respectively. Now 

N (a/u,)(du,/da) 

bears a resemblance to Xsep, and so it is interesting that it increases with 7,/i) in a way 
similar to X¢ep; also its slower rate of increase for values of 7,/i, greater than those in Fig. 1 

} may be significant. 
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distance should decrease by about 20 per cent. as the temperature of the 


‘ 


wall is changed from 0° C. to 100° C, 
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APPENDIX 
The solution when m l 
When m 1, we have N r/(8—r) and P 2/(§—1r), where r is an integer, so 
that the solution applies to main streams such as 
U(x) B,- B, x B,xi-+..., 
u(x) B,xvi+ Bax? -+ B,a+..., 
U;(x) B, x + B,x*+- B,x®+..., 
with wall-temperature distributions expressed in the forms of series of integral powers 
of xt, x, x3, respectively. 
The reason for singling out this case for special mention is that the solution for h, is, 
at least formally, simpler than the more general solution given by (29), for (28) is 


satisfied by h, , Crd tnt » Inr Gi provided that 
r 0 
n n A 
3 d,-l 4g; rv (n— r)917] 7 Zz (n T)otp B n_ rQJi.r> (47) 
rT 0 r 0 
where the extra suffix 1 has been attached to f, g, h to emphasize that m 1. In 


obtaining (47), it is necessary to use the equation 79/,,+(n+3)9i,—T917 = 0. Now, 
starting from the formula 


Jir > e*(E n)'€3- dé, 


it is easy to verify the recurrence relation 


(3r—1)9, -= 3r(g).,- 


—) 
% 
| 
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from which it follows that 
n!' ’ ? 7 1)! n! ( 1)! , 
iP ' —— Jin — Gio (n> 1) 
‘ (n ty! oy (nM )! , (n 1)! 
Hence (47) is satisfied if, for 7 l, 
l ( 1)! . gs! : 
l 3 7 -— dy (n—r)a,B, (7 0, 1,...,2—1), (48) 
< 7 ,-_ ) 
- y! 
und 1 1)! 7 = ap 0, (49) 
hang (7 3)! 
r=] 
which suffice to determine d,,.. (7 _ Serer n). The boundary condition f = Oaty%s = 
is satisfied because f, ,,(00) g 0) 0 so that hy ,(0) 0 for alln. The boundary 
condition > B, t” at ub 0 is satisfied if 
B, > Any, (50) 
¥ r=0 
for all n. 
The first few values of the c’s and d’s obtained from (47)-(49) are 
- O ( s 
dio Og 1s di, Xo Das Cy (1 -4.0%9) P13 
199 SX Po, yy $09B2+3018,, doe Xo B2— FB, Ce (1—$a)B2— to Bi; 
dso 298s, ds) F5% P3+ $a Bo, 3p 3% Bs T $a, By T 20 Bi, 
dus v9 B3 — 40, Bo—402P}, Cz (1—43a9)B,—440, B, 40 B. 


It remains to find x in terms of t, a process which involves determining 7) and lim 7th; , 








n>0 
as before. pince 
n . n n—1 
hi Cnfint D> 4nrgir and lim > dygi, = 3 ¥ (n—r)(apBn_-—Anr) 
r=0 7-0 r=0 r=0 
from (47), it follows that 
i P ‘ n!(—#)! 
lim 7h; c, lim nf fj,, ——— ) :™ 
7-0 n—0 (n 5)! ( 4)! 
Therefore, the equation corresponding to (33) is 
x 3¢[(— 4)!]2[22{(— #)}}4u, p, : 
: * $y—1)/2(y—1) / 2, a = 
t 
é} Q ¢n : - 
is (1 3 S B¢") (> maaicns d§. (51) 
\ pa nee H 
0 ape ’ n=0 . 
For an example, consider the case where the data are 
U,(x) Bzi, 
l ( v)/hy A ( OX)» 
for which f, 0,8, #0. If x Btzige" 3 


; isa non-dimensional measure of position 
along the boundary layer 


u u Fa 
b . aes oe X+44K*4-2¥?+.... 
BT 3 pu 1—3X ~ - 
D 


5092.25 
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Under the assumption that y = 1-4, 


t 
2 = 3i(— 9) Hy py) [ +28, E+ (28, +3 BNE 
0 
x (3c, E+ $c, €?+- Ste, £3 +-...]-# dgé 
2838 ( 4)!)3(5p5 py i} ct) 143 x 
; ; 5 aa 3 Cy _45 98.1 3—2_ 3 mg 273 27 Cs i 
[1 | 27 4 =) 17282 3 Bi— 2B, — C 392 ce 28 Cc IE ——_ 


Comparison of the x- and t-series for 6 now shows that 


B, = C = 2838[(— 4)!)8(5p, p,)- 21 BU(L- gir 


B, S, C?, where S, = 35/(22+9R,), R, (1—#2A)/(1—}A), and B, = S, C, where 
S, (24,3435 | 1688, 7 29109 PR S, | 71820 R? S?) (238 SiR. 2)» 


R, (1—38A)/(1—FA). 
The shearing stress at the wall is 


T) = 3[2%54{(— 4)3}2]-(n, p,)# BH(1— 4A) iat > 


18 115)y 243p 9 3321 p2 C2 116487p ¢ 

x [1+( rR, S, sa )X + (779 Ry S, piiahy Sj 16456 RS; 
169 ¢ 448795),\y2 
68 Sy "6 5824 )X 


and the heat transfer from the wall is given by 
qo 24[54(4)! {(— 4)!}4] (a, p,)47, A(1— A)? Bhat > 
< [1+ (ER, S\—38:)X + (498 R, S)— gts RZ S3— 3429, S,— 18S, + 48292)X24 
The corresponding expressions for incompressible flow are obtained by putting X 
and deleting the terms in 1— 4A. 
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CHARACTERISTIC SURFACES IN IDEAL PLASTICITY 
IN THREE DIMENSIONS 
By J. W. CRAGGS (University College, Dundee)t 
Received 29 May 1952] 


SUMMARY 


Conditions are found for the existence of characteristic surfaces for the equations 
leal plasticity in three dimensions. It is observed that the conditions are every- 
satisfied in the problem of plane strain, so that there is then a complete set 
haracteristics. In more general problems, including those with axial symmetry, 
there are only isolated characteristic surfaces. An example is given of an axially 
symmetric solution with an isolated characteristic surface. 


Introduction 

THE theory of characteristics as applied to the equations of ideal plasticity 
with plane strain has been extensively studied and is now a standard 
method (1 Hill (2) and Symonds (3) have discussed the equations in the 
case of axial symmetry, and shown, by different methods, that these are 
not hyperbolic, and so have no characteristics at a general point. The 
purpose of this paper is to show that characteristic surfaces may neverthe- 
less exist, not only in the case of axial symmetry, but also in general three- 
dimensional problems. Such characteristic surfaces are in general isolated 
and no general methods for the solution of special problems can be obtained 
by their use, but their existence does imply that of plastic rigid solutions 
with discontinuities of stress and velocity components, and such solutions 
ire more easily obtained by relaxation or other methods than the corre- 


sponding solutions for plastic-elastic problems. 


1. The conditions for the existence of a characteristic surface 
Write o,, for the stress tensor and wu; for the velocity vector (7, 7 1, 2,3). 


Then the equations of ideal plasticity are the Mises—Hencky yield condition 


“4 $ »)].2 
07; Fi; 2k 9 (1) 
where o;, represents the reduced stress tensor formed by subtracting the 
mean hydrostatic stress from each of the diagonal elements of o;;, the 
Lévy—Mises flow equation 
l jou cou ’ 
| - Z Ac; 9 (2) 
zZ\o2 Ox; 
and the equilibrium conditions, 
| 
CO; 
i _ (3) 
Cx; 
4 
Now at King’s College, Newcastle on Tyne. 
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A characteristic surface S may be defined as a surface such that, when 
the nine variables o;;, uv; are given at every point of it, the equations (1), (2), 
and (3) do not suffice to determine the normal derivatives of the variables 
across the surface. 

Choose axes at a point O of S with Oz normal to the surface and Ox, Oy, 
Oz forming a rectangular Cartesian set. Then, at O, it is required to find 
the first derivatives with respect to z, assuming the first derivatives with 
respect to 2 and y to be known. Let wu, v, w now represent the velocity 
components along the axes, then 

éw/ez is determined by the compressibility condition ; 
60,,/62, €o,,/€z, and éc,,/éz follow from the equilibrium conditions. 
The flow equations (2) reduce to 


ou ov ; Ow ou, ov 
= Adee, = Ao, -— Ao;,,, and —-4- : 2Ao,,, (4) 
OX cy ; C2 Cy OX i 
cu, Cw ov ow Gc 
and  L == ZAg.., —-+.—— = 2)d,,. (5) 
Oz Cx CZ oy a 


Suppose first that at least one of the stress components on the right-hand 
sides of equations (4) is non-zero at O. Then the corresponding equation 
determines A, and éu/éz, év/éz follow from (5). The incompressibility con- 
dition may therefore be differentiated with respect to z to give 6w/éz?, and 
éA/éz follows from the third of (4). The remaining equations of (4) may 
then be differentiated with respect to z, determining the stress derivatives 


0o.,,./z, Gc,,/0z, and éa 


LI ry 


yy/ez 8o that all the stress and velocity derivatives 
are determinate, and S is not a characteristic surface through O. 
Next suppose that, at O, 


, ’ 


’ > 
ze — Gyy Tz, = Ozy = V. (6) 


C 
Then (4) do not determine A, and (5) reduce to a single relation between 
éu/éz and év/éz. 

The yield condition, differentiated with respect to z, gives a single 
relation between the three stress derivatives @c,,,,/6z, éc,,,/0z, and @o,,,/é2, 
but one velocity component and two stress components have indeterminate 
derivatives with respect to z over S at O. The element of S at O is therefore 
an element of a characteristic surface. 

Now (6) will be satisfied for some elementary surface through O whenever 
the determinant of the reduced stress tensor vanishes at O. Any such point 
may be called a characteristic point. 

Finally, if the characteristic points lie on a finite or infinite surface such 
that the maximum shear stress over the surface is k, the surface is a charac- 


teristic of the equations. 
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2. Special cases 

In problems of plane strain, the determinant of the reduced stress tensor 
vanishes at every point. Every point is therefore a characteristic point, and 
two characteristic surfaces pass through every point of the plane. 

Two other special types of problem have complete sets of characteristic 
suriaces : 

(i) problems in spherical polar coordinates with zero radial velocity, and 
the other two velocity components proportional to the radius, and 

ii) problems in cylindrical polars with the radial velocity zero and the 


circumferential velocity proportional to the radius. 


3. Application to a plastic-rigid problem 
The importance of the existence of characteristic surfaces lies in the fact 
that such a surface may divide a rigid region from a plastic one, when the 
plastic-rigid approximation is made. It follows that the full plastic-elastic 
problem may be simplified, for large plastic flow, by neglect of elastic strains. 
\s an example, consider the spherically symmetric form of the equations 
of motion, with solution independent of the azimuth angle ¢ and the polar 


radius 7. The equilibrium equations are then 

(opp o44,)cot 6 30,9 i) (7) 
und 20, T99— hd) + Fpp COLO 0. (8) 
The flow equations, with w, v for radial and transverse velocity components, 


and primes for derivatives with respect to the polar angle 6, are 


0 yo tu u—+vecoté u’—v 
Be The A tees (9) 
20.—O99—O¢ 2099—F$4—Frr 2F44—Frr—F99  GOzp 
Two possibilities arise from (9). 
i) Take the numera?/ors all zero, then 
W cos 6, v W sin 6, (10) 
and the material flows as a rigid body with axial velocity W. 
ii) Take 2o Ona— (0. 
Substitution in (8 then leads to 
a9 Cot H ). 
7 ’ 
whence ki sin x cosec 6. (11) 
where the constant of integration has been so chosen that @ x is a 


char icteristic surtace 
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The yield condition now is 


$4—%99 = +2k(1—sin2a cosec”@)!, (12) 

and (9) reduce to 
vu utveotd@ u’—v 13 
7056 4s 10,9 ie 


Assume that the radial velocity component is greater than zero, then the 
rate of plastic work will be positive only if the upper sign is taken in (12). 

Substitution of (11) and (12) in (13) leads to a pair of equations for wu, v, 
the solution of which is easily seen to be 


v = {A+ Bf(6)}sin@ (14) 
and u (v’ +» cot A), (15) 


where A, B are constants and 


“ 
" fi 2 sin « 
log f (@) i 13 cot Ga ——— — | dé. (16) 
| sin 6,/(1—sin?« cosec?@) | 


. 


x 


(The integral can be expressed in terms of Jacobian elliptic functions, 
but it was found easier to evaluate it numerically for the particular case 
required. ) 

Now choose A, B such that for @ = 6, > 47, v = 0. Then the equations 
correspond to a material flowing past the outside of a cone of semi-angle 
a—6,. There is a characteristic surface at 6 = a, and for 6 x the material 
moves as a rigid body, with axial velocity W given by the condition that 
the normal component of velocity across the characteristic is continuous. 

The solution has been computed for the particular case « = }7, 6, har. 
This corresponds to the flow pattern when an infinite block of material is 
pressed against a rigid wall by a suitably chosen pressure distribution at 
infinity. The radial and transverse components of velocity are tabulated 
below, and streamlines ABC, A’B’C’ are shown in the diagram which 
illustrates the flow on one side of the axis OS of symmetry. Notice that, 
just as in plastic-rigid solutions with plane strain, there is a discontinuity 


in the tangential component of velocity across the characteristic surface 


which forms the rigid-plastic boundary. 


TABLE 





6 | 30 30 40 50 60 | 7o 80 )0 





u | —o-944. cos @ | 17183 | 0239 | o-r4r | O-105 | O-081 | 0°066 | 0-058 
| 
0°944 sin @ | 0-472 | 0:230 | 0-142 | 0090 | 0°051 0°023 | 0000 
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Fl 














“a 0 


The stress distribution is given by (11) and (12), ogg being arbitrary to 
the extent of an additive constant. This may be chosen such that the 
relation between the shear stress, }/, and the normal pressure on the plane 


boundary fits any desired coefficient of friction. 
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THE VARIATION OF COMPRESSIBLE FLOWS 
By A. R. MANWELL (University College, Swansea) 
| Received 18 October 1951; revised 5 June 1952] 


SUMMARY 

The theory of the perturbation of those compressible flows in which there are 
only two independent variables is developed, with special reference to steady plane 
transonic flow. Expressions are found for variations of both the physical flow and 
of its hodograph and relations obtained between the two sorts of variation. It is then 
proved that even if the Jacobian of the hodsograph of a flow should vanish on the 
boundary the flow need not be critical. Using the expressions for variations in the 
physical plane of flow, boundary-value problems for perturbations can be discussed 
directly in the physical plane. The methods are applied to the case of circulatory 
flow outside circular or nearly circular cylinders. It is proved that instability can 
arise by a resonance process near certain critical speeds, measured on the boundary 
of the flow, and that at these speeds Dirichlet’s problem for the perturbation stream 
function has no unique solution. 


1. Introduction 


It has been conjectured that transonic flows may be unstabie even on the 
perfect fluid theory of compressible flow. In the attempt to demonstrate 
mathematically the existence of a singularity in the hodograph equations 
of transonic flow the behaviour of ‘limit lines’ has been studied in detail, 
cf., for example, refs. (1)—(5), but the ‘limit line’ singularity cannot appear 
in flows of the type usually investigated, as has been shown in (4), (5). 
On the other hand, Frankl (6) has endeavoured to prove instability by a 
general argument concerning the non-existence of solutions for contours 
adjacent to a given flow. Again, Guderley (7) has put forward arguments 
to show that singularities of the flow solutions arising at the end of a local 
supersonic region, due to initially small disturbances, will be accompanied 
by instability of the actual physical flow. Neither of these papers, however, 
gives any precise example of the breakdown of potential flow solutions. 
In the present investigation, in which variational methods are used 
extensively, the basic step is to show that the perturbations of potential or 
stream function satisfy some particularly simple equations in the charac- 
teristic plane, or better, if this is available, the hodograph plane. The 
functions which arise in the perturbation theory are distinct from the 
ordinary hodograph solutions but can be expressed by means of them in 
an elementary way. The equation for perturbation of both physical and 
hodograph flow is derived and it is proved that even if the Jacobian of the 
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transformation from the physical flow plane to hodograph space should 
vanish on a (special) boundary, such a flow need not be critical. The 
perturbation theory is then applied to the case of a transonic vortex, i.e. 
a circulatory flow outside a circular cylinder. It is found that for most 
changes of the boundary there exist corresponding small changes of the 
flow. However, near many speeds of the stream at the boundary there are 
critical perturbations which give proportionately large changes of the flow. 
The situation is closely analogous to resonance in linear vibration theory. 
The critical speeds themselves are characterized by the non-uniqueness of 
Dirichlet’s problem for the perturbation stream function and so for the 
stream function of the compressible flow in the physical plane. The solution 
of Dirichlet’s problem regular within the flow region and taking the value 
zero on the boundary appears in a similar role to that of a proper mode in 
the theory of vibrations 

This investigation is limited to plane transonie flow having a simple 
one-one hodograph representation. It seems to give a more precise 
formulation of the general notion of Frankl’s paper (6) inasmuch as 
potential flow for adjacent yet smooth boundaries is shown to be unstable. 
Here, however, non-uniqueness is proved, not, as in (6), the non- 
existence of solutions. It may be, however. that the type of breakdown 
suggested by Guderley (7) has also an exact formulation. The present 
method seems the most natural since only continuous perturbations of the 
boundary and of the velocity field are required to give instability, whereas 
Guderley considers small discontinuities in velocity. 


2. Basic equations for perturbations 
The equations of compressible potential flow with two independent 


variables can be written in the form 


A,+B, = 0, (2.1) 
Y »)»») 
C,+ D, = 9, (2.2) 
where w, v are functions of the variables (x, y), A, B, C, D are given functions 
of wu, v, and suffixes denote derivatives. For plane flow uw, v are components 
of velocity, and (2.1), (2.2) are equivalent to 
| —* @ 
v | (2.1.1) 
Bb u .) 
( sf — 
ty —, (2.2.1) 
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where ®, are the usual potential function and stream function respec- 
tively, while p = p(p) is pressure and p is the density. Pressure and 
density are determined as functions of the speed by 
| p-!'dp+4(u?+v?) = const. 

The forms (2.1) and (2.2) reappear in the equations of axially symmetric 
flow and, with a slight change of notation, the equations for one-dimensional 
non-linear waves. A linear combination of (2.1), (2.2) such as 

A,+B,+A(C,+ D,) = 0 (2.3) 
is said to be in characteristic form if after expansion it appears that both 
u and v are differentiated in the same direction, say, 


‘dy , 
. 2.4 
a] e 
The condition on A is clearly that 
B,+AD, B,+AD, ial 
pL = — = (2.5) 
A,tAC,  A,+AC, 
so that A satisfies the quadratic equation 
C CB 4 B 
. D i A D B)) AB 0. (2.6) 
u v | u v | u ? | u v 








| |denoting the Jacobian of the quantities enclosed. For plane flow the 
term in A disappears and the roots are given by 


A, A, p-'tane, (2.7) 


where ¢ is defined by c = qsine, 


c being the speed of sound and q the flow speed. Equation (2.4) therefore 
defines two families of curves, one of each family passing through every 
point of the (x,y) plane. To successive members of the system having 
dy/dx = fy let any monotonic set of numbers a be attached and similarly 
a set 8 to the other system. Then a, f are called characteristic coordinates 
and by the mode of construction are in one-one correspondence with the 
(x, y) plane. 

Let (u!,v!) and (u,v) be respectively solutions of equations (2.1), (2.2) 
and let du = ul—u, dv. = v'—v be everywhere small; let the corresponding 
variations of O(x, y). ‘V(x. y) be denoted by x, 7 respectively. Then, in the 
limiting case of infinitesimal velocity changes it follows from equatiors 
(2.1.1), (2.2.1) that 

A,, 8u+A,dv = x,y, B,,du+ B,dv = —x, 


C,,du+C, dv = Ny: D, du+ D,dv = —n, J 
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‘therefore, with dx, dy quite general 


dy (B, du+ B,.dv)dxa+(A, du+A,,dv) dy | (2.9) 
dy (D, du+ D, dv) dx + (C,, du+C, dv) dy 
However, solving equations (2.5) for A, we find 
>) (uA,,—B,) _ A, B,) (2.10) 
(uC, — D,) (uC,— D,) 


so that along a characteristic direction the ratio of the infinitesimal quanti- 
ties in equation (2.9) satisfies 
dy 


A, 2.11 
Zz, (2.11) 


the expressions in du, dv cancelling. The perturbation functions therefore 
satisfy the simultaneous partial differential equations 


Xa A, "Nex 0, (2.12) 
xetAs ng = 0, (2.13) 
and this result is quite general for all quasi-linear systems which can be 


written in the form of (2.1), (2.2). In particular, for plane flow it can be 
shown th it 


( coteq 1 dq \@ 2q9(B), (2.14) 
| coteq 1 dq 6 2f(a), (2.15) 


where f and g are arbitrary functions. If both f(a) and g(8) are monotonic 
in their arguments they may be identified with just a, 8 respectively. 
Such flows will be described as having a non-degenerate hodograph. The 


perturbation equations can then be shown to reduce to 


PIXat Ne = Y, (2.16) 
and WM"q (pq)Xe- (2.17) 


The derivations of equations (2.14)-(2.17) are routine and are omitted 
since equations (2.16) and (2.17) are to be derived by a method which 
makes use of the variation of the hodograph of the varied flow. It may 
be remarked here that «, 8 or qg, 6 refer to the flow before perturbation so 
that y, » retain their identity as perturbations in the physical plane. In 
particular, the conditions corresponding to a fixed boundary are ‘YY = 0 
and ‘1 = 0 so that 7 0 also. This is quite different from the normal 
hodograph method where the new boundary can only be found after a 
complicated integration. It is also noteworthy that with the special choice 
of variables (y, 7) and (q,@) the single system (2.16), (2.17) suffices for the 
whole class of flows having non-degenerate hodographs; du, 6v, however, 
involve the particular hodograph transformation explicitly. 
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3. Hodograph variations 
If !, ¥1 refer to the perturbed flow, regarded as depending on x, y, then 
mo! — M+ y, Yi= P+7 (3.1) 

and the velocities are given by 
d® (uda+v dy)) 
d¥ = p(vdx—udy) y° 


' may therefore be written as ®(q, #)+- x(q, @) and similarly for ‘V1, but this 
is not a hodograph solution since, on account of the perturbation, the 
velocities are changed according to 


em! 
— . " 9 9 
ui = qcos — = U—Xr (v.09) 
Ox 
em! ») 
vl = q'cos@! U—Xy- (0.4) 
Cy . 


Hence, with the understanding that y and 7 are infinitesimal, 
q' cos 6! = q cos 0—(x,9,+x0 9); 
q' sin 6* = q sin 6—(x,4,+ x9 9,)- 


Replacing the derivatives with respect to x, y by those with respect to q. @ 
and solving for dq, 50 


J pq dq Xq Yo+xo¥, (3.5) 
£7 66 Xa Dp Xe ®,, (3.6) 


where the Jacobian / is given by 


g ¥ = (Dy '¥,—®, ¥), (3.7) 
and (3.2) has been used to eliminate derivatives of (a, 7). 
Then (91, 6) (7, 6)-+ x(q, 8) 
O(q), 61) +-4(q!, 6) 
to the first order, where 
(q}, 6) (x—®, 5¢—®4 56)'. 
Using (3.5), (3.6), (3.7) it is found that 
d X—IXz: (3.8) 


By a similar calculation the perturbation of the hodograph stream function 
is found to be 


us N— PXo- (3.9) 
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[It is also possible to work in terms of derivatives of » and thus get 


d = x+p7'n, (3.10) 
hen 
( pt ‘ 
I $= y—-),,. (3.11) 
(PV)q 
omparison of (3.8) 9) with (3.10), (3.11) gives 
9 ner PIX = U. (3.12) 
qq (PT) Xe: (3.13) 
this that is, equations (2.16), (2.17) already given above. It will be noted that 
the and » are the Legendre transforms associated with the hodograph per- 
turbations 4, % respectively (cf. (3), pp. 248-52). A useful check on these 
calculations is afforded by the verification that (taking account of these 
last two equations) ¢ and y& satisfy the usual hodograph equations, viz. 
p*9¢,44 (pq), %e 0, (3.14) 
4) | 
pb Qh: (3. i 5) 
| Equation (3.11) reduces the calculation of 7 to the mere integration of 
standard hodograph solutions according to 
| ¥ i | ‘ > 
pq, | &p-*q-*(pq)_4q}. (3.16) 
By means of the relations of this section, the variation of the Jacobian 
o furs , (PIa wr? — 
5 p=, 2 | V2 + AE ; (3.17) 
\ py 


on the fixed boundary of a flow can be conveniently found. A condition on 
the perturbation is that » = ©, for this boundary remains a streamline 
while the curvature « also remains constant. It can be easily shown 
59 ef. (5)) that ; : a ae 
Kk = —Y/(pgf). (3.18) 
Hence aur R 2 
KOS 8(Y / pq) —8(Dp_/q") 


s 1/71? =. 2 
gil 9 )— 59 (Pl |. 


But 6@ vanishes on the boundary so that 


F } {4 \1 fA2 \1 ‘ 1 
(1) ‘D od! o +(; >) 6g + op . 
e6 oé} : 06 oqee ' \00 


due allowance having been made for both the change of the form of the 
hodograph of the flow and of the image of the point on the boundary under 
consideration. Again from (3.5), (3.6) 





84 = 9Xx4/9,. (3.19) 
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Hence, dropping the accents, 
Kd (q-32D_—q 20 ,0)9Xq ®, —q~*dp, 
which may be conveniently written 
89/F =- (7 *09) — (7X2) 1 by/ Dp. (3.20) 
q 2Dp 0, 

This result shows that if regular perturbations of the velocities are intro- 
duced into a field the change of Y remain small.t Even if % actually 
vanishes on the boundary, in the vicinity of the zero the changes of ¥ in 
the hodograph of the perturbed flow remain proportionately small, and so 
flows with isolated zeros of Y are perturbed into similar ones. 

This section therefore adds to the existing evidence that the vanishing 
of the Jacobian does not explain the breakdown of potential flow. 


4. Perturbation theory for the vortex 

A very simple example of a mixed flow is furnished by the vortex which 
was studied by Guderley (7) by some rather involved calculations in the 
hodograph plane. 

[f we suppose that 0 = (6), VY = V(q), then equations (3.14), (3.15) 


give a solution 


di q 15 dq. (4.2) 

Using equation (3.2), the streamlines are, with a suitable choice of origin, 
x —q-'sin 6, (4.3) 

y = q-'cos8. (4.4) 


The boundary of the flow, which extends to infinity, may be taken to be 
the image of any supersonic circle, say q = q, > Cx, With the usual notation 
for critical speed. Let % = > (a, cosn6+b, sin nO)x,,(q). where the y,, are 
the quantities considered in (8) and so regular at q = 0, which is the 

image of the point at infinity. Then a general perturbation solution is 
n = > (a, cosnO+b, sin n)n, (q), (4.5) 

where the 7,,(q) are derived from the y,,(¢) by means of equation (3.16). 

The perturbation of the boundary is then 

> (a, cos n6+b,, sin n6)n,,(q,) = ¥ (A, cosn6+ B, sin n6). (4.6) 


The A, and B,, being specified, the perturbation of the flow is given by 


nr 


tn Sita Ad EN (4.7) 


<* (A,, cos n6+ B,, sin n6)n,,(q) 
‘ +: cos) And 
nl) 


+ The difficulty due to the denominator vanishing at the sonic line is only apparent. 
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Therefore if g, is nearly equal to a zero of some 7, (q), the term 
11 . 1 In 7 


A,, cos n6+ B,, sin n@ 


gives a correspondingly large change of the flow for a small perturbation 
of the boundary of this form. If q, is equal to a zero of some 7,(q), the 
solution (4.7) has no meaning but the term (A, cos n@+ B, sin n@)n,,(q) has 
the property of vanishing on the whole boundary of the flow, being regular 
within, and yet not vanishing identically. Now it is easily shown by 
evaluating the »,,(q) in terms of the ‘,(q¢) of (8) that the first few do have 
roots in the supersonic region (see Appendix). It is also true that for high 
the m vanish very near the sonic circle. This may be proved by integrating 
the asymptotic formulae (8) for % from the origin up to the first loop of &,,(q) 
in the supersonic region. It is also an immediate consequence of the oscilla- 


tory character of the asymptotic expansions of the corresponding functions 


y,(q). Thus if » g(q)cos n6, x = f(q)sinné equations (3.12). (3.13) give 
ng = pqf’, (4.8) 

n(pq),f = 99: (4.9) 

ind f satisfies (pqf’)’ = n®q-(pq),f. (4.10) 


[tis then a well-known result that for large n. f oscillates rapidly in the whole 
region q > ¢c,, the term p~'(pq), |—q?/c* being negative in that region. 
Therefore, by (4.8), g vanishes frequently in the supersonic region. 

[t follows that, for many speeds of flow on the boundary of the cylinder, 
resonance can occur between certain perturbations of the boundary and 
the flow. Such perturbations give rise to proportionately large changes of 
flow at neighbouring speeds, while at these speeds the functions 7(q, @) may 
be added as desired. Dirichlet’s problem is not unique for a whole sequence 


of boundary speeds q, and, since, for high n, the zeros of 7,,(q) lie close 


together, in a practical sense is never unique. Even for low n < 7 there are 
four roots in the range 0-24 < + < 0-30 where 7 = q?/q?,,, = } corresponds 


to sonic conditions (see \ ppendix ). 

The success of the above investigation into instability of the flow around 
a circular cylinder depends largely on the practicability of constructing 
general perturbation solutions at all flow speeds. Only a single infinity of 
speeds furnish examples of resonance, and of the non-uniqueness of Dirich- 
let’s problem for the perturbation stream function. It seems likely, al- 
though this fact is not proved, that in other cases the perturbation functions 
are unique. To prove the existence of unstable flows it is not necessary to 
construct all the hodograph solutions for a given body as the stream speed, 


or some equivalent parameter, is varied. This was done for the circular 
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cylinder, but is clearly a most difficult matter in general. It is sufficient to 
show that there is both a hodograph solution and a regular perturbation 
function having respectively ‘Y’ = 0, » = 0 on the same (hodograph) 
boundary streamline and that the hodograph corresponds to a true physical 
flow. This can be readily proved for certain families of boundaries differing 
slightly from the circle, the behaviour of the flow at infinity still being that 
of the vortex. 


Let n He Sin 26+-en, sin 46 
be a perturbation function, where ¢ is small, and let the boundary of the flow 
be approximately q xy, Where a, is the smallest root of ,(q) = 0 lying 
in the supersonic range of q. Then to a first approximation the boundary 
hodograph is given by 
q v9 -+ Aq(6), 


where Aq — 2€7,4(a.)cos 20 nol Xo) 


and (see Appendix) 73(a,) and 7,(a,) do not vanish. Then the hodograph 
solution having ‘’ constant to order e on the same boundary is given by 


_ [ pq dq+-e’¥3(q)cos 26, 
provided that €’ = 2ep(xy)4(x2)/{ag Vo(ae)y9(a9)}; 
it being noted (reference (8)) that ‘(a.) + 0. 
A similar family of hodograph solutions and corresponding perturbation 
functions is given by 
7 = N3(q)sin 30+- en %(q)sin 68, 
Y= | pq? dq+e"¥;(q)cos 38, 


the ratio of e’ to « being suitably chosen. 

To complete the proof of the existence of unstable flows it only remains 
to be shown that the hodograph solutions Y can in such cases be mapped 
into physical space without limit lines. In reference (9) the author has 
shown that if any analytic hodograph solution is given with the boundary 
not touching a characteristic in the hodograph plane, then none of the 
streamlines inside the flow can do so and the whole solution corresponds 
to a physical flow. With a slight, but obvious change, this result applies 
in the present case since the boundary streamline is entirely supersonic 
although the flow is of mixed type. It is also clear that the hodograph of 
the flow corresponding to the almost circular cylinder differs only slightly 
from that for the circle and in the latter every streamline cuts the hodo- 
graph characteristics at a finite angle. It has therefore been shown that 
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the resonance effects occur for nearly circular boundaries as well as for 
the exact solutions for the circle. While the analysis may well be difficult, 
it seems likely that all families of transonic flows will contain examples 
of resonance, and of the non-uniqueness of Dirichlet’s problem, provided 
that the parameters on which they depend can be varied sufficiently. 

It conclusion, it may be noted that the above discussion of plane tran- 
sonic compressible flow is carried out entirely within the perfect fluid theory. 
It has not been necessary to make any allowance for boundary-layer effects; 
nor have time-dependent disturbances been introduced to prove instability. 
This phenomenon has been shown to arise as a purely mathematical fact 
associated with the non-uniqueness of Dirichlet’s problem for the perturba- 
tion of the stream function in the physical plane at certain critical speeds. 
[t is also noteworthy that, for all supersonic boundary conditions on the 
cylinder, there are many critical speeds, and that the corresponding proper 
perturbations of the boundary are rapidly oscillating functions. In physical 
terms, roughness of the boundary should give rise to instability for super- 
sonic conditions. This is known to be true in practice. There are, it seems, 
good grounds for expecting the complete explanation of the instability of 
plane transonic compressible flows to be found in a study of the, admittedly 
intractable, problem of the uniqueness of Dirichlet’s problem for the linear 
differential equation satisfied by the perturbation stream function. 

The author gratefully acknowledges that he was first led to apply 
variational methods to the problem of the instability of transonic flow by 
1 suggestion of Professor W. R. Sears. 


APPENDIX 


The perturbation stream function 7 can, according to equation (3.16), be expressed 
terms of the hodograph stream function of ref. (8). Thus 
T 1—7)*5705 | Yi(7)(1—6r)(1—7)-857-85 dr 
0 
I x tor ? 
T) 4(1—7)?*5705 

d for ?  ¢ 1 /7°> and the lower limit is now 7 = }. The 7, have been 
tabulated for n | to 7 and for n 1 by numerical integration of the tables of 
F(z) of ref. (8) and for 7 ? @2,,. from 0 to 0-40 at intervals of 0-02. It is hoped 
to prepare tulle! tables at a later date. 

For convenience of presentation the functions are given as multiples of y,(%) and 
the value of this is given below each column except that 7_,(4) can be chosen arbi- 
trarily and is here taken as zer 
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Table of Reduced Function »,,(7)/,(4) 
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\ FOURIER INTEGRAL SOLUTION FOR THE 
STRESSES IN A SEMI-INFINITE STRIP 
By L. H. MITCHELL 
(Dept. of Civil Engineering, The University, Bristol)t 
ceceived 5 June 1952] 


SUMMARY 


\ solution is given for a stress diffusion problem in a semi-infinite strip. Fourier 
integrals are used and the method is applicable to problems where the boundary 
onditions are given in terms of stresses. Good agreement is found between the 
theoretical solution and a photoelastic investigation. 


Introduction 
In this paper Fourier integrals are used to determine the stresses in a semi- 
infinite strip, the boundary conditions being given in terms of stresses. Both 
Howland (1) and Hopkins (2) have used Fourier integrals for the solution 
of generalized plane stress problems in infinite strips with a variety of 
boundary conditions. Howland’s paper is confined to cases where the 
stresses along the long edges are prescribed, whereas Hopkins treats the 
case of mixed boundary conditions. 

By adding further terms to the Fourier integral stress function given by 
Howland, it is shown here that the boundary conditions of a semi-infinite 
strip can be satisfied if they are expressed in terms of stresses. The particular 


problem solved is one of stress diffusion and is illustrated in Fig. 1. 




















4 y 
r 
~<~—- [ q 

I Thickness t 

0 x 
i —_ | Area A 
<= (—_ ee 

Fic. | 
Notation 


A Area of one stiffener. 
b Half-width of strip, taken as the unit of length. 
B,C, D Functions of m. 
m A parameter. 
P Concentrated load applied to stiffener. 
+ Now at the Aeronautical Research Laboratories, Department of Supply, Melbourne. 
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t Plate thickness. 
(x,y) Rectangular coordinates. 
A Dimensionless parameter A/bt = A/t. 
Direct stresses in the x, y directions. 
o, Direct stress at the end of the stiffeners given by P/A. 
7 Shear stress. 
@(x,y) Airy stress function. 


For convenience the unit of length is taken as the half-width of the strip, 


Details of the solution 

The stresses are to be found in the plate shown in the above figure, which 
also shows the position of the rectangular axes. The coordinates (x, y) which 
have been chosen are dimensionless and are obtained by dividing the true 
distances by b, the half-width of the strip. In effect 6 is taken as the unit 
dimension and consequently all quantities involving length must be modi- 
fied accordingly; for example the area of the stiffeners must be divided by > 
and the result is denoted by A. The stiffeners, of the same material as the 
plate, are assumed to resist longitudinal direct stresses only and it is also 
assumed that this direct stress is constant over any cross-section of a 
stiffener. Thus the stiffeners have been idealized since this assumption 
implies infinite rigidity with respect to shear stresses only. The short side 
of the strip is required to be stress free. 

By means of the Airy stress function 

(B cosh my+ Cy sinh my)sin ma dm, 
0 

which can be shown to satisfy the two-dimensional bi-harmonic equation 
of generalized plane stress, it follows from Howland’s paper that the 
boundary conditions along the two sides of a semi-infinite strip, sym- 
metrically loaded, can be satisfied. However, with this stress function both 
the direct stresses co, and co, are automatically zero along the short side of 
the strip and therefore additional stress functions are required to satisfy 
the stress-free condition (i.e. ¢, = 7 = 0) along this edge. The true trans- 
verse strain distribution at the end can be represented by a cosine Fourier 
series, and the suggested additional stress function ist 


> D, F,(x)cos dnzy. 
n 
It is required that F,(0) = 0 to enable the boundary conditions o, = 0 at 


+ The general form of F,,(x) is (v+H,)e—4""* enabling the two boundary conditions at 
x = 0 to be satisfied. 
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r= 0 to be satisfied and also each term must satisfy the equation of 


generalized plane stress. Such a series is 


LDa?+ S D, xe-*"** cos inzy, 


vt 
where ? |, 3, 5,.... The first term permits o,. obtained from this series, to 
sume other than a non-zero value at y +1. Thus the stress function 
enabling all the boundary conditions to be satisfied is 
D > B. cosh my -+-¢ Yy sinh my )sin max dm+ 


| 1 D, x? | N D, re hnaex cos Try. 


n=1 

The boundary conditions to be satisfied are 

| ()} 4 () 

2) + () 0) 

4 O. y | 

1) + j 0. y | 

’ 

The last condition (4), obtained from consideration of the equilibrium of an 


element of length of the stiffener, does not specify the load applied to the 
strip. By replacing this condition by the condition that the total load 
ried by any cross-section of the strip (including the stiffeners) is constant 
ind equal to the applied load, this difficulty is overcome. It can be shown 
that this latter condition implies the condition (4). With the restriction 
that the terms with suffix other than zero add nothing to the load carried 


by the strip, the conditions (3) and (4) lead to the relationst 
37 By 


D,/m3)(sinh m+-m cosh m+-2mA cosh m-+-Am? sinh m)—(P/tm)sinh m 





m-—+-sinh m cosh m+ 2Am cosh?2m 


(D,/m®)(m sinh m+-Am? cosh m)+-(P/tm)cosh m 





m-—+-sinh m cosh m+ 2Am cosh?m 


The shear stress due to these terms involving B,, Cy, and D, is given by 


rt) 


m(mB, sinh my+C, sinh my +C, my cosh my)cos ma dm. 


f 
* sin mx , 
Obtained with the use of the identity dm $7. 
nt 

















54 L. H. MITCHELL 
For large m and y + 0 we may neglect me-?” and write 
sinh my = cosh my = de”; 
the integrand of this integral then becomes 
_2 [| D(1+Am)—mP/t|(l—y)+AD,+P t 
TT (1+2Am) 


»>—(1l—y)m COs Mx. 





[f the integral is to converge at the point (0,1) then D, = — P/At and the 
integrand for large m approaches 

2P 
th 


(l—y)e-2-™™ cos mx. 


| 


Consequently the shear stresses can be written 


, 
rt 2P 
a | F(m, y) += + (1—y)e“2-™ | cos ma dm, 
Be | a tr 
0 
where F(m,y) > 0 exponentially as m + oo for 0 < y < 1, ie. 
= 2 P l 2 
, _2 (l—y) 
T | F(m, y)cos mz dm+— = -—————, 


a tAx?+(1—y)? 


0 
Hence, there is a discontinuity at the point (0,1) such that 


x. 


lim 7(a, 1) = [ F(m,1) dm 
xr—>0 0 


x 
- 


and lim 7(0,y) = | F(m,1)dm+ 
y—1 


P 
th 


Yio 


0 
Thus, when the end of the strip is free from shear, the peak shear stress 
along the edge stiffener is 2P/ztA which is 2/7 times the direct stress in the 
stiffener. However, this type of discontinuity can be shown to involve a 
rotation of the corner through 90° and for this reason the analytical solution 
cannot represent the stresses in small regions near the corners of an actual 
plate. 

There remains the problem of eliminating the finite shear stress left 
acting over the end. This can be done with the aid of the remaining terms 
of the stress function. The boundary conditions determining B 
D,, are 


, C,,, and 


n 


(2) ow, = 0, y= +1; 


(3) The load carried by the plate and stiffeners at any section is zero. 
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The conditions (2) and (3) lead to the relations 


a7n?m D, sinh m 





B = Batutthecnes Jewtee ne 
. (m?+- 4n?7?)?(m--sinh m cosh m--2mA cosh?2m 
} 
, 7n? D, m cosh m 
m?+-1n222)2(m-+-sinh m cosh m+ 2mA cosh2m) 
} 
Thus the expression for the shear stress at x 0 can be written in the form 


7= > D,Anzrsinjnry+ >Y D, | f,(m.y)dm+P | fi(m.y) dm = 0, 
n=1 © . 


n=] 0 0 


where the infinite integrals are uniformly convergent for |y| < 1. It is 
therefore possible to find the coefficients D, by normal Fourier series 


methods and thus completely determine the desired stress function. 


Convergence of the infinite integrals 

If the above stress function is examined it will be found that certain of 
the integrals do not converge at the lower limit. This is due to the fact that 
the stress function must approximate to {P/(1+-A)t}}a? for large values of 
z and that a Fourier integral cannot represent such behaviour. The non- 
convergent part can be represented by 


F(m.y 


mi” 


VM . 
sin ma dm, 


0 


where F(0,y) 4 0. If this integral is replaced by 


m? m° 


¥ Fi .. 2 F(0 W) F . Wi = 
| mY) _ J) ‘sin me dm+ F(0,y) | tan“; dixdx, 


where? is an arbitrary positive constant, then the modified form of ® is such 
that V¢‘® = 0 and it is easily shown that the boundary conditions are 
satisfied. As the two expressions are identical after differentiating twice, 
both the original and modified forms of the stress function lead to identical 


values for the stresses. 


Example 

For the case where A = | the shear stresses along the edge of the strip 
shown in Fig. 1 have been calculated. Only the terms involving D, and D, 
were used to reduce the end shear stresses but they were not found in the 
manner indicated above. These coefficients were determined by making the 
shear zero at the points (0, 1) and (0,4) as a Fourier series would only con- 
verge slowly near the point y = 1. The expressions for the stresses were 
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checked to some extent by finding their limiting values at x = 0 and 
infinity. These were found by the use of the relations 

: ge sin , . 

lim | k (m)| |ma dm 37 lim m F'(m) 


ee a cos m—>0 
0 


x 


and lim F(m)sin ma dm tz lim m F(m). 


r—>0 ¢ —>o 
0 m 


The first of these theorems is well known. The second can be proved if 
m F(m) is continuous and also monotonic for large m. Both integrals must, 
of course, be convergent. 

The infinite integrals were evaluated by means of the quadrature sug- 
gested by Filon (3) and also used by Howland (1). The result of these 
calculations is shown in Fig. 2 where a comparison is made with the shear 
stresses obtained from a photoelastic investigation carried out by the author, 
The residual shear stress left acting over the end is also shown. It is seen 
that the agreement obtained is good. 


*t 


* 





> 
G 
x and y 
Residual end shear stress 


Fic. 2. Comparison of the theoretical and experimental 


shear stress distribution along y a. 
— theoretical, -—— experimental. 


The author wishes to thank the Chief Scientist, Department of Supply, 
Australia, for permission to publish this paper. 
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THE DIFFRACTION AND REFLECTION OF SHOCK 
WAVES 
By W. CHESTER (Dept. of Mathematics, The University, Bristol) 
Received 17 July 1952; revised 4 December 1952] 


SUMMARY 


The pressure field produced behind a plane shock wave of arbitrary strength 
when it encounters an infinite yawed wedge at incidence is investigated. A linearized 
solution is obtained on the assumption that the shock meets the boundary at near 
glancing incidence. The solution is used to describe the pressure distribution on 
the wedge for different values of the pressure ratio across the shock front, and 
various angles of yaw up 1 certain limiting value. 

\ more accurate investigation is made of the boundary of the perturbed region, 


und the strength of the non-uniform shock which occurs there is calculated. 

1. Introduction 

A THEORETICAL investigation of the flow field produced by the interaction 
of a plane shock wave and a rigid obstacle is difficult, not only because of 
the non-linear nature of the problem, but also because of the fact that the 
non-uniform shock waves which occur imply variations in the entropy of 
the fluid and a loss of the irrotational character of the motion. To simplify 
the problem sufficiently for a theoretical attack to be successful, two 
possible courses are 0, en. The first is to linearize the basic equations of 
motion on the assumption that the incident shock is weak. Since the 
entropy change across a weak shock is of the third order in the shock 
strength, the entropy variations are also effectively eliminated and, in fact, 
the problem belongs to the theory of acoustics. This method was used by 
Sommerfeld (1) to solve the problem of diffraction by a semi-infinite plane, 
and by many later writers. 

The second possibility is to consider an incident plane shock of arbitrary 
strength and linearize the basic equations on the assumption that the 
obstacle produces only small perturbations in the uniform flow behind the 
shock. Although entropy variations are no longer negligible, it appears 
that they can be suppressed from the theoretical investigation by basing 
the analysis on the pressure variations in the fluid. This method was first 
used successfully by Lighthill (2) to solve a two-dimensional problem in 
diffraction, and later by Ting and Ludloff (3), who solved a similar prob- 
lem by a different technique. 

The present investigation belongs to this second category. It is an 
extension of Lighthill's theory (2) and considers the interaction of a plane 
shock wave with an infinite thin wedge at small angles of incidence relative 
to the uniform flow behind the shock, and arbitrary angles of yaw up to 
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58 W. CHESTER 
a certain limiting value depending on the shock strength. Suppose the 
velocity of the shock is U and the wedge is yawed through an angle 8 
(see Fig. 1). Then the point of intersection of the shock front and the 
leading edge of the wedge travels along the leading edge with a velocity 
U/sinB. Let an equal and opposite velocity be superimposed on the whole 
field. If a supersonic velocity of the fluid behind the now stationary shock 
results, we then have a problem of steady flow behind the shock which is 
in many respects similar to the well-known cone-field problem of Buse- 


mann (4). 








Fic. 1. Configuration in the (x’, 2’) plane. 

If the component velocity of the fluid normal to the leading edge is 
subsonic, the region of perturbed flow is bounded by the shock and the 
Mach cone with vertex at the junction of the shock and the leading edge. 
But if the fluid velocity normal to the leading edge is supersonic, the cone 
intersects the wedge and the boundary must be augmented by the tangent 
plane from the leading edge to the cone. This latter boundary will be a 
shock wave or Prandtl—Meyer expansion, and behind it the flow will be 
uniform and parallel to the wedge, at least until the Mach cone is reached. 

The simplifying feature of the problem arises from the lack of a funda- 
mental length in the data defining the problem, implying that the flow 
variables must be constant along straight lines from the vertex of the Mach 
cone. The novel feature is the extra boundary, in the form of the incident 
shock, on which certain boundary conditions must be satisfied. It will be 
shown that the governing equation for the pressure can be reduced to 
Laplace’s equation, and the solution is then obtained by conformal trans- 
formation. 


2. Notation and basic equations 

Behind the shock there will be a region of uniform flow which is not 
affected by the presence of the wedge. In this region the fluid velocity, 
pressure, density, and sonic velocity are denoted by 4q;, p,, p,, and 4 
respectively; ahead of the shock the corresponding quantities are 4p, 9, po: 








and ¢ 


y= 


qi = 

Le 
be tl 
equa 


As 
the 1 
weds 


unife 


This 


whit 


The 
and 
The 
(B < 

H 
intr 
bow 
of t] 
of tl 


shoc 


The 


att 
in t 


101e 
L 
OCK 


h is 








THE DIFFRACTION AND REFLECTION OF SHOCK WAVES 59 


and a). If U is the velocity of the shock, then, taking the adiabatic index 
1-4, the shock transition relations give 


9 


4} BU ay u~*), P} BP U? 145), Py 6H py (1 T 5a6 U “). (1) 
Let M U/a, (> 1) be the Mach number of the shock, and let M, = q,/a, 
be the Mach number of the uniform flow behind the shock. Then, from 
equations (1), 5(M2—1 
M, xn. (2) 
(71 +5) 

As suggested in the previous section, it is convenient to superimpose on 
the whole configuration a velocity U/sinB along the leading edge of the 
wedge so that the shock is reduced to rest. The velocity in the region of 
uniform flow behind the shock is then V,, say, where 

V3 U?/sin*B+-q?—2Uq,. (3) 
This velocity is supersonic provided that 
U?/sin2B 4 qi 2U4q, > az, 


which is equivalent to 





6.M4 
tan®8 << ——. iM a—— g (4) 
(M?—1)(M?-+-5) 
The right-hand side of this last inequality is singular for M = 1 (B < }z) 
and decreases rapidly to a minimum value of 42 (8 < 61°-28) for M? = 3. 


Thereafter it increases monotonically to an asymptotic value of 6 
(8 < 67°-80) as M + oo. 

Henceforth we assume that condition (4) is satisfied. The perturbations 
introduced by the presence of the wedge are then confined to the region 
bounded by the shock front and the Mach cone with vertex at the junction 
of the shock and the wedge leading edge (the point O of Fig. 1). The axis 
of the Mach cone is in the direction of V,, and subtends an angle pu with the 
shock, where 


(M?-+-5)tan B 








tan pu (l1—q,U-)tan B = +— 5 
The semi-angle of the Mach cone is a, where 
a, a, Sin pu (7TM?—1\} . - 
SIN « —— — = | ———__| "sin w. 6 
V; l V1 tong f (6) 


Let (x’, y’,z’) be rectangular Cartesian coordinates in the field with origin 
at the point O and the z’-axis in the direction of V,. If the flow quantities 
in the region of perturbed flow behind the shock are designated by the 
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suffix 2, the equations of conservation of mass, momentum, and entropy 


are then V,.Vp, = —p.V- Ve 


] 
(V,.V)V, : —>VP2 } (7) 
2 


V.. V8. 0 
where s, denotes the entropy. 

The wedge is assumed to lie approximately in the plane y’ = 0; equations 
(7) may then be linearized if the flow variables are assumed to differ by 
small quantities from their values in the region of uniform flow behind the 
shock. If V, = (ug, v2,y,+w,), we then have 





, OPs {CU , Wy , We \ 
AS = -a\—-4+-—5+ 
ez \aa’ ' dy’ ' ez 
- (d) 
, Us l Cp, , OV, l Cp, , OW» 1 Op, O85 
A=; ——35; ho=-- ss ASH Ss G0 
Oz py Cx C2 Py CY Oz py C2 Oz 4 


The left-hand side of the first of these equations, by virtue of the last 
equation, may also be written V, ép,/a? éz’, since the density is a function 
of the pressure and the entropy. 

One further simplification is possible: since there is no fundamental 
length in the data defining the problem, the flow variables are expressible 
as functions of 2’ /z’, y’/z’ only. We thus make the following transformations: 


, , 
x 2 
tape, | ee (9) 
z tana z tana 

lo—/p a u. v. V—w, 

p = Po Pi —1(ps—p,), w: Ra, Pi meen, Dimi wood, 
Pit P11 q, COS q1 COS x 9, SIN a 
(10) 

Equations (8) then yield 

ou cu cp ov ov Cp C o | ou ow 
xr—ty—-: —s x~—+y- —, t—+Y—])|p,p,w = —+-—. 
CXL cy Cx Cx cy cy Ox CY Cx OY 
(11) 


Apart from the equation for w, these are formally identical with the 
equations obtained by Lighthill in the two-dimensional problem of diffrac- 
tion by an edge (L. 8).+ Also, the independent variables defined by (9) are 
easily seen to reduce to (L. 7) if we take a fixed set of coordinates (X, Y, Z) 
for which 


faves: rey, see Loe, 


sin B 7 sin B 








+ This denotes equation (8) of reference (2). 
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where, for the purpose of this transformation, the origin O is considered to 
be moving along the leading edge with velocity U/sin 8. 

We now have 

1 = g,t y" Y ‘ 
+e y= = > (13) 
tan « a,t z tana a,t 


» t) 


We note that the last of equations (11) imply that the difference between 


(a (b) 


Fic. 2 figuration in the (2, y) plane. 


Subsonic flow behind the shock. 
Supersonic flow behind the shock. 


any two of p, p, or w is a function of y/x only. Since these quantities are 
zero on the Mach cone, it follows that they are all equal in the region 
bounded by the Mach cone and the two planes extending from the z’-axis 
to the lines of intersection of the Mach cone with the shock front. The 
entropy variations are also zero in this region, for to the present degree of 
uccuracy 8, is zero on the Mach cone and it is easily deduced from the last 
of equations (8) that s, is a function only of y/x. The motion is thus effec- 
tively irrotational, the reason being that the fluid in this region has not 
crossed the perturbed portion of the incident shock, so that only through 
crossing the Mach cone could entropy variations ensue. Since, as will 
appear later, this boundary can be a shock whose strength is at most of the 
second order, it is not surprising that the entropy variations are negligible 
on a linear theory. 

It is in the remaining portion of the perturbed region, which has the 
incident shock as part of its boundary, that the rotational character of the 
flow becomes evident. The variation of p or w in this region is deducible 
from a knowledge of their values at the shock and a knowledge of the 
pressure field. 


From the pressure one can also calculate the two remaining components 
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of velocity, u and v, from their values on the boundary. This is done by 
integrating along radii vectores in the (x,y)-plane, starting from the 
boundary of the perturbed region and using the first two of equations (11), 


3. The boundary conditions on the shock 
The undisturbed part of the shock front lies in the plane 
2’ /2’ = tanp 


tan pu 


or (14) 


tan a 

Consequently we take the equation of the perturbed shock front to be 
xtana = tanpu+/f(y)sec p 
y’ 
or x’ cos p—2’ sin p— 2 {— = 0, (15) 
z tana 

where f is assumed to be uniformly small. 
The direction cosines of the normal to the shock front are then propor- 


tional to COS LL, f' cota, (sinu+f—yf’), (16) 


and, if n is a unit vector normal to the shock, we have 

n = {cosu — sinw cos u(f—yf’), —f’ cot a, —sinw — cos*u(f—yf"’)} (17) 
to the first order. 

Now the shock transition relations will depend on the normal component 
of the velocity in front of the shock; in fact the relations for a stationary 
shock may be written 


/ 


9 
: az 
l aa 0 5 fy2 1,2) \ 
V.—Vo Vu Ts i) Po = Spi Vin—j45}- (18) 
n 


where V, is the velocity in front of the shock, in the present case U/sin§ 
parallel to the wedge leading edge, and V,, is the component of this velocity 
normal to the shock. 
Since V, = {—U sin(8—p)/sin 8, 0, U cos(8—p)/sin B} and V, = (V,.n)n 
it follows, with the help of (17), that 
7. U{1+-cot B cos n( f—yf’)\n. (19) 
Equations (17), (18), and (19) then yield the following relations: 


5 U cos p 


u a ~ {eos(B-+-)+.M-* cos(B—p)(f—yf’) 
6 g, sin cos pu , 
v f'/sin« 
5 Ucosp , 2 : " : (20) 
Etna a | 
. 6 q, sin Bsin «| mete in(B—p)i(f—uf) 
p A. cot B cos p(f—yf’). 


34, PN 
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Thus, on the shock, we must have 





u = Ap w=’, yg = BS, (21) 
cy cy 
where Ay py cos(B- f)+ M 2 cos(B—) 
2p, l cos a cos B 
| /7.2—1\3cos(8+p)+ M-? cos(B—p) (22) 
2\ M2+5 cos a cos B ; 
QA, py sin(B- ft) M 2 sin(B }t) 
, 2 po f sin a cos B 
| 7M? ‘\* sin(B-+-p) M 2 sin(B— H) 23) 
3| M?+-5 | gina ecos B ‘ = 
R 3 1p; % tan p 3 py Ni secty 
5 py U* sin a cos p 5 Po 
M?—1 
3 'E -}sec*m. (24) 


Equations (22) and (24) reduce to (L. 15) when B 0. 
In conjunction with the basic equations (11), the first and last of equations 
21) yield the following condition on p, to be satisfied when x k (L. 19), 


cp jcep (A + k)y Bky 1 
ca éy | ke ‘ 


(25) 

If the upper surface of the wedge is inclined at an angle 6 to the plane 
y = 0 (8 is defined as a rotation about a line in the plane y = 0 parallel to 
the shock), then the boundary condition on the wedge is v = dseca. In 
particular this is true at the shock, and so we must also have 


B 
dy dp dSeC a, (26) 
CY Y 


. 


? 


where the integral is taken along the shock from the wall to the Mach cone. 


4. The boundary condition on y = 0 

We assume first that the leading edge of the wedge lies outside the 
Mach cone, the condition for which is a < (B pL) or, by (6), 
(7M?—1\3 . : 

7 }*sin p < sin(B LL). 


With the help of equation (5), this gives 
5(M2—1)ecos B 
M?4+-5)(7M?—1)} 


or M, cos B |, as one would expect. 
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In this case the flow patterns above and below the wedge are independent, 
and for the upper half-plane the boundary condition on y = Oisv = 8 seca, 
Since ¢p/ey = 2 év/ex on y = 0, by (11), this also implies ép/éy = 0 on y = 0, 

In the subsonic case (/, cos 8 < 1) it is more convenient to separate the 
symmetrical and antisymmetrical contributions to the inclination, say §, 
and 6, (again these are the contributions to the effective inclination of the 
upper surface relative to the flow behind the shock so that, for example, 
the semi-angle of the wedge would be 5,sec8). We then have two con- 
tributions to the pressure; the first, p,, is symmetrical about y = 0 and 


so cp,/cy = 0 on y = 0 remains true. However, at the leading edge, 


m « m c 
: ° Op rs Ov. 
lim Ps dx — —m Sdx = —mé, sec «, 27) 
y>o J cy J & 
mn c m c 
tan(B L 
where m an(P—#) (28) 
tan a 


and —m is the x-coordinate of the leading edge. 


The second contribution, p,, is antisymmetrical about y = 0, and con- 


tinuous for —1 < 2- m. Hence p, = 0 for l<2x< —m, and 
op,/ey = 0 for —m < x < k. Furthermore, 
-m+C — + C 
: - 1ép ; an) : 
lim Pa dy — lim —“dzx = 6,8eC a. (29) 
y—0-+ J ZC y y—0+4 “ Cx 
1 “1 


5. The boundary condition on the Mach cone 

In the subsonic case the boundary condition on the Mach cone is p = 0, 
but in the supersonic case this is true only from the shock to the point of 
contact of the tangent plane from the wedge leading edge. The latter will 
represent a weak shock wave (or Prandtl—Meyer expansion if the inclination 
is negative). Between this weak shock and the Mach cone the flow is 
uniform and the pressure is, to a linear approximation, 


7: dcosB 
Io = )—-- o—SSSF 
Po = Piv Pr Mi (Mj cos*B—1)! 
5.M, cos 8 m 





or p = . d5 sec «x. (30) 


(M3? cos*8—1)! — (m?—1) 
The last relation follows from the fact that 
M, cos - sin(B—p) 
sin « 
deducible from equations (5) and (6). 
Thus on the Mach cone p is zero only for —m-! < x < k; in the interval: 


m<x < —m-' it is given by (30). 
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6. The solution in terms of the pressure 


[he boundary conditions required are now all expressed in terms of the 
essure. Also it may be deduced from equations (11) that (L. 16) 
| c Cp Cp 
v2, y= + 1)(ee ty (31) 
cy }\ @2 cy 
\\ nstorm: 
» 
Z 2p ; an 
cos 6 y sind (02) 
] o- l p- 
juation (31) reduces to Laplace’s equation (L. 22) with (p,@) as polar 
rdinat lhe Mach cone becomes that part of the circle p l upstream 
the incident shock, the latter being an are of the circle 2p cos 6 k(1-+ p?) 
orthogor it cos 6 = k in the (p, @)-plane. 
ock-boundary condition (28) becomes (L. 28) 


Aktan 6— B cot 6 


oa (33) 
(1—k* sec?@) 
here nd dn are positive tangential and outward normal elements to 
The leading edge of the wedge becomes the point 
0 p ‘1 (1 m*)?! m 
1 conditions (27) and (29) become, respectively, 
! p- cp, 
lim ’ ls dp m0, SEC x 
ed - l p~ po 
" Op m ' 
Or lim | - dp = — ——_—— 9, seC a (34) 
9\1-s 
bor J pou (1 m*)? 
" r l p” cp - 
und lim _ La gq 0, SEC a. (35) 
-_ 20 p o8 


In (34) the interval of integration should include the point representing 


the leading edge; in (33) the integration is taken along a path from p = | 
a point on the upper surface of the wedge. 
We now require a harmonic function which satisfies certain conditions 
on the closed boundary defined by an are of the circle p = 1, an arc of the 
ircle 2p cos@ = k(1-+ p?), and the polar axis 0 = 0,7. 
With 2 pe’ and k (1—k?*)!, the region is transformed into the upper 
half of the z,-plane by the successive transformations (L. 31, 35) 
2k’ | 
Z k ik’ yi — aan Bi (36) 
\'~ —(eFik) 
z, = (x, +1y,) = 4(2°?+Z-*), Z? = z,+(22—1). (37) 


5092.95 


F 











66 W. CHESTER 


The boundaries in the (p, #)-plane defined by the arcs of the curves p = |, 
2pcos@ = k(1+p?) and 6 = 0, 7 now become, respectively, the three 
portions of the real axis defined by —%2 < a, < —l, —1 <2, < 1, and 
1 < x2, < «. The boundary condition (33), which holds for 2, > 1, y, =0 
now becomes (L. 36) 


cp /cp . , Bk Q 
P/#D — 2-a,—1)- .(x,+1)—A(x,—1)]. (38) 
CY4/ CX, l k? 

The condition to be satisfied on the real axis for —o < 2, < —1, which 


now represents the Mach cone, is ¢p/éx, = 0. 
For —1 < x, < 1, y= 0 the boundary condition, in the symmetrical 
and supersonic cases, is 6p/éy, = 0, but in the antisymmetrical subsonic 


problem this is only true for 2) < 2, < 1, where z, is the point corre- 


sponding to the leading edge, i.e. to ¢ —{1—(1—m?)!\/m. Thus 
») . 2 
2(k+-m)? : 
Xo J l = eck ¥. (39) 
(1+km)? 
For —1 < x, < 2, the boundary condition is ép/éx, = 0. 


In addition, relation (34) or (35) must be satisfied in the subsonic case, 
In the z,-plane they become, respectively, 


rot 


. - Op m _ 
lim | Of *dix, = 0, SEC a (40) 
w—0 =, CY", (1 m=)? ; 
; r k'(1+a,)* ep 2) 
and lim | — ae Pode, 5, SEC ax. (41) 
w—0 (1 —2,)* 2'k CY, 


1 
The corresponding condition in the supersonic case arises from the fact 
that p changes discontinuously, by the amount given in (30), at a point on 
the Mach cone represented by ¢ m-+i(1—m-*)!. This point is repre- 
sented in the z,-plane by the same expression as the x, in (39), but now, 
of course, 2, - |. We thus have, by (30), 


; - Op m rs 
lim Me dx, —_—— §8eC a. (42) 
Ao CX, (m*— 1)? 
The boundary conditions have now to be satisfied on the real axis, and 
impose certain conditions on ép/éx, and ép/éy,. 
We consider first the supersonic case. If we put 


Wied a. 44 (43) 
C yy Cc vy 
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then, on y, 0, W must be real for —o < 2, < —1, and wholly imaginary 
for | vy .. 

On the real axis, for x, ‘. 
2+(2,—1)! 
arg W (z,) tan 1} sol 7 eee ——_, 
| Bk(l—k?)-(a,+1)—A(a,—1)J 
1) 
tan 1% -|.teaa-2 = (44) 
Y1 Y2 
V1Yo 2? Bk 23 Bk 
where ance : = A— =e (45) 
YivT Ve 1—k? YT Y2 i— 7 
Using equations (5), (6), (22), and (24), these give 
M2—1 
vs, te (46) 
i = OR+5)i 


In terms of W relation (42) indicates that, near z, = 2p, 


MO SEC a ny 
W ~ A (47) 
ar(m*— 1)#(z,—2p) 
\ function which satisfies all these conditions and in addition is 
integrable at all points in the upper half-plane, except z, = 2p, is (L. 45 


with the sign of 6 changed) 


C8| D(z,—29) 1 











Pw a (48) 
(z?—1)!(z,—29)[ vy, —0(z,— 1) [_ y2o—t(% — 1] 
The constant C is determined by (47) and is given by 
m r=—1\4 
. ae i) Lvs (1—ao) [y+ (1—2%o)*] 
2m(m-+-k)k'(y;+¥3)(¥2+Ys) (49) 
(1 mk)? cos « 
2(m--k)2 ; 
where 3 l—2z, — (90) 


(1-+-mk)?" 
The value of D is obtained from equation (24). Since, on the shock, 


13 


VfXy 
k if 
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we must have 


Bcp dx, 


Sseca = — 
ie CX, 
1 





} (2 oe D(x,;—2Xo)— L(y + Y2)(%,—1)* da, 
J k’\x- 
1 





1} (aj—1)*(%,—2%9)(yi+2,—1)(y$+-2,—1) 
BO8(y, +72) (Dys—1)(yit-ye+y¥s)/"1¥2¥3t+D (59 
—-— —_  — - -- a ~ a2) 
k (vet yaly¥at+vi)(¥i1+72) 


With the help of (49) this gives the following relation, which enables D to 
be calculated: 


(mk+1"" _ Dinitys\yatys)_ vitvetys (53) 
2Bm(m--k) Y172 Y1¥273 rm 
In the subsonic symmetrical case W, satisfies the same boundary conditions, 
but the asymptotic equality corresponding to (47) is now, by virtue of (40), 


1M6,, SEC x i 
= ; (54) 


W,~ - 
a(1—m*)?(z, —Zp) 





near z = 2,. Thus 


VW, = - = 


8 


_ C.84| Dla —%0) — L] tS Ce (55) 
(2{—1)4(%—@9)[y1—1(% — 1) y2—t(2— 1’ 
where C, and D, are still given by (49) and (53). 





The boundary conditions in the subsonic antisymmetrical problem are 
that, on the real axis, w, is real for x, < x and wholly imaginary for 


Xy <x, <1. When x, > 1 relation (44) is to be satisfied. 


in hw 
(21—1)*(z;—2p)*| v1 — (2, — 1) v2—-1(, — 14] 


satisfies all these conditions and is integrable at all points in the upper 





(56) 








half-plane except z; = 2). A singularity in the pressure at z, = 2zj (i.e. at 
the leading edge) seems inevitable: indeed it is expected in a linear theory 
for subsonic flow round a sharp edge. This being so, the boundary condi- 
tions do not define the function completely, and (56) is the solution with 
the weakest singularity at z, = x», subject to being integrable at all other 
points and to being sufficiently general for conditions (24) and (41) to be 
satisfied. Such a singularity in the pressure at a subsonic leading edge (an 
inverse square root) is well known in the linearized theory for steady flow 
(see, for example, reference 5). The fact that W,/5, and W,/5, tend to 
equality as M,cos8 > 1 from below is also encouraging (that this is so 
follows from the fact that 7)» —1 and condition (41) becomes equivalent 
to (40), giving a value of C equal to (49)). Thus, comparing (48) with the 
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um of (55) and (56), and remembering that 6 5, +5,, we see that the 
VW, cos B =]. 

In general, the relations which determine C, and D,, namely (24) and 


yutiol ONntINUOUS at 
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11). give re pectively 











B * (x2+-2)§[ D(x? 2)—1| da -_ 
st (57) 
| J (x+y) (a+ yi)(e +75 
| vy | D (2X, Co) ‘. | dx, 7 
1—a,)'—24k][ y,+(1—2,)' I] y2+ (1—2,)!] 
7 2—2*)'[ D,(a?—ys)+1] dx _ 
y y Pk) (a V1) (2 Y ) 
re 1 sterisk denotes that the finite part of the integral is to be 
6 
7, The pressure distribution on the face of the wedge 
In the plane the upyer face of the wedge is represented by the real axis 
1 (corresponding to —m < 2x k). At a point (x, 0) on the 
wedge we thus have, in the supersonic case and the subsonic symmetrical 
| 4 ; ] : 
/ D(x, Xp) | dx, et (59) 
. Uy) (ey Lo)NY1 + (1 x4) Ye T (1—2,)*| 
1 
nee 4 ero at 2, 1. After some algebra this can be evaluated and 
= es 
yatand 
: coth 1}¥3 = | 
208 i—v3)(¥3—Y3)(2—Ya) \(2—y5)') 
» .. I 
~' SIN®D 
tts ____.. ot 
¥s—Y3)(2—Y3) \(2—y§)*) 
| l y,tand) 2? sind 
'D { tan ae ‘__tan | 
a} i—2) \(y?—2)*) \(yz—2)*} 
| | y,tand 2? sind 
D : . tan J Ya! | tan if = | , 
v¥3—¥3) (Yo— 1) (¥2—2) \(y3—2)') \(y3—2)4J 





(60) 
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“3 1—ke’ 


1—a,\3 c— 2x 

where cos ¢ = ( - =) = “4 (61) 
since, for —1 < x, < l, 

2(k—2x)? 
(1—ka)?" 


The antisymmetrical contribution, on the upper surface of the wedge, is 


%, = 1- (62) 


given by 


ry 




















Opes D,(x;—2% )—1]| dx 
Pa(X) —_ C8, | : | al 1 o)— | . 1 — (63) 
J ( -ay)*(x,—2%p) [y+ (1—2,)* [v2 + (la) ] 
or ~0 
Pa Y1¥2t+¥3 YitVe 
=o ——_ —__—___-_ cot f — —___—__"=__ cosec #4 
2C,8a = v3(Y¥i—3)(¥2—Y3) ¥3(Yi—v3)(Y3—Y3) 
lp + ,; = a - aoe if Yi tan ¢ | -tan if yssiny is 
ee RAI yl GRY lv 
Yi—Y3) (Y2o—Vi)vi-—Y¥3 tf he vi Ne 
ln ye I \ 7 2 tan if Yo tan ¢ \ —tan 1| y3sin pb \ : 
\ *" ¥8—v3) (v2—11) (73-73) \(v3—y8)!) \(y3—y3)}). 
(64) 
» \t 94/()-__ 
where cos & . —_ y 2(k—2) : (65) 
Y3 ¥3(1—ka) 


8. The unyawed wedge 

Equations (13) give the limiting form of the independent variables, and 
all the parameters occurring tend to finite limiting values as 8 > 0. 

The Mach cone now becomes part of a cylinder of radius a, ¢, and its axis, 
in the physical problem, travels downstream with the velocity of the fluid 
behind the shock. The supersonic case now arises when the fluid velocity 
behind the shock is itself supersonic, that is when M, > 1 or, by (2), when 
M > 2-068. The present results for the unyawed supersonic case, and for 
the symmetrical wedge in the subsonic case, agree with those originally 
obtained by Lighthill (2). 

The variation of the quantity 





Po—P,_ _ p{7M*—1)\$ (66) 
5(pi—Po)  d| MP5 | 
on the face of the wedge is graphed in Fig. 3 for 8 = 0 and p,/p,. = 1, 2, 3, 


4-823. These values of the pressure ratio across the incident shock corre- 
spond respectively to M = 1, 1-36277, 1-64751, 2-06809 and to M, = 0, 
0-47245, 0-72739, 1. 

A typical distribution when supersonic flow prevails behind the shock 
is shown in Fig. 6. The scale is adjusted to make the distance between the 
leading edge and the shock the same in all cases. 
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The pressure distributions in the symmetrical and supersonic cases are 
based on calculations kindly made available by Professor Lighthill. 

In the subsonic case the pressure is singular at the leading edge, decreases 
rapidly to a minimum, then rises steadily until the shock is reached. The 





























ia 
1\ | 
1\ 
\ | | 
\ 
\\ = 
ro MA | ieee 
0 € T 
\\\ 
/ 
1:24 ‘a 
/ 
gis 4 
/ 
: A 
0-8} i ecieaiasiiete 
| 
1.4 | 
\ | 
| | 
ii . | 
A 
Leading edge Shock 
Fic. 3. Pressure distribution on the wedge for B 0 and 


various values of the shock strength. 
symmetrical case. 
antisvymmetrical case. 
symmetrical and antisymmetrical contributions are roughly comparable 
for given shock strengths, complete agreement being obtained for the case 
of sonic flow behind the shock. In so far as they differ the symmetrical 
contribution is consistently higher in absolute value except near the edge, 
where the logarithmic singularity is outstripped by the stronger inverse 
square root occurring in the antisymmetrical case. No singularity appears 
in the limiting case p,/pp 1, when the incident shock becomes an acoustic 
pulse. This is a result of the limiting process. If l+e,thenk ~ l—e 
and (62) gives ' 
1l+-x ve 
i t ] ~ de —. (67) 
l—zx 
Further, one can show that, in both the symmetrical and antisymmetrical 
cases, IP } 
CD~ 
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and that C is negligible compared with CD. Consequently the numerators 
of the integrands in (59) and (63) can be replaced by D(x,—2,) (except near 
“, = Xp), and the term whic’ produces the singularity disappears. Since 
Vis Ya. V3 all tend to 2! as « > 0, (59) yields 


9, 2+ #8, ( Ll {4e(1- x))3 dx, ne 
nm J 8.2*| (l—a) } dz 
1 
c P ax 0 7 1 
‘ | —_ l. - 0,{1 Had | (68) 
7 (1+2)2(1 v) a\l—a} 


P, ry (69) 
: 7\1—2x} 


The pressure distributions for p,/p» 1 in Fig. 3 are based on (68) and (69), 

In the case of an acoustic pulse, the pressure distributions are, in fact, 
known for arbitrary values of 5. The extension of (69), for example, when 
6,, is not small would be, in the present notation, 


TD a a7—2 tan | aaa (70) 
= . 2~a 

(see, for example, 7, p. 535). When 6, is small, (69) and (70) agree to the 
first order except at the shock front (x 1), where (69) is singular but 
(70) predicts the well-known doubling of the pressure. Similar conclusions 
are obtained in the symmetrical case (L., p. 469). Thus the limiting distri- 
butions, as described by (68) and (69), are not reliable near the shock. 
This objection, however, should disappear once Mach reflection is well 
established, for then there is neither theoretical nor experimental reason 


to expect unduly large changes at this point. 


9. The effect of yaw 

In Fig. 4 is shown the pressure distributions on the face of the wedge 
for p,/Py = 2,5, = 0, and B = 0°, 40°, and 61°-9, the latter value of 8 being 
the maximum angle of yaw permissible for this particular shock strength. 
The pressure is plotted along a line perpendicular to the axis of the Mach 
cone, and the scale is again adjusted to make the distance between the 
leading edge and the shock the same for all angles of yaw. The computation 
involved in the production of the curves was considerably eased by the fact 


that, for a given shock strength, y,, yo, 


and ys, are all proportional to f, 
C, is proportional to k*, and D, inversely proportional to k?. It is not 
difficult to verify the truth of these assertions from the relevant equations 


(it is also true of C and D in the supersonic case, but not of C, or D,): 
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n (60) the parameter C is absorbed into the coefficients multiplying 


Nn 


Virt 


l 


trigonometrical or hyperbolic functions, then these coefficients 


LW only through the terms (2 y3)t, (v2 2), or (y2 -2)8, 








x + $$ _ 4 _J 
+ “ 
i 4 
~ ant eal yy | 
40 . —T - | 
a SS aaa | 
[a -— -_ = = = | 
A A 
eading edge Shock 
= ition on symmetrical wedge for p,/po 2 
arious angles of yaw 
\odification is necessary in (60) when y, or y, is less than 2+. 
mle en 
| y,tand 
—— tan git ~ 
i—2) (yi—2) 
| y,tand 
id ~ coth-!74 —s 
2—yi) (2—yi) 
I for the other terms. 
reme angle of yaw is approached it is the limiting behaviour 
ist be considered, for the Mach angle tends to Iz. If cosa =e 
juations (14), (28), (50), and (61) we have 
dene ry 
~ tan(8 tL) 
~ 2'‘tan(8—p)-+tan pe 2tV 9 €; 
P 
f ~ 
COs ~™ (tan p x” 2’ )e Le. 


ie of the simple dependence of y,, y2, y3, C, and D on k, it follows 
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that the coefficients in (60) (when multiplied by C) have definite limiting 
values as e > 0. Further 


\(2—y5)') 7. 
and a 
cath = ~ coth-"{1+ 4(—#)e% 
ad 3 


P “ 
~ log(4/e?) —log(73—#*). 

The terms in the last two expressions of (60) behave similarly, and the 

singular terms cancel leaving a finite limiting value for the pressure. 
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Fic. 5. Symmetrical wedge, p,/» 2. Effect of yaw on pressure 


at wedge shock junction. 


It will be observed from Fig. 4 that there is only a slight increase in 
pressure between 8 = 0 and Bf = 40°, but a marked increase at all points 
of the wedge for 8 = 61°-9. This is shown more clearly in Fig. 5, where the 
pressure at the junction of the wedge and shock is plotted against the angle 
of yaw. There is a relatively small increase in pressure up to 8 = 50°, 
beyond which the pressure increases rapidly. 

In Fig. 6 the pressure distribution on the wedge produced by a stronger 
incident shock (p,/p, = 10, M = 2-95200, M, = 1-34463) is shown for the 
two limiting values of 8. When 8 = 0 there is no singularity at the corner, 
for the flow behind the incident shock is now supersonic and is deflected 
discontinuously at the leading edge by a weak shock. In the region of 
uniform flow before the Mach cone is reached the pressure is constant, 
after which it decreases monotonically. However, at the extreme angle 
of yaw (8 = 64°-3) the leading edge lies inside the Mach cone and a 
singularity in the pressure again appears. The transition will occur when 
M, cos 8 = 1, or 8B = 42°, for this particular shock strength. 
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10. Further investigation of the Mach cone 


One of the boundary conditions imposed in deriving the present solution 


is that the pressure should be continuous at the Mach cone. A closer study 


of the pressure field near this boundary shows, however, that the normal 


derivative is singular. For, if ép/én denotes the pressure derivative along 


the inward normal to the circle p = 1 in the (p,@) plane, and if (r,@) are 
polar coordinates in the (v, y) plane, we must have 
cp CPr. ' 7 
p—(p),-. ~ —(1—p) ~ —[2(1—r)}! (71) 
con en 
near } 
Since, as will appear, @p/én is in general finite and non-zero, 
opicr ~ (] r)-* ep/en 
near ? |, which cannot describe the actual behaviour in the fluid. The 


same singularity occurs in Busemann’s cone-field theory (see 5), where 


there is experimental evidence of a shock. In an attempt to reconcile theory 


and experiment in the case of the Busemann cone field, Lighthill (8) has 


evolved an elaborate theory, based on the full non-linear equations of 


motion, by which a more accurate picture of the behaviour near r 


be deduced simply from a knowledge of the linearized solution. 


| may 
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Now the essential difference between the present problem and the Buse. 
mann cone field is that in the latter the whole flow is effect ively irrotationa] 
if terms of the third order of the strength of any shocks which may occur 
are neglected. But it has already been shown that entropy variations in 
the present problem are significant only in the fluid which is swept by the 
non-uniform part of the incident shock, and that the condition of irrota- 
tionality is effectively true of the fluid in the neighbourhood of the Mach 
boundary. Thus the governing equations of flow are the same for both 
problems in the region to be investigated. Since Lighthill’s results depend 
only on the basic equations of motion and the nature of the singularity at 
|, it follows that they may be applied directly in the present problem, 
A detailed description of the procedure will be found in the above-mentioned 
paper and only the final results will be quoted. 

It appears that the actual behaviour of the flow near r | depends on 


the sign of the normal derivative of the pressure there. If 


P2—(P2)r-1 ~ Pr Vi F(A) A—r)! (72) 
asr — 1, then the Mach cone is the linear approximation to a shock wherever 
F(@) > 0. To a second approximation the discontinuity in pressure across 
this shock is "6 

Ap, = fy(y¥+ Pas ee F704), (73 
at(t i1—4}) 


where y here denotes the adiabatic index (8, equations (26) and (39)). The 
strength of the shock is thus O(8?). But when F(@) < 0 the conclusion is 
that, although to a higher approximation there will still be a shock, its 
strength will be at most O(8%) and it will be followed immediately by an 
expansive flow which is relatively stronger (the expansion produces 
pressure changes of order 8?). 

From (71) it follows that 


Cp 1 es 
P2—(Pa)raa ™ UH PrN E [2(] r)*], (74) 
on 
and hence, by (72), that 
F(0) = = roe) (75) 
J i \cn, 


If (75) is substituted in (73) and y is given the value 1-4, one finds that the 
strength of the shock compared with that of the incident shock is 
APs 
Pi-Po 


K?(0) (76) 


if K(?) > 0, where 


‘ > _ ee Ly 
K (6) 3 (2(M N42). (77) 
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find ép/én. This follows at once from (48), (55), 


dy 0. In the supersonic case we have 
Cd} D(x,—2,)—1 | on 
dt  _______. (78) 
-1)4(%9—2,)[ y, + (1—a,)* |] yo +(1—2,)*] 


| contribution to ¢p/én will be a similar expression 


D, 5. The results for the unyawed wedge were 


ntribution is 





CU, 5, D(x, Xo) I "7 
— a ee poms se (79) 
Ce My) lV (1 %3)* yo (1 a,)*| 


been used to investigate the true behaviour at 


limiting cases, namely at zero yaw and at the 


gle of yaw. Of course, in the latter case the Mach 


lach cone degenerates into a plane perpendicular 

flow (in the steady problem). However, Ap, 
1e which provides at least a limiting picture at 
rbed region. 


sider the case of maximum yaw first. Since k 


e have, by virtue of the dependence of y,, ys, y3, 
dz, C.0, D. 
ay a | ee (80) 
CG (ry 1)?(] v,) 
v,)'—2'k} w 24(a?—1)4(1—2,)! (81) 
) ; 9 1 
'~ 0, D,3,( )*. (82) 
l—az, 
1 cone, 
l cos( ‘) 1—cos(@ ‘) 
. “Xx Seas x ; (83) 
I cos(0— x) l cos(@ x) 
hen k + 0, 
1+-sin?6 
t~ —; (84) 
1—sin*°é 
Ce mw — OC, Dd, cos 6 (85) 
ifthe positive square root is taken appropriate to the interval $7 < 6 < 37. 
(85) is substituted in (77) it gives, in the limit as k > 0, 


When 





v 


(C.D \. . 
K(@) = -3i(" "8, 0080, (86) 
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where use has been made of equations (4) and (5). The limiting value of 
C, D,/k is implied, though this quantity is actually independent of k. The 
magnitude of K, for various incident shock pressure ratios may be deduced 
from Table 1. It is singular in the limit p,/p, > 1 (though Ap, itself tends 
to a finite limit). As p,/p, increases, K, decreases rapidly to a minimum 
and thereafter increases steadily. 


TABLE 1 


The value of K, at maximum permissible angle of yaw 








P1/Po l 2 3 10 
K.(@) = a > | > . 
oe 5-963 8-001 12-468 | 26-116 


6,cos@ | 


In the antisymmetrical case, although C, and D, are not strictly propor- 
tional to k® and k-* respectively, as is the case for 


‘and D,, it may be 
deduced from (57) and (58) that the statement is true in an asymptotic 
sense as k > 0. Hence, by (79), 


Pa | dz|C, D,8 


on | df |(1—a,)? 
. (22—1)! 
~ 210, D8, ‘a 
(l—z,) 
~ —C,, D,5,,8in 6 cos 6. (87) 
It follows that, in the limit, 
: . FO: BSS... 
K (@) lim( : ‘) “4 sin 20. (88) 
k->0 k 2 


The range of @ runs from 37 at the incident shock, through z on the 


° 


plane y = 0 to 37 at the shock in the lower half-plane. Thus for a sym- 
metrical wedge with 5,, = 0, the whole of the boundary is a shock, but for 
a semi-infinite flat plate inclined at an angle 6, to the plane y = 0, only 
that part of the boundary in the upper half-plane is a shock if quantities 
as small as 63 are neglected. In the lower half-plane the boundary would 
be a much weaker shock followed immediately by a relatively stronger 
expansion. 

It will be observed that not only is the strength of the shock zero where it 
meets the incident shock, but that Ap,, being proportional to A°(@), increases 
only as the square of the are distance from the triple point. The symmetrical 
contribution is monotonic in $7 < 6 < 7 and attains its maximum at the 
plane y = 0; the antisymmetrical contribution has a maximum at 6 = jr. 
This is by no means the case when 8 = 0, as is shown in Figs. 7 and 8, 
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where A(@) is plotted against the fraction of are distance along the circular 


7 and finishing at the triple point @ = cos-1k 
K (0) is chosen, rather than its square, because the symmetrical and anti- 


symmetrical contributions are then additive). The strength of the reflected 


shock is still zero at the triple point, but even though the figures give the 
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incident shock). 


square root of the pressure discontinuity, it is clear that there is a steep rise 
in strength within 10 per cent. of the are distance from the triple point 
corresponding to an angular distance of roughly 13°). Indeed, when 
);/Po = 2, the maximum strength is attained within 12 per cent. of the 
ire distance from the triple point. It is suggested that this fact reconciles 
theory and experiment to some extent, for it is known that, experimentally, 
‘contact discontinuity appears in the perturbed region starting from the 
triple point, and such discontinuities are associated with three shock con 
igurations. Interferometer photographs also indicate that the strength 
of the contact discontinuity decreases as the triple point is approached 
see, for example, 10, fig. 7), and this would be compatible with a local 


decrease in the strength of the reflected shock. 





80 W. CHESTER 


In the general subsonic case @p/én, and hence the behaviour at the Mach 
cone, will be determined by two contributions which augment each other 
in the upper half-plane (assuming that both 6, and 6, are positive) but tend 
to cancel out in the lower half-plane. One may therefore consider the 


question of where the transition takes place from a shock, to what is 
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effectively an expansion if terms of the third order are ignored. If transition 
is to take place at all ép/én should be negative near the triple point in the 
lower half-plane. As this point is approached x, - — and so, from (78) 
and (79) (with the sign of ép,/en changed), one gets the condition 


} b, F C; D, Cy D,. (89) 


al 
The value of the right-hand side of (89) is given in Table 2 for 8 = 0 and 
various values of the incident shock-pressure ratio p,/p,. It is apparent 
that transition will begin to occur almost immediately after the lower face 


of the wedge becomes positively inclined to the plane y = 0. 
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TABLE 2 


The minimum values of 5,/3, for which a transition point occurs on 
the circular boundary 


DJ 2 3 | 4-823 








C, D/C, D ] 1-068 | 1-049 l 





Table 3 gives the angular distance (in absolute value) from the plane 
j= 0 at which transition occurs when 8 = 0 and 6, = 38,. This may be 
calculated from Fig. 7, or directly from (78) and (79). In the latter case it 
is the root of a cubic. 


TABLE 3 








Anqular distance of pownt of transition from y - (Bp 0,5, — 38,) 
D,/Do 2 3 | 4-823 
6— 180 39°-0 24°-4 16°°-5 | O 





In the supersonic case, when the flow patterns in the upper and lower 
half-planes are independent, transition occurs where the tangent plane 
from the leading edge touches the Mach cone. Not only does Ap, change 
sion but. as shown in Fig. 8, it becomes infinite where the transition takes 
place. This apparent breakdown of the solution can be explained by the 
fact that the tangent plane is a linear approximation to the position of a 
shock or Prandtl—Meyer expansion, and the position would change slightly 

1 second approximation. When the modification is effected it appears 
that the shock (or Prandt|—Meyer expansion), which starts from the leading 
edge, interferes with the Mach cone and prevents 6 from attaining the 
value at which K(@) becomes infinite. The actual behaviour (which is here 
leduced from the analysis of an analogous problem in 8) is as follows. 
With 6 > 0 that part of the Mach cone which lies between the incident 
nd leading-edge shocks is itself a linear approximation to a shock. The 
strength of the latter is in general O(5*) and, except near the transition 
point, deducible from equation (77). The shock increases in strength with @ 
see Fig. 8) and actually attains a strength of the same order as the leading- 
edge shock (namely O(8)) when the two meet. Between the leading-edge 
shock and the wedge there is effectively an expansion of strength O(6*) in 
general, which also increases in strength as the leading-edge shock is 
spproached, but is quickly annihilated once it encounters the latter. 

When 5 < 0 the tangent plane from the leading edge should be replaced 


G 
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Mach cone becomes a shock to a higher approximation (for K(@) is now 
positive). The strength of this shock increases rapidly and becomes ((8) 
when the region of transition is reached. In passing through the Prandtl- 
Meyer expansion the shock is weakened and, when it emerges, becomes 
effectively an expansion whose strength decreases as the incident shock 
is approached. When f = 0 and p,/p, = 10, for example, the strength 
of the shock which replaces part of the Mach cone will be deducible from 
the negative portion of the relevant curve in Fig. 8 (except the part near 
the singularity). 
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(8 THE EFFECT OF RAPID DISTORTION OF A FLUID 


we IN TURBULENT MOTION 
) es 
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noth (Trinity College, Cambridge) 
irom [Received 17 July 1952] 
neal 
SUMMARY 
[his paper is concerned with the calculation of the changes produced in a homo- 
geneous turbulent motion when the fluid is subjected to a superimposed uniform 
distortion. The distortion is assumed to occur so rapidly that the contribution 
) the change in relative position of fluid particles from the turbulence is negligible. 
On this basis G. I. Taylor found the effect of an arbitrary distortion on a purely 
sinusoidal velocity field, and in this paper the effect on homogeneous turbulence 
is obtained by integrating over all such Fourier components. In the case of tur- 
bulence which is initially isotropic, the relative changes in the energies of the three 
7 elocity components are found to be independent of the properties of the turbu- 
lence and are determined numerically. 
Asymptotic results for the important practical case in which one of the three 
951), extension ratios is large compared with unity are presented. For the case of a large 


symmetrical contraction (the common wind-tunnel case), of area ratio c—, the effect 
nany one Fourier component depends on ¢ in a manner which tends asymptotically 
the laws given by Prandtl, but this is not so for the component energies of the 
vhole turbulent motion since the energy after contraction is dominated by a 
lecreasingly small range of wave-numbers. The same asymptotic dependence on c 
s found for all kinds of homogeneous turbulence, although the numerical constants 
vary. When one extension ratio of an otherwise arbitrary distortion becomes large 
mpared with the other two, the asymptotic effect on isotropic turbulence is the 
same, in all important respects, as the effect of a large symmetrical contraction, 
the reason being that the vorticity ultimately is everywhere parallel to the line 
I greatest extension. 
No comparison with experiment is possible, since it can be shown that in none 
f the cases in which measurements have been made is the distortion sufficiently 
rapid for the above linear theory to be valid. Some of the distortions occurring in 
rodynamic practice are rapid enough for the theory to apply, but these will be 
necommon, and a consideration of the non-linear effects is very desirable. 


l. Introduction 

lr has been known empirically for many years that the properties of a 
turbulent motion carried along by a stream can be altered drastically by 
mposing large distortions on the stream. Wind tunnels are usually designed 
nsuch a way that the stream is subjected to a reduction in-cross-sectional 
rea immediately before entering the working section, for the reason that 
such a contraction is known to reduce the ratio of the energy of any turbu- 
lence that may be present to the energy of the stream. A large contraction 
ratio in a wind tunnel is in fact essential if the intensity of the turbulent 


(Quart. Journ. Mech. and Applied Math., Vol. VII, Pt. 1 (1954)] 
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motion is to be reduced below the level at which it affects the characteristics 
of the boundary layers on models at the large Reynolds numbers commonly 


used. Other kinds of stream distortion occur in wind tunnels and ducts and 
pipes, although a contraction in the cross-sectional area of the stream is 
probably the case of greatest practical interest. 

From a different point of view the process of superimposed distortion of 
fluid in turbulent motion may be regarded as occurring in a general turbu- 
lent flow. Large-scale variations of the turbulent velocity produce a 
coordinated strain over a large region of fluid, and so far as small-scale 
velocity variations are concerned this large-scale strain can be regarded as 
a superimposed distortion of the fluid. 

Thus there is some interest in an attempt to determine quantitatively 
the effect of a specified distortion of a region of fluid on the turbulent 
motion of the fluid. The problem is a difficult one, even when very simple 
kinds of distortion and turbulent motion are considered, and it is only on 
a very limited basis that progress can be made. The first work on the 
problem seems to be due to Prandtl (1), who considered the effect of a sym- 
metrical contraction (which reduces the area by a factor c-!) on small steady 
disturbances (which vary with neither time nor distance downstream) 
carried along by a uniform wind-tunnel stream. The lateral velocity com- 
ponents are here determined uniquely by the component of vorticity 
parallel to the stream, and since this vorticity component is increased by 
the factor c' by passage through the contraction, the lateral velocity 
components are increased by the same factor. The energy of these com- 
ponents relative to that of the stream therefore decreases by the factor c—. 
For the longitudinal component Prandtl argued that the pressure drop 
accompanying passage through the contraction is approximately the same 
for all streamlines, so that, to the first order in the disturbance velocity, the 
longitudinal component must be reduced by the factor c-! on passage 
through the contraction. These arguments are obviously valid, but their 
bearing on the effect of the contraction on disturbances which vary with 
distance downstream is uncertain, and no results about the effect on 
turbulence can be inferred with confidence—except perhaps that the 
longitudinal velocity component will be reduced by a contraction to a 
greater extent than the lateral component. 

The first step towards a calculation of the effect of stream distortion on 
turbulence was taken by Taylor (2), who worked out the changes in a small 
disturbance which varies sinusoidally with respect to all three space 
coordinates, for a distortion which is sufficiently rapid. This likewise 
gives an inaccurate guide to the effect of the distortion on turbulence, but 
inasmuch as a sinusoidal disturbance can be regarded as a typical member 
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of the set of Fourier components into which the turbulence can be resolved, 
it is possible to ‘integrate’ Taylor’s result to give the effect of a rapid 
distortion on a genuine turbulent motion. It is the purpose of this paper 
to show how the ‘integration’ may be carried out for the case of isotropic 
turbulence and an arbitrary distortion, and to present the whole problem 
systematically as a preliminary to attempts to treat the case of a distortion 
which is not rapid enough for the following analysis to apply. 

After this paper had been prepared reference (3) came to our notice. 
[t appears that Ribner and Tucker have recently worked out the effect 
of a symmetrical contraction on isotropic turbulence, with methods and 
results similar to those to be described below; these same results are also 
in the text of reference (4) and in the dissertation submitted in 1951 by one 
ofus (I. P) for the degree of Ph.D. of the University of Cambridge. However, 
the present paper goes rather further than any of these references, inasmuch 
as it describes the effect of an arbitrary distortion on isotropic turbulence 
ind the asymptotic effect of a large contraction on homogeneous turbulence, 


ind we feel that it justifies independent publication. 


2. Linearization of the problem for a sufficiently rapid distortion 

The essential meaning of the assumption which will be adopted, following 
Taylor, in order to render the problem linear in the turbulent velocity, is 
that during the process of distortion the displacement history of each fluid 
particle is determined uniquely by the history of the distortion. For this 
wssumption to be valid it is clearly necessary that the influence of the turbu- 
lence on the relative separation of any two fluid particles, whose initial 
separation is given, be negligibly small at all stages of the distortion. It is 
in this sense that we require the distortion to be sufficiently rapid. 

The analysis is best carried out with the aid of Cauchy’s equations for 
the vorticity associated with an element of an incompressible, inviscid fluid 
9), VIZ . 

w (a) mc ‘wi(a) (2.1) 


ca 
J 


the repeated suffix being summed), where w’(a) is the vorticity of the fluid 
element at the beginning of the distortion (time ¢’), when its position is given 


ny 


the vector coordinate a, and w(a) is the vorticity of the same element 
t any stage of the distortion (time ¢), when its position is x. Neglect of 
the effect of viscosity is not in general a justifiable approximation in a 
discussion of the mechanics of turbulent flow in the absence of super- 
imposed stream distortion; nor would it be justifiable in the problem under 
discussion at times long after the distortion commences, for it is well known 
that the effects of inertia and viscous forces on vorticity quickly come into 
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statistical equilibrium. In the problem of distortion of the turbulence, the 
neglect of viscous effects is evidently justified if the process of distortion 
takes place in a time interval small compared with the interval required 
for appreciable viscous decay to occur. That is, if e(t) is the rate of viscous 
dissipation at time ¢, during the distortion, when the energy of the turbu- 
lence} is }u?, we shall require 

t 

| <(t) dt <}u2 ( 

: 


to 
Li] 


for all relevant values of t. An approximate version of (2.2) is obtained by 
using the initial values of the quantities « and u?, and by making use of the 


well-known empirical relation ¢« o (u®)!/1, where / is the length scale of the 
energy-containing eddies of the turbulence, viz. 


: l 
t—t = (2.3) 
Vu’) t=¢ 


If the effect of the distortion on a particular Fourier component of the 
turbulence is to be examined, a condition analogous to (2.2) for that Fourier 
component alone will need to be satisfied. 

teturning to (2.1), contributions to the tensor éx,/¢ ; will arise from the 
velocity field U(a, t) which exists in the iininaiaad ofthe fluid element 
as a result of the superimposed distortion alone (i.e. U is the total velocity 
of the element in the absence of any turbulence) and from the additional 
velocity field u(a,?t) which is a consequence of the turbulence. We have 


t 
aed a= By | (- a +e) dt, (2.4) 
Pa 





oa; oa; ea 


where 6,; is the unit diagonal tensor, and the assumption to be made is that 
the contribution from the second term in the integrand is unimportant for 
t 


- : r OU; »,- , 
the purposes at hand. The quantity | — dt is random, and an appropriate 
oa, 
- 
necessary condition for the assumption to be valid is that 


dt, dt, Oi dt | i dt. (2.5) 


Ca; ca Ca, ca; 
j e 
Pp Pr j Pi 


Since ¢u,/ea; at any time other than ¢’ depends on the effect which the 
distortion has on the turbulence, it is not easy to reduce the condition 


a 
[ [ ee aU, * oU, 


(2.5) to a form involving properties of the turbulence at the initial instant 


T ie. of the velocity fluctuations which are additional to the velocity field produced by 
the distortion in the absence of turbulence. 
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only. A criterion which probably expresses (2.5) correctly until the distor- 
tion alters the order of magnitude of @u;/éa; is obtained by using only the 
' j . OU, Ou; ; ' 
initial value of the integrand, viz. —t—*, where u’ = u(t’). When the 
ca, Ca; 
J J 
turbulence is homogeneous (as will be approximately the case in general, 


so far as velocity derivatives are concerned) the condition thus becomes 


t t 
e(t co ( ou, r OU, 
js | —— df | ——— ey (2.6) 
l ca; ca; 
J J 
i t 
where v is the kinematic viscosity of the fluid. In the simple case of an 


axially symmetrical distortion which multiplies the length of all lines 


parallel to the axis of symmetry by a factor c, the requirement is 


- 2 , 
| | [(c—1)?+ 2(c-4—1)?}, (; 
Vur\lVue/ Jt=t 


which is a more stringent condition than (2.3) in the early stages of the 


bo 
=! 
— 





distortion 

Thus there are two conditions to be satisfied. The first, expressed by 
2.2), requires the duration of the distortion to be so small that viscous 
dissipation does not produce a significant change in the turbulent energy, 
and the second, expressed by (2.5), requires the contribution to the relative 
velocity of two neighbouring particles from the turbulence to be small 
compared with that from the superimposed distortion. When these condi- 
tions are satisfied at all times during the process of distortion, the changes 
in the vorticity of the turbulence are described by the linear equation 


w,(a) = S;,(t)w;(a), (2.8) 
where the distortion tensor S;;(t) is 
t 
* aU, (t) 


S.(t) = 8; - 
ca; 
j 


dt. (2.9) 

t 

The remaining assumption on which the subsequent analysis will be 

based is that @U,(t)/éa, and S;,(t) are independent of a, i.e. that the super- 

imposed distortion is uniform. Then if x and a are referred to the same 
origin, the position of a fluid particle at time ¢ is given by 

Cx; 
x(t) = a;— = a,8,,(t), (2.10) 


' Ca; 
j 


that is, the vector x(¢) is obtained from a by a linear transformation of the 
axes of reference. This assumption is exactly true if the fluid is confined 
within boundaries and the distortion is imposed by a motion of the boun- 
daries such that points on the boundary at any time t are related to points 
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on the boundary at the initial instant ¢’ by a linear transformation, for the 
fluid necessarily then moves in exactly the same way. However, practical 
interest in the effect of distortion on turbulence usually arises in connexion 
with problems in which the boundaries are stationary and the turbulence 
is carried past them, as in a wind-tunnel contraction, and for such problems 
the distortion tensor S;,,(¢) is not exactly uniform. The particles of fluid 
which lie on a lateral plane, far upstream of a contraction, will not lie on 
a plane after being carried downstream (as they should if S;; is to be 
uniform), but will lie on a curved surface with the particles near the centre 
of the tunnel farther downstream than those near the walls. It is difficult 
to make general remarks about the degree of non-uniformity of S;; in such 
cases, but at any rate it is clear that the radius of curvature of the above- 
mentioned. curved surface will be of order of the tunnel width: hence, 
provided the scale of the turbulence is small compared with the tunnel 
width, the distortion tensor will be approximately uniform, although it may 
be necessary to use a different value of S;; at different positions across the 
tunnel. 

The distortion tensor S;; can be written as the sum of a symmetrical 
and an anti-symmetrical second-order tensor, the first representing a pure 
strain and the second representing a rotation. A uniform rotation of the 
whole body of fluid has no effect on the turbulence, except to refer it to 
different axes of reference (since Coriolis forces arising from the presence 
of the turbulent velocities are to be neglected during the distortion), and 
we mey therefore confine attention to symmetrical forms of S;,(t). Later 
we shall choose the principal axes of the pure strain at the time ¢ under 
consideration as axes of reference, so that if e,—1, e,—1, e,—1 are the 
principal extensions (with e,e,e, = | in order to satisfy the incompressi- 
bility condition), the components of the distortion tensor S,,(t) are 

S S50 == € 


1? aa 


é 


n= 2> S33 €3, S,=0 ift?7 J: 


ij 
If T;,(t) is the reciprocal distortion tensor defined by 


TS. = 54, 
the components of 7;; are 


Ty, = le, Ton = ley, Tig = Les, T=0 ifs +3. 
3. Effect of the distortion on the energy spectrum tensor 

The linearization that has been discussed in the preceding section in- 
volves no assumptions about the nature of the turbulent motion, and we 
now turn to this part of the problem. Since the superimposed distortion 
has been assumed to be a uniform pure strain, it is obviously to our advan- 
tage to assume that the turbulence has a similar spatial uniformity; this is 
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in any case a valid approximation for the turbulence in wind-tunnel 
streams, although it will not be valid for fully developed turbulent flow in 
a pipe. We therefore assume the turbulence to be homogeneous. 

Since the problem has been linearized, the standard resolution of the 
turbulent velocity distribution into its Fourier wave-number components 
4) will be useful. We write the turbulent velocity distribution at time ¢’ as 


u’(a) ( e*- dZ'(x), (3.1) 


where \(—1),x is a wave-number vector, and the integration is over all 
wave-number space. The increments dZ’(x) are random (with rather ‘patho- 
logical’ properties), and statistically orthogonal, in view of the homogeneity 
of the turbulence, and the variance of dZ’(x) is proportional to the energy 
spectrum tensor. Equally, the velocity distribution at any time ¢ during 


the distortion can be written as a Fourier integral 


u(a) | e*-9 dZ(x): (3.2) 


since both the superimposed strain and the initial turbulence are spatially 


uniform, it is clear that u(a) will also be a stationary random function of a 
it all values of t, and hence that dZ(x) will have the same general properties 
is dZ'(x 

Now the vector x specifies the position, at time ¢t, of the particle which 
was initially (time t’) at the position a, and hence at any given value of ¢ 
the vector x may be regarded as defining a new coordinate system. For 
coordinate systems obtained by a linear transformation of the Euclidean 
space a, the formulae for vector differential operators remain the same. 
The relation (2.10) between x and a is linear, so that we have, with an 
appropriate change of notation describing functional dependence, 

Cu,(X) OUAX) on 
w,|(X € sik ~ _ €\jk- = Tj; 
CX; ch 
b 
: Gu;(a) 
ind, of course, w,;(a) €iik — . 
t jh ‘ 
ca; 
J 

where e;.,. is the unit alternating tensor. The basic relation (2.8) therefore 


becomes , ) 
u,(X cu ‘ 
U(X) _ pla) & 
Ij mrp um-* 
ca . ca, 


(3.3 


The corresponding relation between the Fourier coefficients of the velocity 
distributions at times ¢ and ?¢’ is obtained by substituting (3.1) and (3.2) 


in (3.3 


€; i.) 1),dZ)(x) € Kk, S dZ,(x). (3.4) 


k mnp nn um 


The coordinate system defined by x provides a more natural framework 
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for Fourier analysis of the velocity distribution at time ft, so that it is 
convenient to define a new wave-number x by the relation 

X.X = x.a, 
that is xi= «KT; and x; = x;S;;, (3.5) 
and to regard dZ as a function of xy. The relation (3.4) then becomes, after 
multiplying by e¢,,, x, (and making use of x; dZ, = 0, which is a consequence 
of the incompressibility condition @u,(x)/éx, = 0), 


te Be op 
dZAX) e iar &mnp Xq*n im dZ,,(x), (3.6) 
x” 





9 


where x? = x’. 
From this relation between the Fourier coefficients ot the turbulent 
velocity distribution before and during the distortion it is possible to deduce 
the effect of the distortion on any of the statistical parameters describing 
the turbulence. Practical interest will usually be concentrated on the 
change in the energy spectrum tensor, and this is the only quantity that 
we shall consider here. The energy spectrum tensor (at the initial instant) 
is defined as (4) po 
;;(%) = lim es 
dx—0 an 


(where dx is the element of volume of wave-number space at position x 


* 
to which the increment dZ(*) corresponds, and dZ is the complex conjugate 
of dZ) and is related to the velocity correlation tensor by 


u,(ajuj(a+r) = [ erm", (xe) dx. 


For the spectrum tensor ®,,(x) at time ¢ we have similar relations, with x 
and x replacing a and x, and the expression for ®,,(x) in terms of ;,(x) 
is found immediately from (3.6) to be 
®; (xX) dy = €5K1 Em np Ejab Euvw X1 Xb X ‘, Ky Siem Pens Dw (%) dx. (3.7) 

This, together with the auxiliary relation (3.5) between x and x, gives 
explicitly the effect of a known distortion on the turbulent energy spectrum. 

The relations (3.6) and (3.7) become more useful when the principal axes 
of the pure strain are used as axes of reference. We then have 

(Ky, Kg, Kg) = (€1 X1> €2 Xa» &3 Xa) 

(so that dx = dy) and (3.6) becomes 





fr 
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1 a Fe 
with two similar equations, where 


dZ,(x) = dZ;(x) — “1 E dZ,(x) +3 dZ,(x) + — aan)|, (3.8) 
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In the case in which 2 ) 1, (3.7) becomes 
l., 2 [/e2 «2\/1 1 2). of. 1\la 
®,,(X) - ®,,(*) + e2 7 (+ ‘\ata-a ss ~<a [Pub) 


2 
1 


2X : 


K5 : ° e To. = 
54 (6 3)(x3 K+ 13) Ml) 
€ 


Wrolnor 


> 9 
KS « . Ca « ‘ , « 
+(e @(t + Sed +K§)D33(%) (3.9) 
€3 X \ 9 


with similar expressions for the two other diagonal components of ®;,(x). 
Since ®,,(y%) is a tensor with Hermitian symmetry which satisfies the con- 
tinuity equation x; D(x 0, it is uniquely determined by its diagonal 
components, and the relation (3.9) is therefore sufficient to determine the 
effect of the distortion on any component of the spectrum tensor. 

[t is at this stage that we can recover the results that have already been 
obtained by Prandtl (1) and Taylor (2). The steady disturbance considered 
by Prandtl does not vary in the direction of the wind-tunnel stream (i.e. in 
the 2,-direction, say) and corresponds to a Fourier component with wave- 
number (0, «s,«,). Putting e, = e, to correspond to an axially symmetrical 


contraction in (3.8) and the two similar equations then gives 


dZ,(x) = ey 1 dZ{(x), dZ.(x) = ef dZ,(x), dZ3(x) = ef dZ5(x), 
with y = (0,e}«s,e}«,), in agreement with the results already quoted. 


The results obtained by Taylor for the effect of the distortion on a purely 
sinusoidal velocity disturbance are identical with (3.8). As Taylor points 








out, if we put e, = eg = ey*, and e, > 0, 
dZ,(x) i ae (3.10) 
: K5-+K3 
‘ K3 1Z,("x)—K AZ, ‘ 2 1Z3(")—Kz d7,/ 
dZ,(y) > eb Be - 1A,\%) dZ3(X) > 4 ks Ky 6 a(X) Kz d () 
° K5-+K5 : 7 K3+K 
(3.11) 


Thus for any one Fourier component the longitudinal and lateral velocity 
components vary as e;* and e} for a sufficiently large symmetrical con- 
traction, which is the same as Prandtl’s result; however, the presence on 
the right-hand sides of (3.10) and (3.11) but not of (3.8), of functions of 
% which have singularities at «3+-«% = 0, shows that this same result may 
not (and, in one case, as will be seen in the next section, does not) hold for 
integrals of the spectrum tensor components over all wave-numbers. 
Although the problem has been solved in principle, the results will be 
useful only if it is possible to integrate the spectrum tensor over all values 
of the wave-number to obtain the effect of the distortion on the energies 
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associated with the different velocity components. It is clear that the 
complexity of the right side of (3.9) will prohibit such integration unless 
;,,(*) has a simple form. For this reason we now assume the turbulence 
to be isotropic initially, and it is again fortunate that this will usually be 
a good approximation in the important case of wind-tunnel turbulence 
immediately upstream of a contraction. When the turbulence is initially 
isotropic we have 
(x78 ;;—K; K;), (3.12) 


where E(x) is (initially) the density of kinetic energy with respect to wave- 


> 


x 
number magnitude (so that | E(x) d« = }u*). In this case (3.9) becomes, 


0 
2 Ke Ke K2 2 
2[ Ke Ss, "4 Q 15 
Ke{— + 7)+4-f—2-t- or |, (3.13) 
ef es : eS e; 


and we find for the density in wave-number space ofcontributions to uw? 


after some reduction, 


®1,(X) 





®5(x) BM) ud(e 


4rk?y? 


nore 
~ 
core 
— 
x 
rope 
wr 


+-e#) + «3(e? +3). (3.14) 


The turbulence is rendered anisotropic by the distortion, so that it is no 
longer possible to express the spectrum tensor in the form (3.12). 

The effect of the distortion on the energy associated with any velocity 
component can now be found by integration with respect to x, the results 
of which are described in the following sections. 


4. The case of a symmetrical contraction 

The effect of a general distortion on the energy of turbulence which is 
initially isotropic can be found, but we shall first consider the simple and 
important case in which 


. se" 


c being the ‘contraction ratio’, or the factor by which the length of a 
material line element in the direction of the axis of symmetry is multiplied. 
Equation (3.13), and the sum of its two partners, now reduce to 








K2, l , .«e+K2 
®1,(X) = xi Py1(%) = 5 E(x) ma . (4.1) 
, cx!- c 342 kK? 9 
Dy0(X) + Pg3(X) = E(x)— =~ (4.2) 
Aor K*y4 
where x? = c~®«7+¢(«3-+ x2) 
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The ratio of the energy of the longitudinal velocity component (i.e. along 
the axis of symmetry) at time ¢ to that before the distortion is 








a 9 2 
1, \KetKs 
Z U(«c) 2——3 dx 
uz J D14(x) dx x* 
by a [oO ] oe oo 
Uy J 11(%) an | E(x) 2 3 dx 
K 
With the use of polar coordinates x, 0, where x, = « cos 6, («3+-«2)! = «sin 8, 


we find that the integrations with respect to « cancel, giving 





T 7 
° “ina r° 
sin?@ dé al 
Ly | | sin°0 d@ 
J (c-*cos*é-+-csin*@)*/ | 
0 0 
1+? 1+-a i. ‘ 
3 | en LOE (4.3) 
203 l—a 
where «? i—c-’. Likewise the ratio of the energy of the two lateral 


velocity components at time ¢ to that before distortion is 











po ext ten3nd x? 
FE (1c) -2 ant Aae 
us+us | [®,.(x)-4 D55(x)] dx | (x) mo ( 
EWE [[D alo) + ybe)] de TG a 
K 
: 0870 F 
( a one ——. |sin 6 dé 
! c8(sin*6+-c-3 cos?) 
; 4 
| (1+-cos?@)sin 6 dé 
0 
l x? 1 +- a) 
- 3 —2i 1 2 y 
aC 1 4° (40 : a lee = (4.4) 


The ratio of the total energy of the turbulence at time ¢ to that before 
distortion is, from (4.3) and (4.4), 





(4.5) 


a (44 +pet+pMs) = $¢+— log 
UU. 40a — Ot 

Numerical values of », and uw, when c > 1, as given by (4.3) and (4.4), 
are exhibited in Fig. 1. , and », are found to be always less than and greater 
than 1 respectively, as might be expected from a rough use of Prandtl’s 
arguments (although Taylor was led to conclude, from the results for a 
single Fourier component, that », would first decrease and later increase 
with increase of c). The increase in the lateral energy dominates the total 
energy of the turbulence, which is always increased by the contraction. 
However, the energy of the turbulence relative to that of the stream which 
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transports it through the contraction in a wind tunnel—and it is the 
relative energy which is significant so far as the effect of turbulence as an 
aerodynamic disturbance is concerned—is altered by the factors 2,/c? and 
H2/c*, both of which are always less than 1 (provided c > 1). In many 
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Contraction ratio C 


Fic. 1. Effect of a symmetrical contraction on the energies of the longi- 
tudinal and lateral components of velocity. 


practical cases c will be large compared with 1, and expansions of py, and 
H in powers of c~* are useful. We find 


Py = $c-*[(1+-c-3+ fc-6+ ... log 4e3— (14 2c-34 7c ...)], 
be = $e[ (1+ $c-3+- 4c +...) (fe... log 4c], 


and for c greater than about 2, the approximation 
Hy = #c-*(log 4c?— 1), jt, = Be (4.6) 
will usually be quite good enough. 

It appears that the distortion factor x, does not vary with c, when c is 
large, in the manner suggested by Prandtl and Taylor (viz. x, 0 c-?). The 
reason, as already hinted at in the preceding section, is that when c is large 
the expression for { ®,,(y%) dy is dominated by a region of wave-number 
space in the neighbourhood of x3+-,3 = 0 (Fourier components with wave- 
number vectors in this region being much less suppressed by the distortion), 
and the extent of this region itself varies with c. Thus the asymptotic 
dependence on c of the effect of the distortion on the energy of any one 
Fourier component is not the same as the behaviour of the distortion factor 
for an integral over all wave-number space. A similar cause gives rise to 
a logarithmic term in the expression for 4, but this is not the leading term 
when c is large. 


The asymptotic dependence on c of the factors x, and yz, can be shown 
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to have the same form for any kind of homogeneous turbulence, whether 
isotropic or not. Werefer back to (3.9) and find that whene, = e, = ey * = c~?, 


A 
®1,(X) = a Pu) (4.7) 


with no assumption about the turbulence. Likewise the two other equations 


of which (3.9) is typical, give 
Ki(K3+-K3) r 
| , , ail 2 " 
Do9(%) + D3(X) = cl Poo(%)+ D5(x) | éF— },(). (4.8) 
xX 
The distortion factors for the longitudinal and lateral energies are now 
obtained by integration over all values of x (or x). Taking spherical polar 


coordinates x, 0,4, where 6 = 0 is the x,-axis, we have for the longitudinal 








energy, eae 
me ; sin 0 
| | | KD, (x)- —— 
a3 3 (c-* cos?@-+-¢ sin76)* 
0 0 0 
By 7 20 
i | x?) ,(%)sin ¢ d¢dédx 
00 0 
P sin 
¥'(@) ey 
P (c~* cos*6-+-c¢ sin?6)* 
0 ’ (4.9) 
| ¥'(@)sin36 dé 
where 
4'(0) — D,,(*)x* dddk, Y'(6)sin3é dé = “,". 
sin*¢ 
0 0 0 


Since },(%) is necessarily an even function ofx (by definition), and vanishes 
with («3+-«2)! (from continuity), ‘'(@) will be finite at 6 = 0 (for isotropic 
turbulence (6) $u,", the dependence on @ being trivial in this case). 
Hence, when c is very large, and the value of the integral in the numerator 
of (4.9) is dominated by small ranges of 6 near 6 = 0 and 6 = 7, we shall 





have = 
(0) sin@ 
‘at ios Bp} ° a) 40 
i? | (c-* cos*é+c¢sin*6)* 
0 
(0) ; 
— -_“¢-2 log c3, (4.10) 


Similarly the expression for the distortion factor for the lateral energy is, 


from (4.8) 


T © 
c a c2 sin®@ cos?@ 
Me C — 4 (0) —; 9D . 9 _¢ 
: ite (c-* cos?6+-c¢ sin?6@)? 


16, (4.11) 
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which is well behaved when c is large; as c > «, 
Hb, > c| 1— ———. | '¥'(6)sin 6 cos?6 dé}. (4.12) 


Thus the asymptotic dependence of y, and p, on c has the same form for 
all kinds of turbulence, although the constants of proportionality vary, 
Somewhat surprisingly, the constants depend only on the energy spectrum 
of the component of velocity parallel to the axis of symmetry. 

There do not appear to be any records of systematic measurements with 
which the results of this section can usefully be compared, although a 
number of incidental measurements demonstrating the damping effect of 
a contraction in a wind tunnel are known. D. C. MacPhail has made a set 
of measurements} of the mean squares of the velocity components w, and u, 
at a number of positions on the axis of a symmetrical contraction in a small 
wind tunnel, but this is not a case in which the conditions (2.3) and (2.7) for 
validity of the linear theory are satisfied and quantitative agreement with 
the theory is not good. Approximate values of / and (u?)! before contraction 
in MacPhail’s experiments were 0-04 ft. and 0-25 ft. sec.—! respectively, so 
that the condition (2.7) requires 


t—t' < 0-02[(c—1)?+ 2(c-!—1)?}! see. 
The linear theory will therefore be applicable, in this case, to positions 
whose axial distance d from the beginning of the contraction is such that 


d 0:02U 


l l 





[(c—1)?-+ 2(e-4#—1)2]}! = 1-7[(c—1)?+- 2(e#—1)?}}, 


where U ( = 3-5 ft. sec.~1) is the wind-speed upstream of the contraction. 
The contraction ratio c is thus required to become large compared with 
unity in a distance of a few multiples of the scale of the turbulence, and 
this was not the case with the contraction used. The value of (u?)!/U is 
here rather large, which makes the condition (2.7) difficult to satisfy, but 
even so a contraction in a wind tunnel will have to be unusually rapid 
before the linear theory applies accurately. 


5. Effect of an arbitrary distortion on the energy of the turbulence 

As was remarked at the beginning of the preceding section, if a field of 
isotropic turbulence is subjected to a rapid distortion, the energy associated 
with each component of velocity in the resulting turbulence may be calcu- 
lated for the general case in which no relation (other than the incom- 
pressibility condition e,e,¢, = 1) exists between the extension ratios. The 


+ Carried out in Cambridge in 1943 and described in reference (6). 
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typical quantity to be calculated, as a function of any two of the variables 


a | 


R a | D14(X) dy 
2 | Dil) dx 


ere D,,(%) has the isotropic form (3.12) and ®,,(x) is given by (3.13). As 


result of the initial isotropy of the turbulence and the homogeneity of the 


listortion, the integrations with respect to the wave-number magnitude x 
incel. and is given by the integral 
} "(( {«e - ce cel\*) face cs. ee * .. = 
BiG J+(3 q) | at—at—3) @S(x) (5.1) 
Smet JIU le a Nea NG a 


er the surface of a sphere of radius « in *-space. 


\fter evaluating the line integral round the circle 
K const., Ky const. 


yy changing the variable of integration to x 3/K., (5.1) becomes, after some 

loebra, 
3] 2 l ( A, B4-Agl 

(5, Bs)’ 


> 
J. ~ 
2 dk, , (6.2 





where 


(5.3) 


Since (5.3) shows the integrand in (5.2) to be a rational function of «, and 
B, B (quartic in «,)’, 


t follows that 4, may be expressed in terms of elementary functions and 


] 


standard elliptic integrals. Thus, with the formal transformation 


Ky 


t XY—, (5.4) 


bore 


nd the notation 2 | 


(5.5) 


aw 


ind leaving aside for the moment the question of the relative magnitudes 


fe, €9, and eg, and the consequent reality of « and k, (5.2) becomes 


5092.25 
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where X and Y are the incomplete elliptic integrals 


x 
a 


dt . 1— k*t?\3 
Zs snitaeaiilondact a iam aieiaas: r a. 5.7) 
(1 —#?)#(1— #2)! ' [ ( l =] ; ™ 


« 


0 0 
Moreover, since the formulae for x, and yz, follow from a simple interchange 
of suffixes, and since numerical tables of X and Y are available. (5.6) 
represents the complete formal solution to the problem. 

However, although the expressions for 1, 2, and x3; must obviously be 
real whatever the relative magnitudes of the extension ratios, it is impor- 
tant for practical purposes that the individual factors in these expressions 
should be real. This is particularly true of the elliptic integrals which are 
only tabulated under the conditions 


Oa <= 7, O<— = I. (5.8) 


For any prescribed set of extension ratios, only one of the expressions for 
the p’s, as given by (5.6), satisfies the condition (5.8) (or its equivalent after 
an interchange of suffixes). Hence alternative formulae must be found for 
two of the ps. 


Now, to fix ideas, we may choose 
e, 5 C3 (5.9) 


without loss of generality. Thus, « and / are real quantities satisfying 
(5.8), and for analytical purposes it is convenient to regard these two 
quantities as the two parameters describing the distortion (in a sym- 
metrical contraction k | and a reduces to its meaning in (4.3)); (5.6) then 
represents the required form of the solution for the energy associated with 
the velocity component in the direction of greatest extension, and it 
remains to express pz, and ps in terms of elliptic integrals with argument « 
and modulus k. 

A cyclic interchange of suffixes in (5.4), (5.5), and (5.6) shows that p, 
involves elliptic integrals X’ and Y’, corresponding to X and Y, with 
argument 


9 9 
a es—es 
and modulus k’2 ain. J 
€3 ey 
‘ k2o2 
Thus x 2 = — 
1— k?a? 
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when arguments and moduli are related in this way, it may be shown that 











x’—} ale k?)X —Y}+ ie2(e? —e3)! 
where ? V ( 
(5.10) 
2-0; 
ps 
f\ 
Fic. 2. Intensity changes in an arbitrary distortion. 
\ similar procedure yields the result 
3/1 e e, *—e, *—e,* ,, . 
bs ar A he = — J ° (5.11) 
bey es—ez — (€2 — €3)(€y — €3) 
In Fig. 2 the information given by the solutions (5.6), (5.10), and (5.11) is 


shown graphically. In order to present all the information on one figure, 


is considered as a function of e, and ¢, without restriction on the relative 


magnitudes of the extension ratios. In this way the above-mentioned three 
forms of the solution always refer to p, for different ranges of values of e, 
ind ¢,. The figure shows the variation of , with e, for constant values 
ot ¢ 


Apart from being in a form suitable for computation, the results (5.6), 
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(5.10), and (5.11) provide a convenient expression for the calculation of 


the change in the total energy of the turbulence due to distortion. Thus, 





5 
3 


= 5 (Py Ha Ms) 


9 


L{ J e2 . 
3 , 2 2)] ma 
oie | (e2— @2) X + (eF €3) ) . (0.12) 
-1¢ 1&3 
This equation, together with (5.6), may be used to recover the results of the 
preceding section. For a symmetrical contraction k? 1, so that 
. l+o : 
X log. . ) x 
‘a l x 


and the relevant expressions for the u’s follow immediately. 

We may also note the expressions for ;.,, 9, fg for a distortion in which 
one of the extension ratios remain unity. This case, which corresponds to 
the flow of a turbulent stream through a passage whose cross-sectional area 
remains constant, has been studied experimentally by Townsend and the 
results are described in the accompanying paper in this journal (7). Putting 

e. b, é 


— ee 
1 €3 c l, 


the results (5.6), (5.10), (5.11), and (5.12) become 





3 c(ct +e? 1) 
es a a (eee 13) 
By val (c2—1)(c4 ni! p. — 
5 c(c*§—c*4-1) 
3 ae me ps —— {(¢?+ 1)y— x}— 1], 5.14 
Me = Fe%Xc2 | (ip WOT | sah 
3 c(c*—c*— 1) 
me ees) pe re nae i 4") 5.15) 
bs Ve 5| - (41) y| (5.15 
and pe = te 2 Cc 1(¢4 1) ba 7-C 1(¢4 1)ty|, (5.16) 
respectively, where 
(l—c 4) 
a= | {l—c#?/(1+c?)}-(1—#)-! dt, (5.17) 
6 
(1—c—4)4 
y = } {l—c?t?/(1-+-c?)}4(1—#?)-+ dt. (5.18) 
0 


Numerical values of 4, 42, and ps, as given by the above formulae, are 
shown in Fig. 3. 
When c is large compared with 1, 
pw, = 3c-l(log 4c—1)+ O(c log c) 


bl, = 3c—-3c-1(log 4c—$)-+- O(c loge) }, (5.19) 


bg = 3c+-3c-1(log 4c—})+ O(c 3 logc) 
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nd when ¢ greater than about 2 the terms written out explicitly will 
usually provide a sufficiently good approximation to the solution. It may 
be noted that, when c is sufficiently large, the values of py and ps, become 
pproximately equal to each other and to the value appropriate to a 
symmetrical contraction about the x-axis (cf. (4.6)).. The explanation is 
when e, is much larger than e, and es, the vortex lines lie nearly 

00 430 
25 
ps 
ul and 
Ps 
2:0 
1:5 
{ 
1:0 
2:0 25 3-0 
C 
) ( nt ty changes in a constant-area distortion. 
} 
parallel to the v,-axis and the dominant contribution to the change in the 
ponents in the 2,- and 2,-directions arises from further stretch- 
Ing ot these vortex lines 
Sol pproximate forms of the general solution are also of interest. 
‘hus, for small distortions, as e,, é,, and es all tend to unity subject to 
the condition (5.9), a tends to zero and k? tends to a real limit between 0 
| d 1. Hence, expressing (5.12) in terms of « and k only, we obtain 
| oom ’ l ao , m 
. | — k?a2)!(1— 1?) f(l—a®)(X—Y)-+Y} 
X 
whicl fter a little algebra, may be expanded in the form 
1+ 2(1—k?+ k4)at+ O(a8), 
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for small values of «. In terms of the extension ratios, the above expansion 
takes the more symmetrical form 


— <9? a — . 
w= 1+dL(e{—e8)?+ (es—e§)? + (e§—e9)?] +... - (5.20) 
The change in total energy is thus of the second order in e, —1, ¢,—1, e,—1, 


and to this order is always an increase. 

On the other hand, expansions for large distortions will usually be of 
more practical value. In this case e, is large and e¢, is small, so that « tends 
to unity. The modulus k may tend to any limit, though this will also be 
unity unless e, is of the order of e,. Considering the latter case first, the 
elliptic integrals X and Y tend to their complete forms 


1 1 


X\(k) = [ U—Re)—eyt dt, ¥4(k) = [ ey) dt, 


and no particular simplification occurs in the expressions for the ,’s. 
However, it may be noted from (5.6), (5.10), and (5.11) that 


Hg > HB +P > $61 Xo(h), (5.22) 


where, from (5.21), l<YV,(k) < 4a 
so that the increase in total energy is very largely determined by e,. As 
€) > €, (i.e. as the distortion becomes an expansion approximately sym- 


metrical about the x-axis), (5.22) becomes 


, 3a 5 98 
as ee 16 (5.23) 
When e, > ey, k tends to unity and, on putting a? = 1—e?, h? = 1—k”, 


where « and k’ are small, we find, after some algebra, that 


{> 


Bt > log4/(« (k’?4 €“) ); - log{4e, (€.+€s)} 


Thus, the relevant expressions for the p’s are 
bey > $U(E3 +63) /ey flog (4e,/(¢.+ es) IS | a) 
, (5.24) 

Me -> fey > 4ey 

As a check on the algebra, we may note that (5.24) reduces to (5.19) when 
by 1, and to (4.6) when e, é;. The results (5.22) and (5.24) show that 
the total energy always increases linearly with e, when e, is large and that 
the constant of proportionality can only vary between quite narrow limits. 
Moreover, half the total energy always resides in the 2,-component of 
velocity, and unless ¢, is of the same order of largeness as ¢, the other half 


resides in the x,-component of velocity. Thus the distribution of energy 


usually tends to become symmetrical about the x-axis, and the symmetrical 
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contraction is more typical of the general case than might at first be 
supposed. 

Analogous asymptotic results for the effect of large distortions on general 
homogeneous turbulence could probably be deduced in the way that 
1.10) and (4.12) were obtained for a large symmetrical contraction. 

While this paper was being prepared, one of us (I. P.) was the holder of 
, Senior Research Award from the Department of Scientific and Industrial 


Research 
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THE UNIFORM DISTORTION OF HOMOGENEOUS 
TURBULENCE 
By A. A. TOWNSEND (Emmanuel College, Cambridge) 
[Received 29 July 1952] 


SUMMARY 

By passing the turbulence found downstream of a uniform grid through a constant- 
area distorting duct, the effect on isotropic turbulence of a suddenly imposed uniform 
mean rate of strain has been studied. A complete theory exists for the limiting case 
of an infinitely rapid strain, and the relevance of this theory to the experimental 
conditions is discussed. It is concluded that the rate of strain is too small to justify 
application of the theory, but that, for small strains, the structure of the energy- 
containing eddies is similar to that produced by an instantaneous strain. For larger 
total strains the turbulence develops an equilibrium structure in which the intensities 
are still changing, and the similarity of this structure to the structure of ordinary 
shear flows is considered. Finally it is found that, during the distortion, the dissi- 
pating eddies are not in a state of local isotropy, although the return to isotropy is 
rapid when the distortion ceases. This is interpreted as meaning that the direct 
influence of the mean rate of strain is appreciable at these Reynolds numbers, but 
that even highly anisotropic turbulence dissipates its energy through nearly isotropic 
eddies. 


Introduction 

THE persistence of the hypothesis that eddy transport of momentum may 
be described conveniently by assigning to the flow a coefficient of eddy 
viscosity is to be attributed to the general similarity of the irregular 
molecular motion which is responsible for the phenomenon of viscosity 
in gases and the irregular velocity pattern characteristic of a turbulent 
flow. In general, such a coefficient, defined as the quotient of the turbulent 
shear stress divided by the rate of shear, varies with position in any par- 
ticular flow and depends on the Reynolds number of the flow, and clearly 
lacks the uniqueness of the coefficient of molecular viscosity. The principal 
reason for this behaviour is that the local turbulent velocities depend on 
transfer of energy from the mean flow acting through the turbulent stresses. 
and naturally this energy transfer varies greatly. The mixing length 
theories attempt to predict the variation by assuming that the turbulent 
intensity, and so the eddy viscosity, are determined by local conditions 
only, but recent work (Townsend (1)) has shown that the turbulent intensity 
at any point is strongly influenced by diffusion of energy from other parts 
of the flow, and this particular hypothesis seems to be without foundation. 
The relevance‘ of the eddy viscosity hypothesis as distinct from methods 


for its determination has been stressed by the observation of Reichardt (2) 


[Quart. Journ. Mech. and Applied Math., Vol. VII, Pt. 1 (1954)] 
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that the mean flow profiles of free turbulent flows, i.e. wakes, jets, and free 


mixing zones, may be described by assigning to any section a constant 


oefticient of eddy viscosity, and the author (3) has shown that, for a 


vlinder wake, the remaining discrepancies are completely removed by 


taking account of the observed intermittency in the outer parts of the wake. 


[he known presence of large-scale mixing motions in these flows suggests 


that the fully turbulent fluid is almost homogeneous in intensity and 


capable of developing a shear stress proportional to its instantaneous rate 


f shear, the constant of proportionality, i.e. the eddy viscosity, then being 
, physically significant property of the turbulent fluid. Such rapidity of 
djustment of turbulent shear stresses to changing conditions is difficult 
to reconcile with the growing number of observations of turbulent motion 
which tend to show turbulence as a relatively slow-moving and stable 
motion, at least in so far as the energy-containing eddies are concerned. 

It is difficult to resolve these doubts by measurements in ordinary shear 
flows. and so an investigation of an extremely simple shear flow was under- 
taken. This flow is the uniform distortion of initially isotropic turbulence. 
carried out by changing the shape but not the cross-sectional area of the 
wind-tunnel duct, and is very nearly a homogeneous flow without the 
complicating factors present in ordinary flows. Somewhat similar experi- 
ments have been reported by MacPhail (4), but his results are difficult 
to interpret, and a repetition using more modern equipment and better 


experimental conditions appeared necessary. 


Notation 


The notation used is conventional for the description of homogeneous 
turbulent flows, using a Cartesian coordinate system and the following 


S\ mbols 


U,VLi components along the axes Ox, Oy, Oz of the mean 
velocity at a point, 
components of the instantaneous turbulent velocity 


fluctuation at a point. 


é components of the instantaneous vorticity fluctuation, 
‘‘ the mean pressure, 
the instantaneous pressure fluctuation, 


the fluid density, 
the kinematic viscosity, 
v. B strain ratios along the principal axes of strain, 


components of the vector wave-number k, 
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F(k), G(k), H(k) the turbulent intensity spectrum functions for u, v, w 
respectively, 


ys, (1) the one-dimensional spectrum function for wu with 
respect to x, 
€ the rate of conversion of turbulent energy into heat, 
L5vu2 
2 
d2 = 
€ 
Cu\? 
d2 = ur (: -l, 
\C x 
L,, the longitudinal integral scale for u, defined by 
ox 


~ oe 
wl, = — | | F(0,m.n) dmdn, 
7 
a 


bs ds, the corresponding scales for v and w. 


For the experimental flow the origin is taken to be in the plane of the 
turbulence-producing grid and on the axis of the duct. Oz lies along this 
axis, and Oy is parallel to the direction of compression of the duct section. 


Effect of an instantaneous distortion 

[t is interesting to inspect the consequences of assuming that the dis- 
tortion of a homogeneous turbulent field is done so rapidly that movement 
of the fluid particles during the distortion is entirely due to the distortion 
and not at all to the turbulent motion. When this is true, viscous and 
inertial transfer terms in the equations of motion are negligible compared 
with the pressure terms, and the problem can be solved. If the distortion 
is a finite strain of principal strain ratios x, 8, y along the axes of reference 
(aPy | if the fluid is incompressible), a particle originally at the point 
Lo. Yo: %p MOVES tO aay, BY, yZp, and the vorticity at the particle changes 
from & , Ho; Cy to ao, By. yoy. Given the original vorticity distribution, the 
final vorticity distribution can be written down, and, from this, the final 
velocity distribution. Taylor (5) did this for a single periodic motion 


in three dimensions, and applied his results to the change in intensity of 


turbulence passing through the contraction of a wind-tunnel. Batchelor 
and Proudman (6) have extended this work to give a complete solution 
for the instantaneous distortion of initially isotropic turbulence. For a 
uniform plane shear, defined by 


x . B B, y *, 
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the spectral density functions for uw, v, w are given by 








es ; . . m2(12 + m2 n2(12+-n2)] ) 
i F k4 (ky) | m2+ n2+(B-*—-1)—® of mo) (B! ye to) | 
; 0 | 

1403 4 b(ke) | n2-+-12-+(B2—1) ng(mi +75) 4 (p21 alma 2 - 
ke ) ke | 
k4H hk d(k,)| l2-+- m2 B2 1) Lolns +1). (g 2 1) mano mo) 
‘ *0O “— 0 0 f k2 | ke J 

where k2 ke + (B-? 1 )m?+ (B? 1)n?. 


\s a consequence of isotropy, the initial spectrum functions are 
Fy = (m§+-n§)b(ko) 
Gy = (no +Mg)b (Ko), 
Hy = (l§+-m§)b(ko). 


t is clear from these equations that the turbulent intensities, defined by 
u* ||| F(l.m,n) dldmdn, etc., 


lepend only on the original intensity and the strain ratio, 8. For this 
rticular type of strain the integration is not simple, but the results for 





smal] total strains are likely to be the most interesting for considering a 
iform rate of strain. Provided B?+ 8-?—2 is small, the intensities are 
u - 8 (B2-+-B-2—2)], (2) 
u vy] (6?+-B-*—2)| (3) 
- B2—1, 5 

u 0+ U9) as 1% 73,(B? B-2 2)}. (4) 

J 
this flow the directions of principal strain are the coordinate axes, unlike 
inary shear flows for which they are at 45° to the conventional axes, 
nd the shear correlation coefficient of ordinary shear flows is in the present 

ced by 
y2— ye? 4 B?>—] , 
/ — oe (2) 
vet we ip -+- i 

the same degree of upproximation. 
‘he integral scales for the three components with respect to the Ox 
ion are determined by the intensities at zero wave-number of the 
ppropriate one-dimensional spectrum functions, which turn out to be 
inchanged by the distortion. For any uniform distortion with no strain in 


Ox direction. these scales are given by 


ur D (6) 


Uo” 
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Since the intensities undergo considerable changes, the one-dimensional] 
spectra must change in shape, and it is found that, for not too small strains 


ys,,(1) develops a maximum at a non-zero wave-number. 


Effects of a finite rate of distortion 

The experimental flow differs from that considered above in several 
respects. The most important difference is that the distortion of any 
particular fluid volume is carried out at a finite rate, and processes neglected 
in the instantaneous theory are now of considerable importance. Further, 
the flow is not completely homogeneous since intensity gradients exist in 
the direction of flow, and the initiation of even a finite rate of strain cannot 
be done instantaneously. Some information may be gained from the 
equations of motion. 

The mean velocity distribution not too near the entrance or exit of a 
distorting duct whose linear dimensions in the Ox and Oy directions are 
given by 4 


A, ¢ ax 


> > . 
B B, e®* 


” U constant 
V aUy >}. 7) 
W al’z 


The equations of mean flow are, in full, 


Cuz cuv  cuw leP 
OX cy C2 pc H 4 
anve cul ov* cuw leP 
azl “yy a —s vi : (S 
Ca Cy C2 p Cy | 
a , | 
ae cuw  evw . cu LeP | 
a? lz | apt = 
Ox cy oz p Cz 


Along the Oxv-axis, symmetry reduces these to 


cuz cuv . cuw i oP 
Ox cy C2 p cx 
LeP | 
0 a (9) 
| 
i 
p cy 
| 
1oP 
0] — 


Pp Ce 
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Simila the equations for the turbulent intensities along the axis are 
eu? l cu lL curv lL oum*w 1 cop —— 
] | af uU fo vuv-u 
! 2 2 7] 2 C2 p CX 
l ov? lL ov=w 1 cop = 
/ U7 4:2 pt voV2v - (10) 
2 CY 2 ©z p cy | 
Cu | Cow? lL ou? | Cp 
{ Al we | ao Ww i +L pwV2w 
? 4 Z oy 2 oc p Cz 
It is possible to neglect the diffusion terms involving gradients in the Ox 
rection since the turbulence level is low. The reasoning is that Cuv?/éx is 
der which is small compared with U(éu?/ex). Similar reason- 


¢ shows that, provided the intensities remain reasonably constant over 


cross-section of the duct, which is actually the case, diffusion terms 


volving gradients in the Oy and Oz directions are also negligible. The 
juations now become 
Cu? 1 op 3 
——U — vuv2u 
: CX p ox 
; v " 1 op " 
1] alr v -peV2r ; (11) 
Pp C2 p cy 
7 ; cu q l cop > 
L | aUue w+ vwV2w 
| a Ca p OZ 
nd the equation for the total intensity, gq? = u?+-v?+-u®, is 
cf ~ 5 wa es 7... 
1p 4 aU (v2@—w? pf uV2utvV20+wVw], (12) 
CA 


since the sum of the pressure terms may be shown to be small by reasoning 
similar to that used above in connexion with the diffusion terms. 

These equations confirm that the experimental flow is nearly homo- 
geneous and that it reproduces the theoretical box-type distortion at a 


finite uniform rate well enough for ordinary purposes. 


Finally, the equation for v?—w?, which is a measure of the turbulent shear 
STress 1S 
} , ( l c ) ei Cp Se 
1] w2)—al (v2-+ w?) _ (v I —wt) Lp(vV2v—wV2w). (13) 
’ C2 p cy C2 


To compare this with the results obtained by assuming an instantaneous 


distortion, the tota! strain ratio B is equated to e-“*-*) where x xy is the 





section at which the distortion begins, and then equation (4) becomes, for 


) small values of a(a—zx 


1) 





a(x 
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“ (v2— w?) 2a(v2+-w®). (14) 


l 
or = - 
2 Ox 


Comparing with equation (13), and remembering that, for instantaneous 
distortion, viscosity effects are negligible, we have 


(v2 w 2?) 3aU (v2 + w?) (15) 
p\ oy oz ; 


valid for small values of a(a—zx,) and high rates of strain. This shows that 
redistribution of energy among the velocity components is an integral part 
of the distorting process, and not a secondary effect, as is sometimes implied. 
Experimental arrangements 

A duct with the velocity distribution of equations (7) is approximately 
one of constant area and rectangular section, such that the horizontal width 


— 
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Fic. 1. Arrangement of wind-tunnel. 


varies as e“"~“)), while the vertical depth varies as e-“"-*», This would be 
exact but for the growth of boundary layers on the walls, so, to prevent an 
acceleration of the flow through the cuct, the sectional area increases from 
6 in. x 24 in. at the entrance to 24-25 in. x 6-25 in. at the exit. The distorting 
duct measures 40 in. along the axis and is constructed of tin-plate with 
perspex sides. It is installed in an open return wind-tunnel, arranged as 
shown in Fig. 1. The air enters through a bell-mouth and passes through 


two 16-mesh wire gauzes and a short settling length before the 4:1 ratio. 
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two-dimensional contraction. The working section consists of a 20 in. long 


parallel flow section, the distorting duct, and a final 20 in. long parallel 


flow recovery section. The air-stream is then turned through a right angle 


by corner vanes and exhausted to the atmosphere by an eight-bladed, axial 


flow fan driven by an electric motor of } hp. 
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Fic. 2. Variations of mean speed along axis. 


The mean velocity distribution along the axis is not completely uniform, 
owing to the sudden transitions between the parallel and distorting sections 
ofthe tunnel. Fig. 2 shows that, except in these regions, the axial velocity 
is reasonably constant. No appreciable change occurs when a turbulence- 
producing grid is inserted at the beginning of the working section. 

Measurements of the turbulence were made using ordinary hot-wire 
nemometer equipment and techniques. The hot wires were supported at 
the end of a long rod running along the tunnel axis, which could be rotated 
through a right angle to compare the intensities of the two transverse 
components of the turbulent velocity or slid along to measure the variations 
of intensity along the axis 

In the absence of a grid the flow is steady, and the turbulent intensities 


the entrance to the distorting section are 


2\h/7 0-20 x 10-2. 2\h [7 0-53 « 10-2. (w?)?/U 0-52 « 10-2. 


The turbulence is composed of extremely intermittent bursts, probably 
riginating in the boundary layers on the walls, which are only 3 in. from 
he centre line. These turbulence levels increase with distance from the 
entrance. and. at the extreme limit of observation, 70-5 in. from the 


entrance. thev aré¢ 
0-37 x 10-2 2)8/ [7 0-69 « 10-2. (w?)3/U 0-71 «x 10-2, 


The increase is again probably due to growth in thickness of the boundary 


ers on the walls. These levels are rather high, but no attempt has been 
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made to correct the final readings for background turbulence. There are 
two good reasons for this. In the first place, it is known that the insertion 
of a grid modifies the wall-boundary layers, and, in the second place, the 
theoretical analysis is equally applicable to the total flow. 

The turbulence-producing grids are of the ordinary biplane type with a 
mesh diameter ratio of 5-35. Two were used, one of periodic spacing 0-5 in, 
and the other 1 in. At the mean air-speed used in the experiments (690 
em. sec.~'), the corresponding mesh Reynolds numbers were 5850 and 11700, 
Experimental results 

The variations in intensity of the components of the turbulent velocity 
are shown in Figs. 3 and 4 for both of the grids. all measurements being 
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Fic. 3. Variations of intensity (half-inch grid). 


made on the axis of the tunnel. For comparison, the intensities to be 
expected if the distorted duct were replaced by a parallel flow section are 
shown, calculated on the basis of the linear law of energy decay (Batchelor 
and Townsend (7)). It is clear that 
(a) the total turbulent intensity decreases more slowly than it would in 
the absence of the distortion, 
(b) while v? decreases less rapidly than in parallel flow, w? decreases 
more rapidly, 
(c) the rate of decrease of u2 is intermediate between that of v2 and w’, 
and is less than it would be in the absence of a distortion. 
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The relations between the intensities are conveniently studied by the 


variation of the non dimensional] quantities 
. yo ue - . - 
k, and A, (3u?/q?)—1. 
Or 
‘ 
| | 
ol | | 
| | 
} | 
‘a 
a vou, FY 
| 
—_— . 
. fe 2 
i PS 
| . 
. 
| 
. ? 2 
é 4 
a 
520. 305 601K 
Fic. 4. \ itions of intensity (one-inch grid). 
‘he variation of these quantities is shown in Fig. 5. and it appears that K, 
is not significantly different from zero. On the other hand, A, increases 
hroughout the distortion and appears to attain an asymptotic value 


The rate of energy dissipation by viscosity depends on the intensities of 


the components of the rate of strain tensor, and three of these were measured. 


In Fig. 6 the variation of A? u?/(cu/ex)? along the duct is shown. As 
indicated in the figure, in parallel flow dA2/dx 10v/U', a rate of increase 


which is significantly greater than that observed here. 


If the dissipating eddies are in a condition of local isotropy, 
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Fic. 5. Anisotropy of velocity components. 


Non-dimensional parameters similar to those used for the turbulent in- 
tensities, i.e. which are zero in isotropic turbulence, 





are given in Figs. 7 and 8. It is obvious that these eddies are not in a 
condition of local isotropy. 

Finally, a series of measurements of spectra of the longitudinal velocity 
component were made at various positions in the distortion. It is noticeable 
that the shape of the spectrum changes considerably, and that the longi- 


tudinal scale (measured as the spectral intensity at small wave numbers 
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divided by the intensity) at first increases as is usual in developing isotropic 


turbulence but later decreases. 
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Fic. 6. Turbulent dissipation parameters. 


These observations may be considered in their relation to several prob- 
ems. The first problem is the most direct, the relevance of the instantaneous 
distortion theory to the initial behaviour of isotropic turbulence suddenly 
subjected to a uniform mean rate of strain, and the nature of the departures 
irom this theory. The second is the relation of the behaviour of this flow 
to the behaviour of more complicated flows of greater practical importance. 
Finally, since laboratory studies of turbulence have been limited to flow 
Reynolds numbers insufficiently high for the satisfaction of the conditions 


r local similarity, it is interesting to inquire whether the observed depar- 


{ 


tures from local isotropy of the dissipating eddies invalidate the useful 


approximation for the local energy dissipation rate 
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Work on free turbulence indicates that the total energy dissipation can be and v 
obtained by integrating 15v(éu/éx)* over the flow, but no complete evidence lent € 
is obtainable on its local validity. behar 
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Discussion: instantaneous distortion theory 
From measurements of the intensities alone it might be concluded that 
the relevance of the instantaneous theory to the actual behaviour of the 
experimental flow is not great, even in the early stages of the distortion. 
This conclusion would be based on the continued decrease in intensity of wh 
all three velocity components, in contrast to the predicted increase of u° Th 
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be ind v?, but it is easily shown that allowances for the direct effect of turbu- 
nce lent energy dissipation are sufficient to reconcile very nearly the initial 
behaviour with theory. If a uniform rate of strain is suddenly imposed on 
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Fic. 9. gitudinal spectra (half-inch grid). 
decaying field of isotropic turbulence, then, to the first order, the changes 
the turbulence due to the uniform distortion and the normal processes 





f turbulent energy transfer and dissipation combine linearly. The reason 
| is that changes in rate of any of these processes only arise from changes in 
the turbulent structure, which to the first order are small. Provided the 


hange in the turbulence is small, for ¢ > ¢,. 
x du? f ) 
dt ; 
=  ad(v?+w?) (16) 
? - Vo Wo ———————E (¢ t,) 
1a dt 
rhe =e. 
* ye uw tal (tou )(t t,) 
of where the derivatives are taken just before the initiation of the distortion. 


The first two relations express the fact that u? and v?+-w? show no sudden 
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change in rate of decay, which is confirmed in Figs. 3 and 4. 


The last 
relation may be written 
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Fic. 10, Longitudinal spectra (one-inch grid). 


while A, is small. In Fig. 5 the variation of K, predicted by the instan- 
taneous distortion theory is shown, computed from the results given by 
Batchelor and Proudman (6). The initial slope of this curve is given by 
equation (17), and this slopé is seen to be perhaps 30-40 per cent. greater 
than the maximum slope of the experimental curves. This discrepancy is 
probably due to the impossibility of commencing the distortion instan- 
taneously, and the comparatively large values of A, which consequently 
are attained before the rate of strain becomes steady and the total strain 
known. 

Further confirmation of the relevance of the instantaneous distortion 
theory is found in the shapes of spectra of the down-stream component. 
A computation using equations (1) and a defining function 
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(which corresponds to an exponential correlation function) shows that, for 
values of the total strain ratio greater than 1-93, the one-dimensional spec- 


trum function develops a maximum at a non-zero wave-number (Fig. 11). 
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11. Theoretical spectra. 
Although the actual spectral intensities do not correspond at all closely, the 


similarity of general shape of the spectra at corresponding values of the 
total distortion is good. Not really distinct from this is the reversal of trend 
of the longitudinal scale, L,,, for large distortions, but a decrease in scale 
is very unusual in developing turbulence and deserves comment. It is con- 

uded that, although turbulent energy transfer and dissipation lead to 
considerable differences between the observed intensities and the intensities 
predicted by the instantaneous distortion theory, in the early stages of a 
iniform distortion the structure of the turbulence is substantially that to 
be expected from an instantaneous distortion. 

The significance of the approximate correspondence between the instan- 
taneous distortion theory and the initial stages of the uniform distortion of 
isotropic turbulence lies in the conclusion that the structure of turbulence 
under shear is substantially independent of eddy interaction and transport 
and the corollary that the energy-containing eddies of turbulence are 
relatively stable in structure and slow moving. Their stability is directly 
confirmed by the comparatively slow approach to isotropy of the turbulence 
alter the distortion ceases (Fig. 5) and by previous observations of the slow 
return to isotropy of the axi-symmetric turbulence obtained by passing 
initially isotropic turbulence through a fine wire gauze (Townsend (8)). 

The only previous detailed study of the effect of distortion on iso- 
tropic turbulence is that of MacPhail (4), but his results show several 
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peculiar features and do not agree with those presented here. The relevant 
measurements were obtained by passing grid turbulence through a dis. 
torted passage resembling the one described above. It changed section 
slightly more rapidly (a 2:1 distortion in 15 in.), but had only half the total 
distortion and was of smaller cross-sectional area. The turbulence was 
produced by a biplane grid of periodic spacing 2 in., placed 26 in. upstream 
of the distorting section. Over the first half of the distorting section the 
behaviour of the turbulence was qualitatively similar to that observed here. 
but. farther downstream, the anisotropy diminished and at the end of the 
distortion it was nearly zero. In the parallel recovery section the anisotropy 
increased once more. MacPhail accounted for these results by supposing 
the eddies to precess under the action of couples due to the distortion, but 
the meaning of such an effect in fluid motion is doubtful. It is more probable 
that the results are due to the large scale of the turbulence relative to the 
duct dimensions, the influence of the boundary layers on the walls, and per- 


haps the closeness of the turbulence-producing grid to the distorting duct. 


Discussion: turbulence ia equilibrium under constant rate of shear 

For total strain ratios greater than 2 the shear coefficient K, varies com- 
paratively little, and it is possible that it approaches an asymptotic value 
near 0-42. While the individual intensities may also approach asymptotic 
values for large strain ratios, it is clear that they do so more slowly, and the 
interesting point is that the turbulence appears to develop an equilibrium 
structure after a strain insufficient to establish a time invariant equilibrium 
state. This result is consistent with previous analyses of the energy balance 
in turbulent shear flows, which show that energy-production processes are 
usually too slow to account for the development of Reynolds stresses, and 
that equilibrium, if it exists, must be structural. By this is implied that the 
final configuration of the turbulence is substantially independent of energy 
supply from the mean flow and is determined by realignment of vorticity 
under the strain imposed by the mean flow. 

This establishment of an equilibrium structure for which the shear co- 
efficient is constant may be described formally as an equilibrium between 
the orienting effect of the distortion and the tendency to isotropy of the 
turbulence. This description fails to account for the very slow return to 
isotropy after the distortion ceases and the rather more reliable observations 
of a similarly slow return to isotropy of axi-symmetric turbulence. Taken 
all together, it seems certain that, once set up, Reynolds stresses disappear 


very slowly in a uniform mean flow, and it follows that the competing 
process must also be part of the distortion. It is likely that the opposition 


to complete alignment is due to rapidly increasing energy transfer from 
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‘vant eddies whose vorticity is in the direction of positive rate of strain, probably 
. dis a transfer to eddies one size smaller and containing less energy. The cause 
ction | of this accelerated transfer could be the continued crowding together of the 
tota favoured eddies with consequent increases in the velocity gradients, while 
was | eddies whose directions are not such as to gain energy will in general move 
rean apart. This description is consistent with the slower rate of increase of Aj, 
L the during distortion implying maintenance of a high rate of energy dissipation. 
here. | Complete equilibrium in which the turbulent energy remains constant 
[ the is not attained in these experiments, and very little can be said about it. 
rOpy It must satisfy the energy equation 
sing | . 
' (v= wal € (18) 
able ‘ u> L5v 
ind consequently A2 aie (19) 
the y2—_ ape al 
pel 


Under the conditions of these experiments. A? would be 0-28 cm.? Since Az 
isin the range 0-1—0-2 em.®, final equilibrium is clearly impossible. If, before 
ear entering the distortion. the turbulence were allowed to decay more, final 
equilibrium might be observed, but this is not practicable with the existing 


hue apparatus 


the Discussion: equilibrium structure in ordinary shear flows 

The comparatively rapid establishment of an equilibrium structure in 
turbulence acted on by a simple shear immediately suggests that the struc- 
ture may be similar in all fully developed turbulent shear flows, as is pos- 
tulated in the mixing-length theories. The simplest structure parameter is 


the shear coefficient, A,, where the components are along the principal direc- 


he | tions of the Reynolds stress, but a survey of published measurements shows 
2) that, while in any particular flow this parameter varies little over the regions 

of maximum shear, from one flow to another it varies between 0-40 and 0-65 

Corrsin (9): Liepmann and Laufer (10); Townsend (3); Dryden (11); 
0- Townsend (12); Laufer (13)). An examination of the values shows that, 
en generally speaking, the higher values are found in flows for which the mean 
i flow vorticity is everywhere of the same sign (boundary layers and free 
to mixing layers) and the lower values in flows of zero total vorticity (jets and 
ns wakes). It might be thought that, in the latter class of flows, the equilibrium 
on structure may never become fully developed owing to large-scale mixing 
re] between regions of positive and negative vorticity, but this is unlikely for 
ig | two reasons. In the first place, the value for the uniform rate of strain flow 
mn is much closer to the jet and wake values than to the other values, although 


mixing could hardly affect the issue, and, secondly, there is experimental 
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evidence that the large-scale mixing in a wake is small between the two 
‘halves’ though fairly vigorous inside them (Townsend (3)). 
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Fic. 12. Variations of K... 


Although the notion of general similarity must be discarded, it is interest- 
ing that a less absolute structure parameter shows a very regular behaviour, 
When describing ordinary shear flows it is usual to use a coordinate system 
with Ox in the direction of mean flow and Oy in the direction of the mean 
velocity gradient. Denoting velocities in this system by the subscript 1, a 
convenient shear coefficient is 

K,, U,V; ua, 
This coefficient is found to be reasonably constant in magnitude over the 
greater part of the flow, and variation between one flow and another is 
comparatively small (Fig. 12). In addition, the mean value is close to the 
value of K, found for uniform rate of strain. Reasons for the constancy of 
K,, are difficult to discover, but, if this result is generally valid, it is clearly 
of considerable practical value. 

The applicability of the eddy viscosity hypothesis to fully developed free 
turbulence must now be regarded as accidental, and the result of two factors, 
each of which tends to give a distribution of Reynolds stress of the type 


expected on the hypothesis. In jets and wakes there is sufficient transport 
of turbulent fluid in the central core that fully developed volumes of both 
signs of shear stress (arising from both halves of the flow) may alternate, 
more of one sign predominating farther from the centre of the flow. In the 
outer regions the effective shear stress is reduced, partly by the intermit- 
tency of the turbulence, and partly because the newly entrained fluid has 
not been distorted sufficiently to develop its equilibrium structure. 
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Discussion: structure of the dissipating, eddies 

Fig. 8 shows that, during the distortion, there is considerable and growing 
ynisotropy of that part of the motion that contributes to the intensities 
of the velocity derivatives, a result apparently in conflict with previous 
evidence that the energy-dissipating eddies of free turbulence are isotropic 
in structure (Liepmann and Laufer (10); Townsend (8), (14)). The depar- 
tures from isotropy observed here are large, especially at the lower of the 
twogrid Reynolds numbers, and are probably caused by the direct influence 
if the mean rate of strain on the dissipating eddies. If it is supposed that 
transfer of energy from larger eddies to smaller eddies is through the rate 
of strain in the larger ones (as in Heisenberg’s formulation of the transfer 
function (15)), then isotropy of the dissipating eddies demands that the 

ean rate of strain should be at least one order of magnitude less than the 
rate of strain imposed on the dissipating eddies by the rest of the turbulence. 
Since these larger eddies do not contribute appreciably to the dissipation, 
their rate of strain is at least one order of magnitude less than that of the 
dissipating eddies. It follows that the root mean square rate of strain in 


the dissipating eddies must be two orders of magnitude greater than the 


nean rate of strain if local isotropy of the dissipating eddies is to be possible. 
In the present experiments the mean rate of strain is 9-4 sec.—!, while the 
values of {(@u/éx)?|! in the distorting section are roughly 20 sec.-! and 
50 sec.-!, neither of which are even one order of magnitude greater. 
fter the distortion ceases there is, by the standards of the energy- 
mtaining eddies, a rapid approach to isotropy of the dissipating eddies. 
This behaviour is consistent with the notion that it is the direct influence 
ofthe mean rate of strain that leads to the anisotropy and not the anisotropy 
of the energy-containing eddies. The conclusion that, even at comparatively 
low Reynolds numbers, the energy-containing eddies of turbulence, though 
themselves highly anisotropic, dissipate their energy through nearly iso- 
tropic eddies one size smaller, means that the turbulent transfer process 
transfers energy from the larger scale motion to all components of the 


smaller scale motion simultaneously. A distinction between the pressure- 
velocity and triple velocity correlations in the equations of motion as the 
agents of energy interchange between the velocity components and between 
different eddy sizes respectively is not physically sound. 

That the dissipating eddies are anisotropic in uniform shear makes it very 
probable that they are also anisotropic in other shear flows of comparable 
Reynolds numbers. The evidence for isotropy in other shear flows rests on 
measurements of the intensities of the velocity derivatives where the axes 
of reference are the practical axes and not in the directions of principal rate 
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of strain. Since any anisotropy may be concealed unless the derivatives 
refer to the principal axes, this type of evidence cannot be considered con. 
clusive or even very relevant. 

Assumption of isotropy of the dissipating eddies has often been used to 
calculate the local turbulent energy dissipation from measurements of the 
intensity of a single velocity derivative. If the assumption is untrue, the 
question arises whether it is still possible to write the turbulent energy 


dissipation in the form ——~ 
= Ou, sei 
€ 1 5y (20) 
CX, 


instead of the accurate expression 


; 2 2 ‘ 2 »\2 »\2 any.\2 
(ty + (4) 4 (ce) + (=) (4) +(2) : 
CX, CY; C24 OX, CY, C2, 
/éw,\2 Cw,\? ow,\2 
ccc 
OX, CY; C24 


There is a certain amount of evidence that. in free turbulence. 


€ V 





( 4)’ (=) ( =] (: =) ( a)" 
ca —~ fon P al on — ‘Sioa “Oo |? 
CX) o OX, a OX), v4 CYUy ~ C44) 


but there is no reliable information about the remaining four terms (Liep- 
mann and Laufer (10); Townsend (14)). The validity of the short expres- 
sion is suggested by the internal consistency of energy balances drawn up by 
using it, but a more exact test would be useful. One may be obtained from 
the energy balance of this flow. Since there is negligible energy transport, 
the energy equation becomes 

ip 


aU (v?—w?)+e 0. (21) 
CX 
To compare ¢«. computed from this equation, with the values computed 
from equation (20), A? defined by 
, l5vu? 
Ne 
€ 


9 


is compared with A, in Fig. 6. The correspondence is reasonably close, and 
suggests that the approximation is valid while Ox, is in a direction of zero 
rate of strain, although local isotropy cannot be assumed. 
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APPENDIX 
[ trum function f the turbulent velocity components after distortion, 
n by equations (1), can be integrated to give the ratio of the initial and 
final tur t intensities. Batchelor and Proudman (6) have done this, expressing 
final 1 t in terms of elliptic functions, but the results quoted above may be 
1 dir . With t hange of variables, 
ka, m k(] r2)tcos 8, No k(1 —x?)tsin 6 
rrangement of the integrand, 
( Bx? (1+ Ba? 
Fes Bay ESD mt 
x 4 a—becos26 (a bcos 26)? 
i 6?+-B-2—2, 
a 1+4(1—2?)B, 
b 4(1—a?)(B?—B-?). 
. Oo! i re p j 
| Br? az?(1+- Ba? 
tJ} 1 | | — am a —* da 
Li] Ba*(l1—a?)\* {14+ Ba*(1—2?)'! 
4 B—l)a 4 B( B—1)x*— 4} Bx? 
tI} 1+ | +e — da 
: ] Bx*(1 v“) 
] | 
For B $ t integrand may be represented by a power series, and integration 
ls t 
Sarl 
(l+a35B 
3 s , 
S ng tl ul of variables, 
] «(1+ Ba? 
| kid Lk | | _—_—_— — aan dxd@ 
i J J la—beos 26° (a—bcos 26)? 


ax*(1- 
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Expansion as a power series and integration gives 
y?tw? — ——(1+8B-...). 


In the same way 


*—B 





1 27 
2 [ | k v+4B (14+ Ba?)(1+3B2?) 


52 __ ap? Adi(ke) dk sat -| dad 
i, | k*p(k) tB t— bcos 20 (a—b cos 26)? — 
0 


=—3 9 


1 
8x1 [ 1—2z?+43B 


BB 3 


The leading terms of the series are then 








a (1+ Bat(1—2?)}t” (1+ BaX(1—2*)}! 


a(1+ Ba?)(1- A) hh 
7 ax. 


- i 
us,,(0) | | F(0,m,n) dmdn = —wL,, 


7 
co —o 


and y,,(0) is easily shown to be invariant during instantaneous distortion. With the 


change of variables , 
. —- m, = keos8, Ny = ksin8 


ad , 1+(B-4— 1)cos*é + (B4— 1)sin40 
ub, (¢ bes diec 
ul )) | | (ih )k (B 2 cos B?sin?0)? lk 10 





0 0 


XL 


= 27 | d(k)k dk, 
0 
which is independent of f. 
The longitudinal scale of the transverse component, v, is similarly defined by 


oO x 


Ye oo 
us,,(0) | | G(0,m,n) dmdn = =v? L,. 


Using the same change of variables, and after some rearrangement of the integrand, 


rr 48? sin? 
.(O d(k)k? ——_—__— — dad 
Pol ) | | p(k )k [B? |. 8-2 (B? B 2)eos 26]? ned 





0 0 
a | d(k)k® dk, 
0 
which is also independent of 8. The proof that y%,,(0) is independent of £ is exactly 
similar. 
The one-dimensional spectrum function for wu is defined by 


x x 


(1) - f [ F(l,m,n) dmdn. 


=—a = @ 
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Changing the variables to k, @ defined by 


Mo k?—[2)teos 6, Ng (k?—[?)tsin 0, 





; — Br 12(k? +- Bl*) ’ 
(1) | | keg(k [1+ B’c0s 20 ~ (a 8 cos apy) 249, 
wher a’ k2+4B(k?—P), 


b }(p2 B-2)(k2 12), 


after some rearrangement. Integrating with respect to 0, 





‘ Bi? I2(k?-+ oa | 
] Dae kedtk | cea _ AAT rans ond | . 
| Peon} Bk? BlyA (+ BPR? BE) - 








' 
Assuming that the defining function is 
2u2 p 
$(k) = ——F 
| ar(k? +- p?)' 
! corresponding to an exponential correlation function), and substituting 
ou: PAS 
it is found that 
| a 1 
Lay =a L+28/p* 1+3B(1—2)+}3B*x(1—2) dx 
wb ——se + 8, | 2 BR 7 2/p2)3\" 
aw (1+/?/p?)? p* | {1+ Ba(1—2z)}! (a+ 1?/p?) 
0 
This integral has been computed for several values of B with the results shown in 


Fig. 11. 





CORRIGENDUM 
THE STABILITY OF A COMPRESSED ELASTIC RING 
AND OF A FLEXIBLE HEAVY STRUCTURE SPREAD 
BY A SYSTEM OF ELASTIC RINGS* 
By G. W. H. STEVENS (Royal Aircraft Establishment, Farnborough) 
[Received 17 April 1953] 


Ix a private communication to the author, Arthur Boresi, University of 
Illinois, has shown that the criterion for the stability when the loads act 
towards the centre is incorrect in the above paper due to an error of 
algebraic manipulation. Equation (20.2) should read 


(k2—1) El r Too 


P—2 Be. 


{0 


The statement preceding this equation, referring to equation (12), is correct 
and the sign convention of the paper remains consistent. 

This slip has little bearing on the discussion of the hooped skirt because 
the assumption was made that the loads act normally. The comment 
following the table in section 5.2 of the paper is, however, incorrect because 
the value of q,73/HJ at each transition point is unaffected by the 
direction of loading. The revised criterion shows that if tangential forces 
were significant in the example of the hooped skirt then the transition 
ratios r,/r, would all be lower. 


r G. W. H. Stevens, QVuart. J. Mech. and Applied Math. 5 (1952), 221 -36. 


{Quart. Journ. Mech. and Applied Math., Vol. VII, Pt. 1 (1954)] 
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